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ABSTRACT 

This project involved the design and evaluation of a 
radioisotope on-stream coal ash analyser for use in fine 
coal preparation circuits. The instrument was designed to 
predict the solids concentration and ash content of a slurry 
fFlowstream while accommodating changes in the mineral matter 
composition, without recalibration. 

The work included the development of a methodology for 
establishing the economic merit and technical feasibility of 
regulatory control systems, under the umbrella of process 
optimization. Some plant experimental work was conducted to 
justify the project and to provide the parameters necessary 
hor both the equipment and the experimental designs. 

A review of the literature revealed that a limited 
amount of instrumentation development had occurred in this 
area, mostly in monitoring the ash content of coarse dry 
coal. A number of candidate designs were considered and the 
one studied by Kawatra on a bench scale was adopted for 
further investigation on the prototype scale. The theory of 
instrument operation was developed in some detail to provide 
both design information as well as semi-empirical 
calibration models. 

A prototype instrument was constructed and calibration 
Studies performed in the laboratory. The experimental 
results on simulated froth flotation product flowstreams met 
the evalutaion criteria, with one exception. The latter was 


the direct result of a design problem and should be easily 
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corrected in the next generation of instrument development. 
The laboratory results were sufficiently encouraging to 
warrant a plant trial of the prototype. Although 
' experimental problems partially masked the true potential of 
the device, it was still possible to obtain a successful 
instrument calibration. This exercise served to demonstrate 
that an on-stream ash analyser of the type studied could be 
used in process monitoring and in the process control of 


fine coal preparation circuits. 
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1. INTRODUCTION 


1.1 INDUSTRY OVERVIEW 

Coal ee (coal washing, coal cleaning, coal 
beneficiation) is the art and science of removing inert 
inorganic material from the run-of-mine raw coal. The 
principal reasons for undertaking this step in the mining 
sequence are: 
a) to comply with market specifications on product quality, 
b) to minimize the cost of transporting the coal to distant 

markets. 
According to statistics compiled in January, 1981(1), 
Canadian coal production was estimated at 36.5 million 
tonnes per year. Of this, 19.5 million tonnes (53%) were 
prepared to some extent prior to shipment. This coal is 
destined primarily for international and distant domestic 
metallurgical and thermal coal markets. Recent 
forecasts(2,3,4) predict a substantial increase in coal 
production over the next two decades. Although the current 
economic downturn has slowed growth in both areas of 
utilization, significant growth is still expected to occur, 
even in the short term(4). In view of these figures and the 
expected trends, coal production and specifically coal 
preparation plays an important role in the Country's 
economy, particularly in the international trade balance. 

Most of the Canadian coal preparation plants are 


relatively new having been commissioned in the period since 
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1970. However, since this represents a re-emergence of coal 
after an extended period of stagnation, the industry lags 
the other process industries in terms of design and 
operating poset ieee as well as in research and development. 
While much of the newer design and operating philosophy has 
been borrowed from overseas experience (eg. UK- National 
Coal Board) this is not a panacea. Canadian coals possess 
some unique preparation problems, in particular the 
requirement to treat large amounts of fine coal, for which 
borrowed technology is largely inapplicable. 

In the process industries, understanding process 
mechanism (through modelling) and installing regulatory 
control systems are ‘pre-requisites to optimal ‘plant 
operation. The coal producers recognize this as an important 
area and are concerned with improving upon existing methods, 


especially in the treatment of fine coal. 


1.2 AN INTRODUCTION TO COAL PREPARATION | 

It is assumed that the reader may not be familiar with 
coal preparation. In order that he may more easily 
understand some of the material which follows, a cursory 
treatment of the subject is provided. 

The overall objective of coal preparation is to treat 
the material from the mine in such a manner that the organic 
material is effectively separated from the extraneous | 
mineral matter producing a maximum amount of clean coal’ 


‘It is implicit that maximum production ensures maximum 
protit. 
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which complies with the sales contract specifications. 

It iS necessary to define some of the jargon: 
COAL: A naturally occurring, rocklike derivative of forest 
type vegetation which contains 250% by weight of organic 
material. 
MINERAL MATTER: The inorganic (usually solid) portion of the 
coal which consists primarily of silaceous material. 
Material which is colloidally dispersed in the coal organics 
matrix, aS a result of the depositional environment, and 
which cannot be separated by physical means, is called 
Inherent Mineral Matter. Material which is loosely bound to 
the organic material arising from dilution during mining or 
interbedded seams, and which can effectively be removed by 
physical means, is called Extraneous Mineral Matter. 
ASH: The residue which remains after complete combustion of 
the organic portion of the coal. Since there is some 
volatile matter associated with the mineral matter (eg. 
water of hydration), the ash fraction of the coal is always 
less than the mineral matter fraction. However, fora 
particular deposit the Mineral Matter to Ash Ratio is 
approximately constant (=1.1). 
YIELD: The ratio of clean coal produced to the feed ona 
weight basis for a particular separation (eg. unit, circuit, 
plant). | 
CLEAN COAL: The organic rich product of a separation. 
REFUSE(REJECT): The mineral matter rich product of a 


separation. 
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Historically, coal cleaning was based on exploiting the 
difference between the specific gravity of the organic 
fraction (typically 1.25 < sg < 1.45) and the mineral matter 
fraction (typically 2.3 < sg < 2.8). (It must be recognized 
that as a result of the incomplete liberation of the 
extraneous mineral matter there is a continuous distribution 
of solid material with specific gravity.) Unit operations 
Such as the jig, table and, more recently, dense medium 
vessels and dense medium cyclones were used to effect 
separation. In general, separation efficiency decreases with 
decreasing particle size. The inability of these devices to 
effectively treat fine coal coupled with increased need to 
do so has necessitated the development or application of new 
process technology. This includes devices which will 
Separate on the basis of specific gravity (eg. automedium 
cyclone) or on the basis of surface properties (eg. froth 
flotation) 

To obtain maximum separation efficiency in a plant, the 
feed material is usually classified on the basis of size, 
into two or three streams. Typically a coarse cut is made at 
12.5mm and a fine cut at 0.6mm, resulting in three separate 
feed streams. This design allows optimal equipment selection 
with respect to maximizing separation efficiency, assuming 
that economic justification exists. A general flowsheet of a 
Standard western Canadian coal preparation plant is included 


in figure 1.1 to illustrate these design concepts. 
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Figure 1.1 A Typical Western Canadian Coal Preparation Plant 


Flowsheet. 
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1.3 PROJECT HISTORY 

Western Canadian coals are very friable, particularly 
the metallurgical coals in the Rocky Mountain region. This 
results in a very high percentage (typically 225%) of fine 
material, nominally <0.6mm, in the preparation plant feed. 
In this respect these plants are unique for they must wash 
very substantial tonnages of this fine coal in an efficient 
manner, uSing processes which are relatively new to the 
field. Unlike coarse coal washing, this task is quite 
difficult because of the general lack of process 
understanding and the absence of the necessary process 
control equipment, in particular, on-stream analytical 
devices. 

The beneficiation of fine coal is an area which 
requires a considerable amount of research and development 
work. Some of this work has been undertaken at the plant 
level but this is usually in response to an operating crisis 
such as very poor separation efficiency. This leaves the 
majority of the more fundamental research to be performed by 
publicly funded research organizations and the universities. 
In recognition of this need for research in fine coal 
preparation, L.R. Plitt? conceived a "Coal Beneficiation 
Project." In response to his proposal the Department of 
Energy and Natural Resources of the Province of Alberta 
provided a capital grant of $250,000 in support of this 
work. It was envisaged that the work would proceed on two 


2 Chairman, Department of Mineral Engineering, University of 
Alberta. 
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fronts, namely: 

(1) laboratory and plant analysis of the process unit 
operations (eg. froth flotation, automedium 
cyclones, classifying cyclones, screens, etc...) 

2:) instrumentation and control of fine coal separation 
processes. 

Within the objectives of the project it was deemed that the 

development of an On-Stream Coal Ash Analyser (OSCAA) should 

be given the highest priority. The reasoning was that the 
successful development of such a device would provide the 
basic tool for process control. Furthermore, it would 
facilitate process studies in the plant by eliminating what 
would otherwise be a horrendous sampling and analytical 
overhead. 

The OSCAA project was the subject of this thesis 
investigation. It officially began in September, 1978. 
Funding for the operating expenses was obtained from the 
Alberta/Canada Energy Resources Research Fund (A/CERRF) at 
the following levels: 

YEAR 1 $33,000 
YEAR 2 $40,000 

The goals of the OSCAA project were: (1) to develop a 

prototype instrument in the form of a laboratory pilot plant 

system, (2) to calibrate (or otherwise study) the equipment 
performance in the laboratory, and, (3) to test the 


prototype in the operating plant environment. 
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Under the conditions defined by the Alberta/Canada 
Energy Resources Research Fund (A/CERRF), an Alberta based 
producing company was to be chosen to participate in this 
project. McIntyre Mines Ltd., Grande Cache, Alberta, was 
selected because of their expressed interest in the work as 
well as the strong personal contacts between this group and 


the university personnel involved. 
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2. OPTIMIZATION IN COAL PREPARATION 


2.1 INTRODUCTION 

Before proceeding with the experimental portion of the 
thesis it was necessary to justify the work. While the 
project was intuitively justified, experimental evidence was 
required to allow for an independent and unbiased assessment 
to be made. The problem was that there existed a void in the 
technology, namely a methodology for evaluating the impact 
of process optimization. In short, a description of what 
measurements had to be made and how they were to be 
analysed. 

Process optimization requires that the various circuits 
within a plant be manipulated in some organized pattern such 
that the optimization objective is ultimately met. Implicit 
in this manipulation is process stability. Since stability 
is rarely inherent some form of process control is 
frequently required. In such cases process sensors are 
necessary to allow control strategy development and 
implementation as well as in the subsequent process 
manipulation. In this context, optimization includes two 
broad areas; the short term - process stabilization, and the 
long term - process manipulation. 

The material which follows represents an effort to 
develop a methodology to assist in process analysis as well 
as in the a-priori justification of process instrumentation. 


With respect to the longer term activity the steady state 
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optimization criteria are defined and their practicality 
discussed. For the shorter term activity, methods of 
assessing the economic merit and technical feasibility of 
process control systems are presented. This provides the 
theory which the operator requires to make an independent, 
and not intuitive, evaluation of the plant performance. 

The author recognizes that this material is not 
Strictly a part of the development of the OSCAA prototype 
instrument. However, it is his opinion that such a treatise 
was absolutely essential as the foundation for project 
justification. Moreover, this chapter by itself represents a 
scientificrcontribution ito the ‘industry “and, asa 


consequence, should be included in the thesis. 


2.2 OPTIMIZATION: OBJECTIVES AND OBSTACLES 

Once a washery has been commissioned and most of the 
initial mechanical bottlenecks have been removed, it remains 
to optimize the preparation plant system. Optimization in 
coal preparation is simply making the maximum economic use 
of the available coal reserves. This translates to 
maximizing clean coal production from the reserves subject 
to the constraints (quality and quantity) of the sales 
contract(s) and the external regulatory processes. It would 
be desirable to have the flexibility of rapid tuning of the 
plant to take advantage of spot markets as they present 


themselves. 
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The ideal situation for the coal preparation plant 
engineer is to have a totally homogeneous plant feed 
(composition, size, and density distributions) and a well 
controlled plant. Under these conditions the individual 
cleaning circuits within the plant may be tuned (e.g. 
EVOP(5), SSDEVOP(6)) to provide for maximum profitability 
over the lifetime of the sales contract. However, because of 
the varying nature of the coal deposit quality, the 
reliability characteristics of the mining system, and the 
lack of process control equipment in the plant, this is very 
seldom possible. More typically the engineer is faced with 
defining optimum operating policy on a daily basis given: 

a) the upstream feed disturbances caused by variations in 
seam quality and mine production, 

b) poor mining methods, 

c) the downstream constraints of pie product 
Specifications, 

d) recent production history, 

e) his knowledge of process behaviour. 

Considering the stochastic nature of (a), the problems 

associated with obtaining a reliable material balance in (d) 

and the general lack of process understanding in the 

industry, this is a difficult task. Another serious problem 

arises in (b). Since coal preparation unit costs are often 

significantly lower than the mining unit costs, it is not 

uncommon for the mine to off-load much of its "selection" 


reponsibilities to the plant. The problem is that coal which 
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is expensive to mine usually requires a high washery 
efficiency. This serves to further complicate the quality 


eontrol -function, 


2.3 THE THEORETICAL STEADY STATE OPTIMIZATION STRATEGY 
In principle, it is easy to define the steady state 
optimization objective for a coal preparation plant. 
"Maximize the yield of clean coal from the available 
reserves while maintaining the the specified product 
quality." 
Mathematically, 
Maximize: 
ae es sj (ra: ks) 


SUbRECE LO? 


De De BE a8 by x. alee?) 
atk 


where: Y =objective function (plant yield). 
M,;,; =the fraction of raw coal from block j (which 
has homogeneous properties) in the deposit which 
reports to cleaning circuit i in the plant with 


ZIM ; ;=1. 

Yi =the fractional yield of M,,; to the clean coal 
product. 

A; ,; =the value of some cumulative attribute in the 
Clean. coal product’ (e.g> ash, sulfur,’ calorific 
value). 


Ao =the quality specification for the attribute. 
The attribute of greatest interest in Canadian operations is 
the ash content. Sulfur is of lesser importance because of 
its relatively low concentration. Furthermore, for most 
deposits, calorific value can be inferred from a functional 


(usually linear) relationship with the ash content. 
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Therefore, ash content is generally the attribute with which 
the operator's characterize the various solids flowstreams. 
For this reason equations (2.1) and (2.2) will be solved for 
the ash constraint. However, any number of constraints could 
be included one of which would determine the extent to which 
plant yield may be maximized. 

For any attribute, such as cumulative ash, which 
increases monotonically with increasing yield, maximum yield 
will be obtained when the constraint of equation (2.2) is a 
Strict equality. By nature, the physical and physiochemical 
Separation operations employed in beneficiation exhibit 
operating yield-ash curves which are monotonically 
ineLeasing.in theyregiuonyof interest. Inva properly operated 
plant the overall yield-ash curve will also be monotonically 
increasing in the region of interest. An example of a 
yield-ash curve for a simulated three circuit thermal coal 
preparation plant is presented in Figure 2.1 to illustrate 
this monotonicity. It is possible that some particular 
combination of these units will produce an overall plant 
product curve which does not meet this condition (e.g. see 
(7)-Figure 18-18). This is the result of unusual operation 
and should be infrequently encountered in practice. 
Therefore, in normal operation, the plant yield-ash curve is 
such that the strict equality of equation (2.2) will ensure 
maximum yield. 

The system of equations (2.1) and (2.2) can be solved 


analytically using the method of Lagrangian Multipliers(8). 
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Figure 2.1 The Yield-Ash Curve for a Simulated Thermal Coal 


Preparation Plant. 
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Incorporating the constraint into the objective function 
gives the Lagrangian: 


Maximize: 


aiag pie? 7B is es mi x th Maa =, p> Meer aa 
where: <A =the nisi Multiplier. 

Every circuit has some control variable(s) which can be 
altered to produce a change in Y;; (A;;). Equation (2.3) can 
be maximized with respect to the Y,, which then defines the 
control variable settings and A;;. Taking the general 
partial derivative, setting the result equal to zero and 
rearranging gives: 


eee iy. = (1+AAQ) 


Vaud r 
| 


pane 24) 


The interesting result is that the right hand side of 

equation (2.4) is a constant, regardless of the coal block 
or cleaning circuit being considered. Now the ash content, 
a', of an incremental yield element AY of coal is defined 


ass 


dak. -5) 


where: Y =yield. 

=ash. 

' =the ash content of the incremental yield 
element, hereinafter referred to as incremental 
ash. 


o > 


In the limit as) AYs0*, equation (2.5) reduces to the 


differential equation: 
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a =6 (YA) ea t6) 


From equations (2.4) and (2.6) the following important 
conclusion can be drawn; 

To maximize production the coal must be washed to a 

constant incremental ash level regardless of the 

source or separation method. The actual value of 

this incremental ash is defined by Aj» and by the 

circuit yield-ash curves. 
This argument, albeit restricted to perfect separators’, was 
advanced by Bird and Marshall in 1931(9). The result is 
intuitive if, for example, one considers the specific case 
of a two circuit plant. If these circuits are not operating 
at the same’ incremental ash level then it is theoretically 
possible to increase the plant yield while maintaining the 
Same product ash quality. If circuit 1 has the lower 
incremental ash then suppose an increment in production is 
obtained from this circuit with an equivalent decrement in 
Circuit 2. The result is that the overall yield will remain 
unchanged but the product ash will decrease as a result of 
the quality differences in the production elements. When 
there is no difference in the incremental ash levels there 
is no further driving force for change. This difference is 
exploited until the yield is maximized at the product ash 
Specification. 


7In perfect separators the incremental ash is equivalent to 
the "elementary ash", a term which is frequently encountered 
in washability studies. 
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To provide an appreciation of the importance of abiding 
by the steady state optimization criteria, Figure 2.2 is 
presented. This figure shows yield isopleths at a constant 
- plant product ash quality for the simulated thermal coal 
plant mentioned in connection with Figure 2.1. (The axes 
represent the control variables for the two dense medium 
circuits employed in the simulation, i.e. the specific 
gravity of separation.) As expected, the optimal (maximum) 
yield corresponds to the solution of equation (2.4) for the 
particular value of A, chosen (11.5%). While the response 
surface is relatively flat in the neighbourhood of the 
optimum it is apparent that yield decrements on the order of 
~1% could be easily incurred with slightly off-optimum 
operation. Even with thermal coal, losses of this nature 


represent a significant economic penalty. 


2.4 CONSTRAINTS ON STEADY STATE OPTIMIZATION 

Having defined the optimum operating scenario it is 
important to look at the technical practicality of 
implementation. It would be very difficult to measure 
incremental ash on-line. Probably the easiest method of 
implementing this strategy would tavolve model inference. 
This requires an accurate knowledge of the washability 
characteristics of the coal as well as the separation 
characteristics (partition curves) of the various unit 
operations in the flowsheet. (A further requirement is 


Stable operation through close regulation of the process 
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Figure 2.2 Yield Isopleths for a Given Product Ash in the 


Simulated Thermal Coal Plant. 
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control variables.) Implementing such a technique would 
require a very extensive investigative effort. 

It has been argued that for density separators, which 
have near ideal partition curves, the incremental ash is 
strongly correlated with specific gravity. Therefore, fixing 
the specific gravity of separation will fix the incremental 
ash. Unfortunately, the ash-gravity correlation is often not 
that good, and the partition curves may deviate from near 
ideality depending upon changes in the operating conditions. 
Specifically, rigorous control of the separation density 
cannot be guaranteed under manual control. In Summary, 
although the plant steady state optimization criteria are 
known, it 1s seldom possible to conform to them unless both 
the plant feed and plant operation remain stable over fairly 
long periods of time. 

This raises another practical constraint on the 
theoretical steady state optimization strategy. That is, 
while the overall ash specification will be met in the long 
term, in the short term plant product quality may deviate 
considerably from Aj. This depends upon the washability 
characteristics of the particular coal block being treated 
at the time. This kind of deviation may not be tolerable and 
therefore it would be impractical to optimize production 
over the entire deposit. There is a time constraint (period) 
over which this methodology should be implemented, which is 
considerably shorter than the lifetime of the sales 


contract. The reason is that the buyer is expecting a 
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consistent product for coke blending or steam raising 
purposes. That is, each shipload of coal must meet contract 
Specifications or be subject to penalty. For a typical cargo 
ship with a payload of 100,000 tonnes and a "typical" coal 
plant producing 2.5 million tonnes of clean coal per year, 
this time period is about 15 days. Furthermore, there exists 
an additional time constraint because of both unexpected 
upstream variations in the plant feed* as well as the 
methods used for storing coal at the terminal. In the latter 
case, stacking will result in additional product blending 
which may be deleterious to the anticipated quality of the 
shipment. The result is that the true optimization period is 
the time required to fill one unit train (capacity ~10,000 
tonnes) or approximately 1.5 days. In view of this time 
scale, and while longer term optimization remains important, 
the emphasis shifts to the shorter term optimization 
activity. 

Considering this time constraint and the practical 
difficulty of implementing the constant incremental ash 
control philosophy, the operators attempt to maintain a 
constant plant clean coal product ash, A,. (This target 
value is frequently slightly below the specification, Ao, to 


leave room for error. This is costly in terms of lost 


production.) For multi-circuit plants’, the operators set 


* If, for some reason, the mining plan is altered, then the 
optimization exercise will be invalid if the coal to be 
cleaned has different washability properties from that which 
was expected. 

> Most of the newer Canadian coal preparation plants feature 
multi-circuit flowsheets. 
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individual circuit ash control targets which, based on their 
knowledge of the plant feed characteristics, will provide 
the desired quality of clean coal at _a maximum yield for the 


' daily production. 


2.5 EXTERNAL OPTIMIZATION ASSISTANCE 
There are a variety of means by which the coal 
preparation quality control function can be made easier, 
both from within and outside of the plant. Considering the 
latter case, there are three principal methods of reducing 
the "disturbance load" on quality control in the washery. 
a) Blending and Feed Forward Control Through Mine 
Production Scheduling. 
Provided that a reasonable model (preferrably from 
geostatistical analysis) of the deposit exists, the mine 
is capable of providing quality information to the 
preparation plant in advance of raw coal delivery. This 
information may be used to set the plant for optimal 
cleaning of a particular blend‘. 
b) Pre-Process Homogenization’ 
The use of large raw coal feed homogenization operations 
(e.g. stacker and reclaimer) prior to the preparation 


‘Some interesting work has been done in France where the the 
parameter set which describes the washability of the coal 
has been modelled using a geostatistical approach. In this 
case the mine could supply the plant with the washability 
curves for the coal they will be receiving which would 
simplify the control problem. (see J. Int. Ass. Math. Geol., 
Vol « =13', “NOP "4 F196 FSG. 93213) 

7These methods are much more rigorous in their attempt to 
blend the coal than the standard product surge stockpiles 


' common to all coal preparation plants. 
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Stage will yield a more consistent plant feed. A more 
consistent feed reduces the demand for process control 
action and will improve yield. 
‘¢) Post-Process Homogenization 
This method involves the homogenization of plant 
product, as described by Butcher and Beninger(10). In 
essence, the implications on preparation pant operation 
are identical to those in (b), however, the total solids 
to be handled is significantly reduced in comparison. 
Of the three methods, only (a) iS practiced to any 
extent in North America. Method (b) is popular in Europe and 
method (c) has been investigated by various operating 


concerns but has never been implemented on a large scale. 


2.6 PROBLEMS WITH CURRENT PROCESS CONTROL PRACTICE 

It was stated above that the typical operator process 
control strategy was to establish target clean coal product 
ash values for, the various cleaning-circuits. Using "grab 
samples" with turn-around times of from 2 h to 24h, 
depending upon the plant, the control variables are modified 
in an effort to maintain the product ash close to set point. 
In many instances this feedback control loop is ineffective 
for the following reasons: 
a) sample frequency is too low, 
b) sample is not representative, 
ea, ine human controller does not have a quantitative 


understanding of the process, i.e. the relationship 
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between control variables and process behaviour. 
In these cases only long term (i.e. daily or weekly) control 
action is specified depending upon trends in the the 
cumulative plant product quality. 

The task of maintaining clean coal product quality 
while maximizing yield from a preparation plant is not easy 
under present conditions. There is strong motivation to 


improve upon current methods. 


2.7 OPTIMIZATION THROUGH PROCESS CONTROL: ECONOMIC 
IMPLICATIONS 
Choosing a target ash value in each cleaning circuit to 
achieve optimum yield is, in essence, the operators's 
physical solution to equations (2.1) and (2.2), under the 
constraints previously discussed. However, even though 
control set points are established, it can be easily 
demonstrated in most plants that there is a significant 
variation in product ash. Furthermore, the mean values may 
not correspond to the targets. These variations are the 
combined effects of process disturbances and ineffective 
process control. The disturbances may be classified as: 
(te) Intrinsic Disturbances 
Those disturbances resulting from changes in the 
physical and/or physiochemical properties of the 
coal (e.g. washability characteristics, size 
consistipntloatabidi ty) etc. is) 


(20) Extrinsic Disturbances 
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Those disturbances resulting from changes in the 
normal process equipment operation, the control 
variable settings and the "gross" character of the 
feed stream (e.g. mass flowrate, slurry solids 
content, media density, etc...). 
Changes which arise from intrinsic disturbances require that 
the target value for the circuit be re-assessed. This 
requires detailed knowledge of both the new properties as 
well as of process behaviour. It is important that the 
operator identify and react to changes of this type, 
however, this is not an easy task. In fact, solutions to 
control problems of this type can only be evaluated when 
successful strategies and systems have been developed for 
extrinsic disturbances. 

Extrinsic disturbances are important since, in these 
cases, deviations from set point are undesirable. In 
general, each washing circuit within a coal preparation 
plant is operated independently, reflecting their physical 
independence. (e.g. see Figure 1.1) Each circuit should be 
controlled to maintain target for the sake of operational 
Simplicity. The discussion which follows relates to the 
economic implications of extrinsic disturbances under this 
set of circumstances. 

In the discussion, it is implicity assumed that the 
target ash set points for the individual circuits are set to 
provide equivalent incremental ash in the clean coal 


products. Clearly, any deviation from this condition will 
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result in a production loss. Because of the complexity of 
the control problem it is expected that this production loss 
is beyond the limits of the control systems currently 

’ envisaged for coal preparation plants. For this reason this 
type of production loss is excluded from further 
consideration. 

One can consider two contributions to production loss 
which result from variation in product ash due to extrinsic 
disturbances. The first arises simply because of the 
variation itself (direct production loss) and the second 
arises from conservative set point specification (indirect 
production loss) as a result of the variation. They are 


treated separately below. 


2.7.1 Direct Production Loss 

There are a number of mathematical means of 
demonstrating that production from a process, operating on a 
given yield-ash curve, is maximized if the process is 
Maintained at the set point. To minimize the number of 
restrictive assumptions which are required for analytical 
solution and to demonstrate the interrelationship of the 
variables, a simple technique is presented. 

Rather than operate exclusively at set point, coal 
cleaning units, or circuits, or plants, show a distribution 
of operating regimes about the overall mean. This can be 
seen by sampling a product stream over time and analysing 


for the solids ash content. The distribution can be thought 
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of as balanced periods of undercleaning and overcleaning 
such that the overall mean ash content is maintained. This 
is the approach used by Cammack and Balint(11). Although 
arriving at the correct conclusion their algebraic proof was 
unsatisfactory, in the opinion of the author. An alternate 
solution is offered here. 

The ash balance for the off-average operating periods 


is written as: 


= bint YAAP te +h Act 
A= ok ee EY a 3 


poet ey et Mee aie) 


where: (A,Y) =mean operating point on the yield-ash curve. 
(A,,Y,) =overcleaning (i.e. A,<A) operating point. 
(Az, Yo) sundercleaning (1.6. A->A) operating point. 
t.,t, =time of undercleaning and overcleaning 
respectively. 
F =solids mass flowrate to the cleaning process. 


Assuming that in the neighbourhood of (A,Y), the yield-ash 


curve can be approximated by a linear function of the form: 


Y=Y+m (A-A) a DNs ae 
where: 
m=dY 
dA (A,¥) bee zr) 


the tonnage loss, L, over the combined period (t,+t,) is 


calculated as: 
Peariv(tette ev bev ee } Nae Decl 
pel iy. Daeg ay Ae 


Substitution of the appropriate form of equation (2.8) into 


equations (2.7) and (2.10) with rearrangement allows 
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equation (2.10) to be written in the form: 


L=Fm(A-A,) t,, [2- Er] eet) 
Y+m(A_—A) 

This equation constrains L to be non-negative. The 

production loss is observed to increase with’; 

a) an increase in the slope, m, of the yield-ash curve at 
the mean operating point, 

b) the magnitude of the differences between the off-average 
ash values and the mean, 

c) the duration of the off-average operation. 

On the basis of this analysis the process control or 
optimization objective is to minimize the effects of (b) and 
(c) on process system performance. 

To translate these observations into the production 
environment and to illustrate the effects of (b) in 
particular, consider the curves shown in Figures 2.3(A) and 
2.3(B). This simplified example indicates the impact of 
variance reduction on direct production loss. Case 1 on 
Figure 2.3(A) shows a separation with a highly varying 
instantaneous product ash. Combining this instantaneous ash 
curve with the yield ash curve of Figure 2.3(B) allows the 
cumulative production statistics to be calculated, (A,Y). 
Observe that the cumulative product ash (as t7#) has a value 
of A,=9.21% which is greater than the mean value of the 


instantaneous ash, €@,=9.0%. This is expected since the slope 


‘These findings have been verified using numerical methods. 
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Figure 2.3 An Illustration of Direct Production Loss Through 


Quality Variation. 
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of the yield-ash curve at the mean operating point is 
non-zero. Now, assuming that the variance in the 
instantaneous product ash can be reduced by ~94% gives the 
results shown for ase 28 Here, the mean instantanteous 
product ash value, 42, has been adjusted to make the 
cumulative ash values equal, A,=A,. The important 
observation is that the yield increase obtained by reducing 
the variance is significant (2.3% for the example). 
Moreover, as the variance is reduced and, aS a consequence, 
A7xa*, the potential for improving the yield decreases very 
rapidly. However, experience has shown that many western 
Canadian coal processing circuits have a rather substantial 
variance in the instantaneous product ash and there is 
considerable economic incentive to pursue variance reduction 
through improved process control. An example of this will be 


given in sub-section 2.7.3.2. 


2.7.2 Indirect Production Loss 

Producers generally specify the plant product ash 
value, A,, to be somewhat lower than the sales contract 
value Ao. The magnitude of the difference, A,-A,, 1s 
governed by the variance of A, in the shipments. The 
objective is to reduce the probability of delivering coal 
with an ash vane greater than Ag. Methods of reducing this 
variance such as those described above will allow A, to be 
increased while maintaining the low probability of over 


specification shipments. This is illustrated in Figure 2.4 
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where, for the sake of simplicity, the ash quality of the 
clean coal shipments has been assumed to be normally 
distributed. (Wells (see (12) Chap. 3) has used this in 

~- analysing the economic payback of control systems.) 
Increasing A, will result in an increase in plant yield, Y,, 


and hence reduce the indirect production losses. 


2.7.3 Methods of Assessing the Economic Impact of Improved 
Process Control 
Before committing the required capital and operating 
funds to a process control project, management has to be 
convinced that; 
eye athe technology, todo the job either exists or can be 
developed with the resources at hand, 

b) the project makes good economic sense (e.g. acceptable 
discounted cash flow rate of return). 

This sub-section deals with some facets of (b). 

In order to perform the economic analysis it is 
necessary to estimate the direct and indirect production 
changes which are expected from the reduction in product ash 
variance. This knowledge also allows the operator to 
establish priorities within the plant. There are two methods 


which may be used to quantify the production change. 


2.7.3.1 Direct Estimation of Yield Increase 
Based on the experience of other plants which have 
undertaken similar projects, or, on good engineering 


judgement, one estimates the yield increase expected 
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Figure 2.4 An Illustration of Indirect Production Loss. 
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from improved process control. The additional production 
can be calculated from the feed rate to the process 
System, and, after correction for the increased 
operating and capital costs, the before tax cash flow is 
estimated. In the absence of any other information 
Wells(12) suggests that a yield increase of ~1% is a 


good conservative estimate. 


2.7.3.2 Estimation of Yield Increase Through Variance 
Reduction 

This method represents a more realistic means of 
estimating the yield increase discussed above. It is a 
generalization of the concepts introduced in the example 
presented in Figure 2.3. (The "best" technique for 
quantifying the reduction is treated in detail in the 
next section.) It was been observed in connection with 
equation (2.11), point (b), that variance reduction has 
the overall effect of reducing direct production loss 
and allowing for the reduction of indirect production 
losses. The former is generally much easier to calculate 
than the latter. The reason is that the contribution of 
the product ash variation of a circuit to the overall 
plant product ash variation is a function of the period 
required to produce a shipment. That is, it is required 
that the variation in ash quality per shipment period 
from the circuit be known. This depends upon the 
frequency of the disturbances which enter the circuit 


and their effects on product quality. This kind of 
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information is only available through extensive 


Sampling. It is suggested that economic calculations be 
based on direct losses, recognizing that there will be 
some further benefit which may accrue from indirect loss 


reduction. 

To estimate direct production change from variance 
reduction requires the following: 

a) A knowledge of the typical distribution of product 
ash quality in time. This is usually obtained by 
high frequency sampling (i.e. >3 samples per h)over 
an extended period (>100 samples). 

b) A knowledge of the yield-ash curve over the region 
of interest. This can be obtained through 
experimentation, by mathematical models and digital 
Simulation, or by good engineering judgement. 

There are a number of theoretical sampling distributions 

which may be used, in preference to the experimental 

histogram, to model the data from (a), although the 
normal approximation is commonly employed’. The fact 
that the ash data are not independent is not of concern 

Since it will not affect the calculations. The 

distribution is simply a convenient way of summarizing 

the observed data. If f(a) is the probability density 
function describing the distribution of the measured 
product ash content such that: 


>The validity of this assumption should be be checked using 
a x* test. 
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"instantaneous" product ash. 
arithmetic mean product ash derived from the 
Sampling campaign. 
o(a) =the standard deviation of a, also derived 
from the sample data. 

N{f(a);a,o(a)} =a normally distributed variate, 
f(a), with mean, @, and standard deviation, 
a(a). 
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the cumulative product ash from the system is calculated 


as: 
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where: A =cumulative product ash. 
Y(a) =the functional relationship between yield 
and ash. 
Now, if one speculates on the fractional variance 
reduction, r, then a new distribution function f'(a) is 


defined where: 


Nif*(a)sa'. /ro(a)J 
Wiech Oest rs So 


Re-solving equation (2.13) under the conditions a'=a 
provides a new value for the overall product ash, A', 
which will be less than A. To satisfy the condition that 
A'=A requires an iterative solution for a', and, since 


this a' will be greater than a, the average yield 
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(denominator in eqn. (2.13)) will increase. The yield 
difference is the estimate required in the previous 
sub-section. A word of caution is necessary. It is very 
important that the relationship Y(a) be as accurate as 
possible since the yield difference is very Sensitive to 
this function. 

To illustrate, a local metallurgical coal plant 
treats 2.8 million tonnes of raw coal per year. There 
are two circuits in the plant, one for the coarse coal 
(>0.6mm) and one for the fine coal. The fine material is 
prepared entirely by froth flotation and represents, on 
average, 35% of the plant feed. A sampling campaign on 
the froth product from the fines circuit gave the 
following results a=8.75%, and o(a)=1.5%. Composite 
Sampling over a bank of flotation cells provided the 
estimate of the functional form of Y(a) as: 

70 .1+7 56 (a-8.625) -1.24(a-8.625)* + a<9.75 
va) na 24) 


80.0+0.1(a-9.75) ; a>9.75 


Assuming that f(a) is normally distributed and solving 
equation (2.13) by numerical methods yields; A=8.99% and 
Y=68.3%. If a variance. reduction of 25% (i.e. r=0.75) 
can be reasonably expected, the iterative solution of 
equation (2.13) for the case A'=A=8.99% yields; a'=8.82% 
and Y'=69.6%. The expected yield increase is 1.3%. At 


the production levels specified and using a coal price 


LEY | Bol? ites an id its 
) Soe 1 se +05 Bete iwi Re ade 

yoey es ot Qe atasen Ro iii 
ab etrieds 26 ot Veg ine 4b. 
of gui Kites wexew a). Pasovet® [oye 
; Ml 4 a 7 a4 a 
Maa HOD BOLT Nd lbin te ‘ede Sins a. 
eran? «whey, ded lags yeux? sh Sea at 


eo S29RaS ase ol ns Re. 


fo pl Sages 2 Agsmins: 
ar % Ss =4 


osha ott jamb toubieagy tapas 
er Lengo s «<2? ; a ia y4 i bie Fs a i ear, ches giluaer ‘end a th 


old. Qetsyor elf a ntl Pee Ate! nage ne s pat ls na 
sam (A)% +e fr: tad a8 30.9 be i? 
rie ot 
Y z Bee, 6, ay oN } 
Let 6 S ) >a es ; ! ’ « | a Pee ri m 
at had od | et a 


griviow duie hedudistath gone, et id We. oistta privy, es. 
Mine 206.84R f2B lo iy aN ban. ever Ts tet SY sor séupe, 
(20.@e4) a. 7). 22K Yo houses, opléi ai 5, is i ‘ 7 

oe bie Te ‘1 

20 Aclguics coip Gadi att! ote Biaiik sengaan Teac 


PSB.£= ad sob i ping fay vio. Hee ae, ah 102) hee ra) sdboaupe, oy ss 


‘yee! s 


6 WAC. t eh eesarsoi§ rae i Beta, oat rs erso Hiale 


36 


of $50/tonne f.0.b the mine gate, the increase in yearly 
cash flow before tax is $637,000. This neglects any 
additional operating and capital costs which may be 
incurred, although these are expected to be small. The 
economics suggest that there is good reason to improve 
upon process control in the fine coal preparation 


eircousit. 


2.8 AN IMPROVED METHOD OF ESTIMATING THE POTENTIAL REDUCTION 

IN ASH QUALITY VARIANCE 

In the previous sub-section a method of estimating the 
economic impact of reducing the product ash variance was 
desevssedigit nequiredithatra weasonabilei estimatesof 
variance reduction be made. However, in the absence of 
detailed knowledge of the time trajectory of the ash content 
the resulting estimates can provide grossly misleading 
results. For example, if the ash "Signal" contains primarily 
very high frequency "noise" then little, if any, variance 
reduction can be effected by conventional feedback control 
systems. (In this case the investigator must establish what 
the "generator" is for this noise and attempt to deal with 
it first.) This section provides a method of making 
reasonable estimates of the potential variance reduction 
from extrinsic disturbances using automatic control systems. 
In addition, the information extracted can be used to 
speculate on the technical feasibilty of such control 


systems. It is discussed in a "recipe" format to keep the 
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treatment as brief as possible. 

The technique to be described is part of the general 
subject of "Spectral Analysis." There are a multitude of 
books on the ene Cece hue the one by Rayner(13) is 
particularly readable, even for the neophyte. In essence, 
the method involves fitting a Fourier series to suitably 
transformed data and then estimating the variance components 
according to frequency. Low frequency contributions of less 
than one cycle per hour are expected to be virtually 
eliminated by automatic feedback control systems. Thus 
variance reduction can be be more accurately estimated and, 
in addition, the amenability of the process to control can 
be evaluated. Another feature of the variance versus 
frequency spectrum is the identification of "noise" 
generators, conspicuous by peaks in the spectrum. This may 
provide the operator with important "clues" in the search 
for physical causes of process upset. 

The method is taken almost directly from Rayner's book. 
Since this represents a condensation of the underlying 
theory and because additional comments are required relative 
to this application, the method is given in detail rather 
than simply by reference. To illustrate the results of 
calculating a variance spectrum and to provide a test case 


for the reader an example is included. 


Step (1):Gathering the Raw Data 
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For the purposes of this analysis the operator is 
required to sample the process stream of interest at regular 
intervals for an extended period of time. For statistical 
‘integrity of the data at least 100 samples (see Wells-Chap. 
3) must be collected. Ideally the entire exercise would be 
repeated several times although the variance spectrum is 
often observed to be more stable than the Fourier 
coefficient spectrum. For coal preparation units, or 
circuits, or plants, sampling frequencies of greater than 1 
Sample per h, typically 3 samples per h, are recommended. In 
addition, it is desirable to make some high frequency cuts 
(say >6 samples per h) prior to finalizing the sampling 
schedule. This allows an estimate of the contribution from 
the very high frequency component and is useful in defining 
the sampling interval. High frequency contributions which 
are present but not accounted for will manifest themselves 
through the phenomenon of aliasing (Superimposed on lower 
frequencies) and provide an incorrect spectrum. In coal 
preparation this component is often observed to be 
negligible. 

The data must be contiguous in time and obtained at 
regular intervals. If the latter can not be strictly adhered 
to in all samples, interpolation may be used. If mechanical 
interruptions are experienced, both the control as well as 
the process variables must be set to the previous levels 
after start-up and sufficient run-in time allowed before 


restarting the sampling. If the ash level varies 
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Significantly over the interruption the analysis can not be 
applied. If there are too many mechanical interruptions then 
perhaps the biggest problem in plant operation is not 
‘process stabilization. 

Symbolically, the data points are referred to as the 
set of (j,a,;); j=0,1,...,d-1, where a total of d samples 


were taken. 


Step (2) Remove the Mean and Long Term Trends 

The mean and long term linear trend are removed from 
the ash data to reduce them to a perturbation measure. This 
is required to suppress their influences from appearing in 
the zero frequency component which can also blur the higher 
frequencies. This step is accomplished by first regressing 


a; on j uSing the simple linear model: 
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where; bo,b;, =the unweighted least squares estimate of 
a the intercept and slope respectively. 
a; =the predicted value of ash at "time" j. 
If the regression is significant, i.e. if an "F-test" shows 
that the true slope is different from zero, then the raw ash 


data is transformed as: 
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where: a,;' =perturbation quantity. 


Step (3) Applying the Window Function (Tapering) 

The perturbation quantity from step (2) must be 
multiplied by a window function to render the theoretical 
development of the variance equations valid. The one 
Suggested by Rayner is Tukey's cosine bell. Here, the first 
10% (an arbitrary value in the range 5%-25%) of the data are 
multiplied by a rising cosine bell and the final 10% by a 
falling cosine bell. Let g=d/10, then this function will 


have the form: 


0 Std Costas G) 5.50 =) <1¢ 
idl ible ES A oe a 


i etts (ae to) 
0.5{1-cos(m[d-j] /g)}; d-g< j <d-1 
where: h,; =window (filter) function at time j. 
and the new variable to be used in the subsequent 
calculations is aa;, where: 
aa.=a.'h, 
aes eee ea hieten ios) 


This second transformation has the effect of reducing the 
overall variance somewhat but should not alter the variance 
spectrum appreciably. The result of this operation is that 
the estimate of variance at a particular frequency is not a 


discrete estimate, rather it is the estimate for a band 
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centred at that frequency. 


Step (4) Adding Zeroes 

Zeroes may be added to the end of the data record of 
aa;. This is done to adjust the frequency bands in the 
variance histogram to the desired widths and central 
frequencies (see Step (6)). This "cosmetic" alteration 
reduces the overall variance by a predictable amount but 
leaves the spectrum unaffected. 

The number of zeroes to be added is made equal to n-d 


such that a total of n entries are in the aa; data record. 


Step (5) Compute the Fourier Coefficients 
The data record from step (4) can be fitted with a 


Fourier series of the form: 
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where the Fourier coefficients are calculated from the 


expressions'°: 
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'°Other, faster methods may be employed to expedite these 
calculations. 
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data points. 


Step (6) Calculation of the Variance Spectrum 
For periodic or non-periodic data'' the variance 
contributions from the various frequencies (k=0,1,...,kmax 


cycles per sampling period) are calculated as: 


. Z 
V = 0] k=0 
OT Lea 
ee4s- 
V, =| k x | tasks - 
[ 2 amas Mee (22a) 
Viemax [ ise é OOS Aye 
where: V, =Vvariance contribution from frequency k. 


These values are collected according to the standard rules 
of histogram generation and displayed in this form as the 
variance spectrum. 

Table 2.1 and Figure 2.5(A) contain some hypothetical plant 
data for the previously discussed flotation circuit 
(sub=section. 2.7.3.2). It is hypothetical in the sense that 
it represents two contiguous records obtained on consecutive 
days and "pieced" together for this example. Samples of the 
flotation concentrate were taken every 20 minutes and the 


solids analysed for ash. In step (2), the regression was 


''In the theoretical development of this method the data is 
assumed to be non-periodic. There are few truly periodic 
functions in nature. 
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Table 2.1 Hypothetical Data for Spectral Analysis 
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insignificant and the raw data were transformed using 
equation (2.17). In step (4), twenty zeroes were added to 
the record for a total of n=100 points. This was done to 
have the class boundaries in the spectrum histogram integer 
valued. The variance spectrum is shown in Figure 2.5(B) and 


given in Table 2.2. 


Table 2.2 The Variance Spectrum for the Hypothetical Data 


CENTRAL VARIANCE PERCENT 
FREQUENCY COMPONENT CONTRIBUTION 


The estimate of variance reduction, assuming that components 
with a frequency lower than about one per one and one-half 
hours can be eliminated, is 92%. If this were real data, 
then, since the sampling period was quite short, i.e. <27 h, 
it would probably be fair to say that the observed 
trajectory is predominantly the effect of extrinsic 
disturbances. Thus, at least on the basis of this analysis, 
process control would have a very significant effect in 
Stabilizing the operation and reducing direct production 


losses. 
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Figure 2.5 The Hypothetical Data and the Resulting Variance 


Spectrum. 
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2.9 SUMMARY 

It has been shown that, in theory, it is desirable to 
have each circuit within a coal preparation plant operating 
at steady state with a constant incremental ash in the clean 
coal products, subject to an overall product ash constraint. 
In practice this is extremely difficult to do and it is 
observed that the operator specifies control set points on 
product ash, as an indirect measure of incremental ash, in 
his attempt to optimize plant performance. However, the lack 
of effective plant control systems results in Significant 
variations in product quality and hence both direct and 
indirect production losses are experienced. While some of 
this variation is "genuine", in the sense that it is caused 
by intrinsic disturbances, the balance is a result of 
extrinsic disturbances and must be eliminated. Methods have 
been presented for quantitatively assessing the technical 
feasibility and economic impact of control systems. There is 
considerable potential for improving process performance 
through optimization. This should lend some impetus to 
improving process control systems, principally through the 
development of sensors (instrumentation) and through 


increased process understanding. 
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3. PRELIMINARY IN-PLANT TESTING 


3.1 THE DATA ACQUISITION PROGRAM 

Before Asche aac the laboratory experimental program to 
evaluate the OSCAA prototype it was necessary to define some 
of the relevant operating parameters. Since McIntyre Mines 
Ltd. had been chosen to take part in the plant trial phase 
of the project it was decided that these values should be 
determined in this plant. A plant data acquisition campaign 
was conducted at the mine site in August, 1979. The detailed 
results of this work are reported elsewhere(14) anda 
summary of the relevant information is given below. For 
background information and because it will discussed again 
in Chapter 8, a brief description of the McIntyre operation 
is given. 

The McIntyre Mines Ltd. coal preparation plant is 
designed to produce approximately 2 million tonnes of clean 
coal per year. The product is shipped primarily to Japan for 
use in steel manufacture. The plant consists of two washing 
circuits. A heavy medium cyclone circuit is used to prepare 
the coarse coal (>0.6mm) and a froth flotation circuit 
treats the fine fraction. The fine material in the plant 
feed averages 35% by weight and excursions to 40% are not 
uncommon. A schematic of the flotation circuit is presented 
in Figure 3.1. It consists of four parallel and identical 
banks of four WEMCO #120 flotation cells. Process control is 


practiced on two levels. For short term modifications to the 
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Figure 3.1 The Flotation Section at McIntyre Mines (after 
McIntyre(15)). 
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control variables the operators use a visual interpretation 
of process behaviour. In this regard, pulp level seems to be 
the principal manipulated variable. Reagent addition rates 
' are modified on le from the quality control laboratory 
personnel. Decisions on these rates are made on the basis of 
daily composite sample ash values for the clean coal product 
in conjunction with historical trends in product quality. 
The target ash value for the flotation circuit is 8.5%. It 
appears that, as a result of the ill-defined nature of the 
flotation process, fine tuning on plant product quality is 
done principally by adjustments in the coarse coal circuit. 
When compared to other metallurgical coal plants 
treating coal in the same general geological setting, 
McIntyre's fine coal circuit is uniquely simple. This single 
unit operation approach is possible because the "oily" 
nature of their coal gives it a high inherent floatability. 
However, blocks of oxidized coal which are encountered from 
time to time are essentially lost to tailings under this 
processing scheme. In other plants, the coarser fraction of 
the fine material (typically -.6mm < X < .15mm) is treated 
uSing the automedium cyclone which effectively separates on 


the basis of particle specific gravity (see Figure 1.1). 


3.1.1 Fine Coal Preparation Circuit Behaviour 
The purpose of this portion of the data acquisition 
campaign was to observe the flotation circuit ash quality 


both on a long term (daily composite) and a short term (=20 
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minute intervals) basis. The results were used to define the 
need for improved process control as well as to establish 
calibration ranges for the laboratory investigation of 

' OSCAA. The "long term" data were obtained from the Quality 
Control Lab. at McIntyre while for the latter, a particular 
bank of flotation cells was sampled. 

Figure 3.2 illustrates the variation in the fine coal 
product ash over a 30 day production period in the summer of 
1979. The quality control target of 8.5% is shown as the 
dashed line. The samples from which the data were derived 
represent a 24 hour composite with a total of 3 cuts, one 
per shift. it isyevident by inspection that there 1s a 
Significant fluctuation in quality around the target value. 
The two moments of interest calculated from the data are; 
mean ash=8.9%, and the variance=0.55%*. Not only are there 
fluctuations in quality, but the long term observed mean 
does not correspond to the set point. It is impossible to 
determine whether these deviations are due to intrinsic or 
extrinsic (disturbances however the tect that a target 
value is defined and very seldom met is a comment on the 
control problem using existing methods. For example, there 
are several consecutive-day periods where a substantially 
constant offset was experienced with no apparent success in 
modifying circuit performance. These results demonstrate a 
need for improved process control methods in the McIntyre 


fine coal preparation circuit. 
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Figure 3.2 The Daily Variation in The Fine Circuit Product 
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High frequency sampling was carried out on the number 
four flotation cell bank for an 8 hour period on consecutive 
days. (The lack of man-power precluded a sampling exercise 
which would have provided the data record required for the 
analysis method described in section 2.8) The objective was 
to establish the frequency and amplitude of the fluctuations 
of selected stream attributes (%solids, %ash) around the 
cell bank. Initially, samples were obtained for the feed, 
concentrate and tailings streams, however, the extreme 
difficulty in filtering the latter required that it be 
excluded from the high frequency testwork. The results are 
presented in Figures 3.3 through 3.6. For approximately two 
hours on August 9, the sampling interval was reduced from 20 
minutes to 10 minutes. The purpose was to identify the very 
high frequency components. (The error bars represent the 95% 
confidence limits assuming a normally distributed variate. 
The variance was a pooled estimate from duplicate cuts.) The 
plant throughput was essentially constant on both days with 
the exception of a mechanical interruption at noon on August 
a. 

It is readily apparent from the figures illustrating 
the attribute variation for the feed stream that any trend 
is effectively masked by the sampling and analysis error. 
This relatively large error arises because of the method 
used to sample the feed stream. (A dipping bucket is 
submerged in the distribution box (Figure 3.1) to obtain the 


Sample.) Despite any apparent trend in the measured 
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Figure 3.3 High Frequency Sampling Results for the Flotation 


Feed-August 8, 1979. 
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Figure 3.4 High Frequency Sampling Results for the Flotation 


Concentrate-August 8, 1979. 
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Figure 3.5 High Frequency Sampling Results for the Flotation 


Feed-August 9, 1979. 
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Figure 3.6 High Frequency Sampling Results for the Flotation 


Concentrate-August 9, 1979. 
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variables it was suspected that the deviations were 
primarily the result of sampling error. 

To test this suspicion the following calculations were 
done. Each Seat eeaon each day was assumed to be a 
normally distributed variate. Using Chauvenet's Criterion, 
outliers were rejected. Only in one of the two data sets 
were observations excluded and in this instance only two 
points were rejected. Using the data subset from outlier 
rejection the normality assumption was verified uSing a 
x*test. Subsequently, the measured variable variance was 
tested for homogeneity against the sampling and analysis 
variance obtained from the duplicate cuts. In all cases the 
calculated "F" statistic was significantly lower than the 
tabulated value at the 95% confidence level. Thus, it was 
concluded that, for the feed stream, the % solids and % ash 
values were substantially constant over the high frequency 
Sampling periods. 

The measured attributes in the concentrate stream do 
exhibit a considerable variation in time with minimal 
contribution from the high frequency components. Since 
product quality is the principal control variable, the 
discussion will be restricted to the % ash data of Figures 
3.4(B) and 3.6(B). Again, as in the previous data, the mean 
operating level did not correspond to the control target. 
The quality trend is one which makes the process potentially 
amenable to process control. Furthermore, since these test 


periods were relatively short, there is a high probability 
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that the deviations are due primarily to extrinsic 
disturbances. 

From the figures the following typical operating 
conditions were defined. These were used together with other 
available information in setting the factor levels for the 
laboratory experimental program(see Chap. 6). 


Table 3.1 The High Frequency Statistics for %Ash and 
*Solids. 


STREAM % SOLIDS % ASH 


FEED 8.9%+3.4% 1%+4.9% 
CONC. 25.5%+5% bf ieee 3% 


The intervals quoted represent the 99% confidence limits for 


a normally distributed variate. It is recognized that the 
mean values are probably not good overall estimates of the 
true mean (See Figure 3.2). However, the variances are 
probably reasonable estimates since this statistic is 
generally more homogeneous. | 

Even though this data was obtained from McIntyre Mines 
Ltd., based upon discussions with other operators the 
process behaviour illustrated above is typical and is by no 


means an extreme case. 


3.1.2 The Yield-Ash Curve in Coal Flotation 
A knowledge of the yield-ash relationship is important 
in estimating the economic impact of improved process 


control in coal preparation. To gain an appreciation of the 
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nature of this curve in coal flotation a "kinetic" study was 
done. (It is implicit that other independent variables are 
expected to have an effect on the relationship somewhat 
Similar to that of flotation time.) Composite samples were 
taken from the number four flotation cell bank on August 9 
from 4:30 pm to 5:30 pm. In addition to the cell bank feed, 
concentrate and tailings streams, samples of the concentrate 
and tailings streams from each cell were obtained. Each 
sample was analysed for %solids and %ash. Using a general 
least squares solution to the overdetermined material 
balance problem (MATBAL 2), the measured data was adjusted 
to provide a consistent material balance(16). Using the 
adjusted data the cumulative yield-cumulative ash curve was 
constructed and is shown in Figure 3.7(A). 

The data points on the curve correspond to the 
different sequential groupings of the flotation cells. For 
comparative purposes the target ash and expected yield 
values are included. It is apparent that under the operating 
conditions which prevailed during sampling, the final two 
flotation cells were serving to contaminate the concentrate 
product, relative to the target. This suggests that air flow 
may be a useful control variable in a process control 
scheme. 

Using a quadratic approximation to the curve in Figure 
3.7(A) and the appropriate numerical methods, it was 
possible to construct the incremental ash versus cumulative 


ash curve. This is shown in Figure 3.7(B). It is quite clear 
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Figure 3.7 The Yield-Ash Curve and Incremental Ash in Coal 


Flotation. 
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that if the target of 8.5% ash represents the optimal 
operation relative to incremental ash, the actual operating 


point is grossly in error. 


3.1.3 Nominal Residence Time of the Flotation Circuit 

To assist in defining the "turn-around" time for an 
OSCAA system set of attribute analyses, it was necessary to 
determine the nominal (or mean) residence time of a 
flotation cell bank. This parameter is a measure of the time 
scale of the process and is frequently used in process 
control loop tuning. 

For the purposes of this experiment the following 
assumptions were made; 

a) The residence time distribution (RTD) of the solid 
particles and the water are identical in the cell 
bank. This allows the overall nominal residence time 
to be inferred from experimental data obtained on 
water. 

b) The flotation cell bank functions as a single 
input-single output (SISO) system. This allows 
Sampling on the flotation tailings stream only, 
neglecting water losses to the concentrate stream. 

The standard impulse tracer test was performed on the 

number four flotation cell bank at McIntyre Mines Ltd. 

NaCl (specifically Na*) was used as the tracer. The 

experimental RTD profile is shown in Figure 3.8 where 


Na* concentration is plotted as function of the time 
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after the impulse addition to the feedstream. 

The estimation of the nominal residence time, 7, 
from the experimental profile is straightforward and is 
only briefly documented in the discussion below. For 
more background the reader is referred to any basic 
chemical engineering textbook dealing with process 
analysis. 

If a quantity of tracer, M, is added as an impulse 
to the feed end of a constant holdup SISO system 
operating at steady state with a constant volumetric 
Enpoughnput, Gd, then 1f'c(t) isthe concentration of 
tracer in the vessel discharge stream at time t, the 
fraction of material: having Spent a time, t,; in the 


vessel is given by: 
E(t) dt=c(t)qdt 
M Pers cals) 
where: E(t)dt =fraction of feed stream material 
spending a time t in the vessel. 
c(t) =tracer concentration in the vessel 
discharge stream. 
q =volumetric flow rate through the vessel. 
dt =infinitesimal time increment. 


M =quantity of the traced species added in the 
impulse. 


The function E(t) is referred to as the exit age 
distribution function and has the character of a 


probability density function with; 
+00 
E(t)dt=1 eesti) 
0 
The nominal residence time is computed by the standard 


methods of statistics as: 
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Figure 3.8 The RTD for Water in the Number Four Flotation 


Cell Bank. 
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+00 


r= ftE(tyat 


0 
where: 7 =nominal (mean) residence time. 


BOE BEN, 


Applying the appropriate numerical methods to the data 
of Figure 3.8 the solution of equation (3.3) gives an 
estimate of 7r=4.8 minutes, under the experimental flow 
conditions. These conditions, specifically the 
volumetric flow rate, were observed by the operators to 
be typical of plant operation. 

Control system designers would prefer to have 
process measurements available with a sampling interval 
of roughly one-tenth of the nominal residence time. 
However, in this case, given the practical constraints 
of the sampling and measurement system, the nominal 
residence time and the product quality variations of 
Pigures o.4(B) and 93 6(6) (sit is evidentiithat an 
interval of up to 20 minutes could be tolerated without 
seriously impairing the control system performance. Thus 
an OSCAA system should be capable of providing a set of 
analyses on an interval of from 5 to 20 minutes to be 
useful in the process control of the McIntyre flotation 


Circuit. 


3.2 SUMMARY 
There is a Substantial day-to-day variation in the 
product ash for the McIntyre fine coal preparation circuit. 


These daily variations arise from relatively low frequency 
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variations in the flotation concentrate product quality. In 
principle it should be possible to practice effective 
process control around this circuit, provided that an OSCAA 
type sensor is Judi deb hey The response time of this 
analytical device should be in the range of from five to 


twenty minutes. 
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4. A REVIEW OF ON-STREAM ASH ANALYSIS TECHNIQUES 


4.1 THE FUNDAMENTALS OF ON-LINE ASH MEASUREMENT 

The principal methods for the instrumental 
determination of ash which have been considered so far are 
based either on the measurement of specific gravity, on the 
effective atomic number of the coal measured by 6°, X, or y 
radiation, or on neutron activation. These methods are 
non-destructive and rapid (relative to the more conventional 
off-line methods), but are indirect and hence assume a 
strong correlation between ash and mineral matter as well as 
between mineral matter and the parameters being measured. 

With respect to the nuclear techniques, Clayton(17) has 
provided a good general overview of the applications in the 
coal industry which deals, in part, with the on-line 


determination of ash content. 


4.2 PHYSICAL MODELS OF THE PROCESS FLOWSTREAM 

In most of the analytical methods employed to estimate 
the solids ash content, some physical model of the process 
flow stream is either implicitly or explicitly adopted for 
the purposes of calibration. The model generally assumes 
that the process flow stream consists of from two to four 
phases, depending upon the exact nature of the measurement 
system and the application. Associated with each phase are a 
set of physical and chemical properties which can be 


employed in constructing the calibration model. 
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To illustrate, consider the simplest case of a dry 
solids flow stream. In a baSic measurement system the solids 
would be assumed to consist of a coal organics phase and a 
-mineral matter phase. Each phase would have its own mean 
chemical composition, density, etc... These phases and their 
properties would be the basis of a calibration model (as 
will be illustrated in the next section). Clearly, the 
models are a Simplification of the real system, but a 
necessary one since there are usually only several (<3) 
independent measurements made on the flow stream. Since, in 
all models, the mass fraction of the phases must sum to 
unity, one less equation (or independent measurement) than 
phases is required to exactly compute flow stream 
composition in terms of the model phases. Thus, for the 
example, a single measurement is required. Any change in the 
petrography or mineralogy of either of the phases will 
invalidate the physical model and require that the 
characteristic properties be re-determined. 

Finally, consider the most complex case of a fine coal 
Slurry flow stream where an iron correction is required to 
circumvent frequent re-calibration (i.e. estimation of the 
phase properties). Here the flow stream would be assumed to 
consist of four phases, water, coal organics, iron-free 


mineral matter, and some characteristic iron mineral. 
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4.3 DENSITY MEASUREMENT 

There is a Significant density difference between the 
organic portion of the coal (S.G.=1.35) and the mineral 
“matter (S.G.=2.6). Assuming a two phase physical model of 
the coal solids, with the further assumption that the 
mineral matter to ash ratio is constant, permits the 


construction of a calibration model based on measured solids 


density: 
_ Winog [| ia 
Rma = 
Pmm °o} Ls padi al dhe) 
where: a= ash content of the coal solids. 
Rma=mineral matter to ash ratio for the coal 
solids. 


Ps = measured coal solids density. 
Po =coOal organics density. 
Pmm =Mineral matter density. 

The model parameters (phase properties) can be estimated 
from calibration data by standard curve fitting techniques. 
In some cases, where the range of a is small, empirical 
approximations of the form a=f(p,), (f£(p,) linear or 
quadratic), are used. In essence these are simply truncated 
Taylor Series expansions of the more general model given by 
equation (4.1). Of course, changes in petrography or 
mineralogy, which can and do occur quite frequently, 
invalidate the assumed physical model since the phase 
properties change. 

Three types of coal-ash monitoring devices based on 


density measurement have been developed. The first is simply 


batch weighing using water or some organic liquid as the 
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displacement medium. This method, while not an on-stream 
technique in the strictest sense, can be quickly performed 
by the operator to provide an ash content for process 
control purposes. The second method, which has been 
described by Crosland(18), measures the bulk density of an 
essentially dry coal stream. Assuming a constant void 
fraction, the coal density is calculable. The third type of 
monitor involves measurements made at the feed end of a coal 
preparation plant, where the solids are slurried with water 
prior to wet processing. This device has been described by 
Mikhail(19) and, to date, has only been tested ina pilot 
plant setting. This unit requires a mass flow meaSurement of 
both solids and water flow to the plant. Subsequently a 
Slurry density meaSurement is made and the solids density 
back-calculated. 

These devices have the advantage of requiring minimal 
sample preparation, however, they are very dependent upon a 
Stable ash-density correlation. In addition, they cannot be 
applied in systems where the feed is already in the form of 


a Slurry. 


4,4 NEUTRON ACTIVATION 

This kind of method, because of its complexity and 
expense, has seen only limited application. Those mies 
which are in operation deal exclusively with monitoring the 
feed to thermal power plants. Rogers et. al.(20) have 


described a particular unit known commercially as the CONAC. 
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In principle, a uniform bed of coal is irradiated with fast 
neutrons most of which experience deccelerating collisions 
with the nuclei of the elements in the coal solids. When the 
neutrons are reduced to a sufficiently low energy they are 
captured by a nucleus creating an unstable isotope. The 
instability of the nucleus results in the prompt emission of 
a y-ray(s). These gamma rays have a characteristic energy 
and their intensity provides the means of estimating 
chemical composition. 

A similar device has been proposed by Sowerby(21) which 
dispenses with the sample metering requirement and employs a 
separate y-ray backscatter measurement to correct for the 


coal sample bulk density and thickness. 


4.5 TRANSMISSION AND BACKSCATTER OF ELECTROMAGNETIC 

RADIATION 

These methods employ atomic number inference to 
estimate the ash content of the coal solids. Im tacth. ltteers 
the mass absorption coefficient which is measured but which 
is related to the atomic number. There is a significant 
difference between the atomic number of the coal organics 
(Z=6) and that of the mineral matter (Z=12). Assuming a two 
phase physical model of the coal solids, changes in 
composition are reflected by changes in the measured 
effective atomic number of the coal. Since the derivation of 
the calibration model pre-supposes a knowledge of both the 


instrument design as well as the means by which 
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electromagnetic radiation interacts with matter, discussion 
on this topic is deferred to the next chapter. As with the 
density methods, the mineral matter to ash ratio is assumed 
to be constant. Of course, changes in the chemical 
composition of these phases, particularly heavy elements in 
the mineral matter, will alter the effective atomic number 
of the component and thus invalidate the assumed physical 
model. 

This general method of coal-ash analysis appears to 
have been the most widely investigated to date. 
Specifically, those methods which employ y-ray and/or X-ray 
irradiation have received the greatest attention because 
they are Superior to the 6” techniques. Calculation of the 
ash content of coal through correlation with the atomic 
number is more reliable and applicable than the density 
methods and much less expensive than neutron activation. The 
major difference in the implementation is the geometrical 
arrangement of the irradiation/detection equipment of the 
ash-analysis subsystem. 

Figure 4.1 illustrates the differences between 
transmission and backscatter geometry. In transmission, the 
radiation source and detector are located on opposite sides 
of the sample. As a result, the majority of the radiation 
which reaches the detector are primary photons which have 
passed through the sample without interaction. In 
backscatter, the source and detector are located on the same 


Side of the sample. As a result, all of the radiation which 
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reaches the detector is the result of an interaction of the 
primary photons with the sample. With regard to classifying 
the ash analyser, the device geometry is taken from that 
part of the instrument which is sensitive to atomic number 
(or mass absorption), i.e. the ash analysis sub-system. 

These methods can be further subdivided on the basis of 
application. One class of instruments is applied to coarse, 
essentially dry solids while the other is applied to fine 
@peacticie slurries. The existence of high concentrations of 
water, which has both a density and an effective atomic 
number similar to the coal organics, complicates the 
measurement task appreciably. A further complication which 
can arise in slurry applications is the presence of varying 
quantities of entrained air in the pulp. This will distort 
the pulp density measurement which is usually made with a 
sub-system in the transmission geometry. AS a result, much 
of the development work has been done on the dry solids 
analysers and only recently have these ideas been 
transferred to the slurry measurement problem. 

The underlying theory of operation of all of these 
on-stream analysis devices is presented in the ee chapter 
where the physics of the OSCAA prototype are described in 
detail. 

There is a relatively large number of off-line 
implementations of these techniques discussed in the 


published literature as well(for example:22-27). 
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4.5.1 Transmission Methods 

In general, transmission methods are well suited to 
this analysis problem because they are very sensitive to 
changes in atomic HUNBET 4 They are limited in flexibility, 
however, because additional radiation measurements, which 
may be necessary for compensation purposes, require 


duplication of the basic transmission measurement equipment. 


4.5.1.1 Dry Solids Applications 

Three classes of devices have been successfully 
applied to essentially dry solids, at least on the 
laboratory prototype scale. The first class involves a 
Single transmission measurement. Starnes and Clark(28) 
and Nagy and Varga(29) describe systems which involve 
fairly elaborate sample handling (comminution) and 
Sample metering equipment control to ensure sample 
homogeneity and constant mass per unit area 
respectively. The single, low energy photon, attenuation 
measurement provides an estimate of the mass absorption 
coefficient (atomic number) which is correlated to the 
ash content for calibration purposes. The second class 
involves a dual transmission system as decribed by 
Fushimi(30). In this case, a high energy photon 
attenuation measurement is also obtained to correct for 
the sample bulk density. This obviates the difficult 
task of controlling mass per unit area, however, size 
reduction is required to provide uniformity of field. 


(This device operates in a batch mode.) The third class, 
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described by Fookes et. al.(31), employs three 
transmission measurements. Their unit 1S virtually the 
Same aS Fushimi's but includes a second, low energy 
photon, attenuation measurement to partially compensate 
for possible mineralogical variation in the mineral 


matter. 


4.5.1.2 Slurry Applications 

Two devices have been proposed/applied to the ash 
monitoring problem for coal slurries. Starnes and Clarke 
altered the design of their dry solids unit to 
accommodate fine particle slurries. They included a slow 
Speed agitator in the supply tank to facilitate 
deaeration. In addition, a pulp density gauge was added 
to the system and the flow cell modified to allow for 
both transmission and fluorescence radiation 
measurements, the latter for compensation purposes. The 
author knows of no industrial applications of their 
device in this form. Lyman(32), using Fushimi's system 
aS a basis, included a high pressure pump to alleviate 
any problems which might arise from slurry aeration. 
Lyman's device, known as Automatic Slurry High Speed 
Coal Analyser (ASHSCAN), has been extensively field 
tested in Australia and is apparently ready for 
commercial distribution. Lyman has acknowledged the 
potential problem of compositional variance in the long 
term accuracy of ASHSCAN. He has suggested that a system 


Similar to the one described by Fookes et. al. would 
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offer a solution to this problem. Since ASHSCAN is 
typical of instruments based on transmission geometry, 
and because it is intended for the same application as 
OSCAA, a eoiensrane representation is included in Figure 
4.2. 

Interestingly, while Starnes and Clarke and Lyman 
make provision for slurry deaeration, work by Churchill 
et. al.(33) on lead ore slurries did not mention serious 


aeration problems. 


4.5.2 Backscatter Methods 

In general, while backscatter methods are less 
sensitive to changes in atomic number than are transmission 
methods, they appear to have been much more widely studied 
and developed, specifically for monitoring dry solids. A 
major reason for this is the relative ease of compensation 
through fluorescence meaSurements aS opposed to the 
duplication of equipment. This inherent flexibility of the 
backscatter geometry in slurry applications makes it a 


logical candidate for further development. 


4.5.2.1 Dry Solids Applications 

The first commercial unit of this class was the 
CENDREX developed by the Dutch State Mines(34). Since 
this unit was designed for on-stream application using 
low energy X-rays, a rather elaborate (ref. 28,29) 
Sample presentation system was required. CENDREX 


utilizes a conventional tube/target assembly to produce 
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Figure 4.2 A Schematic of Lyman's ASHSCAN Device(32). 
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the radiation. Crosland et. al. describe the evolution 
of a similar device which utilizes a radioisotope 
source. This work led to the development of the SORTEX 
unit which has been described by Cammack and Balint(11). 
Very similar work by Hardt(35) led to the development of 
another coal-ash monitor which is also commercially 
available. In all three systems, the intensity of the 
backscattered radiation is the measure of the coal 
solids mass absorption coefficient. The major 
differences among these three devices lie in the sample 
presentation sub-system, the energy of the radiation 
employed and the method of (iron) compensation. In all 
cases they are claimed to have been successfully applied 
to dry coal streams yielding the required accuracy in 
acceptable meaSurement time. 

Fookes et. al. have suggested that a dual source 
excitation system would provide good results in 
backscatter geometry as well. In this case the second > 
measurement would have the same function as fluorescence 
in compensation. They do not provide any experimental 


Support for this statement. 


A O26 2eoLUrrys ADDLicCations 

Thusfar, only one implementation of the backscatter 
geometry for monitoring coal-ash in fine particle 
Slurries has been suggested. Kawatra and Dalton(36) 
describe the successful calibration of a laboratory 


bench scale unit. As in ASHSCAN, two measurement 
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Sub-systems are incorporated. One to measure the pulp 
density and the other to meaSure backscatter and 
fluorescence radiation. This is the system upon which 
the OSCAA design was based'?. A schematic of the 
laboratory scale test rig used by Kawatra and Dalton is 


presented in Figure 4.3. 


4.6 IMPORTANT DESIGN CRITERIA IN SLURRY MEASUREMENT SYSTEMS 
In conducting this literature review in the light of 
the goals of the OSCAA project, several important design 
criteria for slurry application units were frequently 
encountered. These are briefly discussed below with 
reference to the analytical devices concerned. Special 
emphasis is given to their impact on the design of an OSCAA 


type system. 


4.6.1 Air Sensitivity 

In ASHSCAN, Lyman has eliminated the entrained air 
problem by pressurizing the system to the extent that the 
air bubbles collapse. On the other hand, Kawatra has not 
made an explicit attempt to deaerate. Since OSCAA, like 
ASHSCAN, requires an accurate pulp density measurement, the 
current prototype may need design modification to handle the 
air problem. This could include anything from using 
mild/strong centrifugal forces to the pressurization 


‘2Tt has come to the author's attention very. recently that 
Kawatra et al.(37) are essentially duplicating the OSCAA 
study but their work is not sufficiently completed to allow 
comparison of results. 
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Figure 4.3 A Schematic of the OSCAA Laboratory Prototype of 


Kawatra and Dalton(36). 


oe us 


i Jibs 


| 
wit 
; ; ut 
‘i WT na = 
i op ~ hy ee 
Roi. wil 
As 
ats i 
ot : 
et 
oe 
a? v 
, ‘er 
tie 
(er 


- sini TS ripe yy 
Mp imene Ph a Er ey, Ge Te ee 


SO) Ra ee La, jn 
‘er wer eye ene 


| “souag vriensel} are 
mareye-Guey] Rig er at 
a) He ae the i 
| mn or bead . ; = |} im otra 
0 naman signe. aestision’ hy sel | pie “ 

el te sais * 098 |: 


81 


technique employed in ASHSCAN. 


4.6.2 Sampling Characteristics 

Lyman's device utilizes transmission measurements with 
turbulent flow in near vertical pipe runs. Under these 
conditions essentially all particles in the flow stream are 
"Sampled" in the analysis. In the backscatter geometry, a 
flow cell equipped with a thin mylar viewing window is used. 
As a result only a "thin slice" of the pulp in the 
neighbourhood of this window is sampled. Proper design of 
the flow cell and the correct choice of the incident 
radiation energy will enhance the "Sampling" capability of 
such a system. This is the approach employed in the 
on-stream X-Ray Fluorescence (XRF) analysers common in 


mineral processing plants. 


4.6.3 Heavy Element Correction 

Changes in the concentration of the heavy elements in 
the mineral matter can have a devastating effect on the 
performance of Lyman's analyser. Theoretically, for every 
element whose concentration changes significantly another 
transmission measurement would be required. However, in the 
backscatter geometry, fluorescence measurements can be used 
to compensate for such changes with little modification to 
the system hardware, apart from the detection/spectral 


resolution equipment. 
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4.6.4 Compatibility/Flexibility 

In an operating plant, it is supposed that a multitude 
of streams will require monitoring. For a device like 
" ASHSCAN this will necessitate a considerable duplication of 
equipment and extensive stream multiplexing. The OSCAA 
System is quite compatible with existing on-stream XRF 
hardware. In these devices both multiplexing and duplication 
are minimized. It is expected that the manufacturers of such 
equipment could quickly produce a commercial version of 


OSCAA with excellent product support. 


4.7 SUMMARY 

There has been a substantial amount of development work 
in the area of on-stream coal ash analysis, in particular 
with those methods utilizing electromagnetic radiation. Most 
of this work has concentrated on the analysis of coarse dry 
solids. Unfortunately, since most of the preparation plants 
employ wet separation methods, this class of device is not 
well suited to feedback control applications. Furthermore, 
the mechanical problems associated with the sampling and 
Sample presentation sub-systems make them difficult to 
operate. Research into devices for use in fine particle 
Slurries has begun only recently. Lyman has developed the 
ASHSCAN system to the point where it is apparently ready for 
commercial distribution. The other geometrical variant in 
this area, OSCAA, has only recently been proven. Given the 


important design criteria for such devices, the backscatter 
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geometry would be the preferred implementation. However, it 
is possible that in some plants, the entrained air problem 
may be severe enough to warrant appropriate design 


‘modifications to the prototype described in this thesis. 
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5. THE PHYSICS OF THE OSCAA MEASUREMENT SYSTEM 


5.1 INTRODUCTION 

This chapter describes the physical principles upon 
which the OSCAA system is based and from which certain 
design parameters were deduced. In addition, the possible 
calibration equations for OSCAA are presented. 

In order that the material which follows be put into 
context, it is necessary to briefly discuss the OSCAA unit. 
(A more detailed description is given in Chapter 6.) OSCAA 
consists of three sub-systems, one for sample presentation, 
one for pulp density measurement and one for backscattered 
and fluorescent radiation measurements, together with their 
associated support equipment. The pulp density sub-system 
employs a Gd-153 radioisotope with a principal photon energy 
of ~100 kev, in the transmission geometry. The ash analysis 
(backscatter/fluorescence) sub-system employs a Cd-109 
radioisotope with a principal photon energy of =~22 kev. In 
this case, the slurry sample is irradiated through a thin 
mylar window. The resultant detectable backscattered 
incident and fluorescent radiation return through the same 
window. 

Since y-ray and X-ray photons are physically 


indistinguishable, the same laws of physics apply to both'?®. 


'3y-rays and X-rays are distinguished on the basis of their 
origin. y-rays result from nuclear transformations while 
X-rays result from electron rearrangement. 
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5.1.1 Practical Considerations in Measurement System Design 
and Operation 

This chapter is intended to develop the theory 
' underlying the OSCAA system and there has been no attempt to 
include any discussion on practical considerations which may 
influence system operation. While the performance of the 
detector and resolution equipment are critical factors in 
ensuring successful operation, these topics are sufficiently 
general that much has been written about them in the 
Standard textbooks. For example, Bertin(52) provides a good 
description of detector structure and operation, including 
choking, shifts, escape peaks, etc...(Chap. 6), as well as 
of resolution equipment (Chap. 8) and counting statistics 
(Chap. 11). The philosophy adopted by the author is that 
before one can begin to design and operate a system such as 
OSCAA, the theory of operation must be understood. The 
practical considerations become important in the final 
design and operation stages. Since the latter have been 
treated by many authors, material of this is nature excluded 


from the thesis to avoid needless duplication. 


5.2 THE INTERACTION OF PHOTONS WITH MATTER 
The mathematical model of the physical principles upon 
which nuclear transmission measurements rest is given by 


Lambert's (Beer's) law: 
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where: Ip =intensity of radiation in a narrow colinear 
beam which is incident upon the absorbing 
medium. 


I =the intensity of the incident colinear beam 
after transmission through the absorbing 
medium. 

p =the density of the absorbing medium. 

w =the mass absorption coetficient of the 
absorbing medium. 

x =the distance which the photon beam travels 
through the absorbing medium. 

This law is observed to hold for a narrow, perfectly 
collimated beam of mono-energetic photons. 

A parameter of particular interest in equation (5.1) is 
the mass absorption coefficient, uw. The magnitude of this 
Parameter is dependent upon the incident photon energy (E) 
and the effective atomic number (Z) of the absorbing medium. 
Miss isiillustrated in Bigure S.t where wu is plotted against 
Z for three photon energies. The range of Z chosen 
represents the range of major elements found in coal solids; 
hydrogen to iron. The effective atomic numbers of the coal 
organics (Z=6), and the mineral matter or ash (Z=12), are 
also indicated on Figure 5.1. 

In the photon energy range of interest in this work, 
i.e. ES700 kev, the mass absorption coefficient is the sum 
of the coefficients for three mutually exclusive events, 
namely; 

(i) photoelectric absorption (rT) 
(ii) Compton scattering (om) 


(iii) Rayleigh scattering (dnm) 


. Mathematically, 
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Figure 5.1 The Mass Absorption Coefficient as a FAInc elon: 


Photon Energy and Atomic Number. 
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These processes are briefly described below and are 


pechematicaily alijustrated™in Figure 5.2. 


Photoelectric Effect 

"An incident photon disappears and an atomic electron 
(usually from the proximity of the nucleus) leaves its atom 
at high speed, having absorbed the photon energy. This 
effect predominates for lower energy gamma rays, especially 
in high Z materials."(38) The electron rearrangement which 
follows this event gives rise to the phenomenon of X-ray 
fluorescence. When an electron in a higher energy orbital 
moves to fill the vacancy left by the photoelectron an X-ray 


of characteristic energy is emitted. 


Compton (Inelastic) Scattering 


"A photon is scattered inelastically and.an atomic electron 
recoils out of an atom, much as though it had been initially 
free. The energy taken up by the electron depends primarily 
on its recoil momentum. This effect predominates for gamma 
rays of 1-5 Mev in high Z materials, and even more greatly 
and over a much wider energy range in low Z materials. 

The recoil electron flies off predominantly in the 
direction of the incident photon, the more nearly so the 
higher the fraction of photon energy it carries away and the 


higher the photon energy. It never recoils backwards."(38) 
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Figure 5.2 The Absorption Processes for Low Energy Photons 


(E<i Mev). 
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Rayleigh (Elastic) Scattering 


"Gamma rays scattered by small angles impart only a small 
recoil, especially when their energy is low. The recoil is 
then absorbed by the whole atom or molecule so that the 
Scattering action of different atomic electrons combines 
coherently. Coherent combination enhances the probability of 
this effect so that it prevails greatly at low energies. 
This scattering is more likely for high Z than for low Z 


mMacerial sene3e)) 


Narrow beam meaSurements yield only u, but the coefficients 
for the individual processes are uSually provided by the 
theory or by. special experiments. Numerical values for these 
coefficients are available in tabular form (for example: 


C39 ))R3 


5.3 THE PRINCIPLES OF PULP DENSITY SUB-SYSTEM OPERATION 

Nuclear pulp density measurement has a relatively long 
history in both mineral processing and coal preparation 
applications. This type of pulp density gauge has a number 
of desirable features which include: 

a Pn Sieers tindirect:; 

b) it is mechanically reliable and geometrically 

flexible, 
c) it "samples" a large portion of the pulp stream. 
The installation of a nuclear pulp density gauge 


requires that the following factors be taken into 
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accouncs: 
a) avoid a buildup of one phase (e.g. solids) in the 
radiation pathway. 
b) ensure that the Slurry pipeline is running full. 
c) minimize entrained air in the pulp. 
ad) minimize pipeline vibration to avoid associated 
intrumentation problems. 
e) ensure adequate shielding to avoid radiation 
contamination. 
The first two requirements are usually met by mounting 
the gauge on a vertical section of the slurry pipeline 
where the direction of flow is upward. Point (c) is 
often dictated by the nature of the operation. However, 
this can be become an important factor in high accuracy 
measurement systems. The final points are partly 
compensated for when off-the-shelf instrumentation is 
used. Vendor equipment is generally quite rugged and 
must pass strict safety inspection. sa 
The two most common industrial gauge geometries are 
schematically illustrated in Figure 5.3. The reason for 
the different configurations arises from sensitivity 


considerations which are discussed later. 


5.3.1 Applying Lambert’s Law 
In slurry applications the absorbing medium is a 
multicomponent mixture, by definition. Thus, while equation 


(5.1) remains valid, provided that the implicit geometrical 
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Figure 5.3 Common Geometries of Industrial Nuclear Pulp 


Density Gauges. 
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constraints are met, the sample mass absorption coefficient 
must reflect the elemental composition since w=f(E,Z). It 
can be shown that the mass absorption coefficient for any 
phase is simply a weighted average of the coefficients for 
the elemental Species present. Using the same approach, if 
the mass absorption coefficient of each phase is known, the 
sample mass absorption coefficient is simply the weighted 


average value, i1.e.: 


u= DD Uw 


7 al aoe oess) 


where: wu, =the mass absorption coefficient for component i 
in the sample. 

w, =mass fraction of component i in the sample. 
Inspection of Figure 5.1 shows that by choosing a 
radioisotope of sufficiently high photon energy (e.g. E2100 
kev for coal), then uw is substantially independent of Z and 
can be considered constant. That iS, mu;=u2=...=u. Under 
these conditions and for a fixed path length, x, the only 
variable parameter in the exponential of equation (5.1) is 
p, the pulp density. 

This is the basis for the design of nuclear pulp 
density gauges where variations in the measured value of I 
are related, in a logarthmic manner, to changes in the 
sample pulp density. Linearization of equation (5.1) with 
respect to p permits gauge calibration with a minimum of two 
(I,o) points. 

Equations’ (5.1)7-(5.3) anos (5.4) aré-sufficient to 


construct a calibration model for Lyman's ASHSCAN decive. 
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iD 
o= 


2a Wi/y. sea 4) 
where: p; =the density of component i in the sample. 
This is useful since it serves as an introduction to the 
construction of semi-empirical calibration models based on 
nuclear measurement techniques. In addition, it provides a 
vehicle for theoretically studying the performance 
characteristics of ASHSCAN. From Figure 4.2 it 1S apparent 
that this device employs a pulp density gauge (Cs-137, E=660 
kev> uw effectively constant) to estimate p, and a lower 
energy photon transmission gauge (Am-241, E=60 kev> u 
dependent upon composition) to estimate pu. Denoting the two 


sub-systems as (1) and (2) respectively,then: 


Dey) 


aaene| nUNCe Ss 
Fo /-u a) (2) 


p=log, | 


pu (2)=Log, peed | 
5 2) IY-x (2) SEG) 


where: I(j) =the measured intensity of transmitted photons 
for sub-system j. 
Io(j) =the incident photon intensity for sub-system 
ae 
u(j) =the mass aborption coefficient corresponding 
to sub-system j. 
x(j) =the radiation path length for sub-system j. 


Since the terms on the right hand sides of these equations 


are either measured or known from the theory or experiments, 
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they can be used as predictive relationships. Equations 


(5.5) and (5.4) can be combined to form: 


Dw./. =-u(1)x(1) 
aaa tog | ie Ge | LAs Copal) 


Ty (1) 


Similarly, dividing equation (5.6) by (5.5) and.setting the 


result equal to equation (5.3) yields: 


2) 
_ H(1)x(1) log f yt (2)] 
dein? 28 66 (5.8) 


Equations (5.7) and (5.8) provide two simultaneous linear 
equations in the w;. ASSuming a three phase model for the 
Sluncyy (1) water, (2) coal organics and (3) mineral matter, 
and knowing that Zw,;=1, allows construction of the 


calibration model. In vector-matrix notation: 
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The stream attributes are computed from the solution for the 


right hand side of equation (5.9) as: 


s=100%(w.+w 5 5.1.6:) 
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where: s =% solids by weight in the slurry. 

a =% ash by weight in the solids. 

Rma =mineral matter to ash ratio for the coal. 

Using this calibration model, the sensitivity of 

ASHSCAN to changes in the mineral matter mineralogy, in 
particular the iron concentration, was theoretically 
evaluated (1.e. us waS varied). It was shown that for the 
McIntyre flotation circuit product streams, variations in 
the iron concentration of the mineral matter of +0.7% will 
contribute a 5% relative error in the estimated ash content. 
(This was also observed for an OSCAA type system where iron 
compensation was excluded.) Since iron concentration does 


change, ASHSCAN would be of limited utility unless 


‘frequently recalibrated. 


5.3.2 Deviations from Lambert's Law 

While the development of the theory underlying the 
nuclear pulp density gauge is reasonably straightforward 
using Lambert's law, the restrictions: for strict compliance 
are seldom met in practice. There are potentially three 
effects which can serve to complicate the measurement. Under 
"normal" circumstances they will seldom render the 
previously developed equations invalid, but the reader 


should at least be aware of their existence. 


savela eds o} haga ee eee 
: SSryue edt ad aie ee 
: Fad os eiel hes OE i 9 ee 7 a} Bs OF: r e * at ) - : 
ic TRAE CLE: ne ban nena | 
| : er ioe 
: 7 
nk wees eran tn pene Le 90 Rt: eae ; agua “7 
Lenivensad sei eaw 6 pe | sextoatingl tigoit wit, 3 


ett 403 texte vote dew 1F Unelisee aay ae ak} 


Poca SOV ees OY wang 5 Hsiao: nal saiold 


Litw 2) 02) 30) “aden Dawa tee os ‘to nok tergay: 

j ; n _ a 
Csioyh Hae Serene pee ot2'ink Aeate) evisaten 66 e 
a 


{Siw ine dee Says baba: fe a0, Boviaade oats 


es eke LE ce ‘ soak, (.Bebti fsx mee Rey. 
Via “pat bale dp, wid binow: Madan: 


ast taney itn 

Tu 2 a 

aie cls yey 

| vet | eh asia ‘eoad ‘sigatetede Be | 

ot? prey (reba yxoulsmdt "Se  peraheaet adaatiaw | es. 


nt 


i 


; i 


biewseaadp laid y idansnee ta) i epue tiene qlna sastou 


Ssrea! a G0 ' ] ~ ri - 23 iliated ‘olsdusy ods ent, “” +2edmeat po an ” | 
; \ Ve oo 

eared vijs 9408 ane conte 9okdgeag ak, ‘yom mobiee' © 5 

= 4! ey — 7 


106n0: dees ser Hy saan jqmon 6 av'iae FBS Hoidw ype?! ig 
eG: sebas: wth ie er 99a. aechatemosts ‘tamaen® 


nie 


o “eGae? eds tud Didi: cAI aupe: bogoteven yleuot¥eig 


smote. +heKa bree 9 howe od. ees 6: btueda, a 
ee . = . bia): 


oH, 


Sasi zed amie bration 

In some special instances the radioisotope chosen 
for the particular gauging problem will emit photons at 
two or more distinct energy levels. For example, Gd-153 
emits photons with energies of E,~100 kev and E2,=40 kev. 
In this case the lower energy photons would be 
attenuated to a greater extent by the absorbing medium. 
This phenomenon is known as filtration. Accordingly, as 
the pulp density and hence the pux term increase, the 
total photon attenuation will occur by smaller and 
smaller fractions. The linearized form of equation (5.1) 
will show an increasing slope as uw apparently decreases 
with increasing pw PE Ett vation -induceducurvature,is 
apparent from the calibration exercise, and remains 
unresolved, then the calibration model will require a 
series of terms on the right hand Side, Similar to the 
right hand side of equation (5.1). That is, one for each 
photon energy level which is contributing to the 


measured intensity. 


5.3.2.2 Build-up 


The: build-up factor, B, 1s defined as; 


p-Observed mcensity at the Detector 


Primary Intensity Expected at the Detector 


_ intensity due to Scattered Radiation 


Primary Intensity 
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Build-up results from the multiple scattering of 
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radiation in the absorber such that some of the photons 
report to the detector which otherwise, under narrow 
beam conditions, would not. Relaxation of source and 
particularly detector collimation will cause the 
build-up factor to increase. The multiple scattering 
process in "broad beam geometry" is illustrated in 
' Pigure 5.4(A). Figure 5.4(B) illustrates the variation 
of B with beam geometry. These results were obtained by 
Wyckoff et. al.(40) using concrete as the absorbing 
medium. (Note that in Figure 5.4(B), for a particular 
abcissa value, the build-up is simply the ratio of the 
corresponding untransformed ordinate values.) Since 
Compton Scattering 16) thevprincipal attenuation’ process 
for the photon energy levels encountered in nuclear pulp 
density gauges, build-up would be expected to affect 
measurement. This is particularly true given the broad 
beam geometry employed in industrial gauges (see Figure 
Seo ye 

Fano(41) has derived equations which estimate 
build-up factors for various absorbers. However, the 
assumed system geometry for the published solutions and 
the complexity of the calculations make them essentially 
useless in this application. This has prompted other 


workers to adopt simple empirical models for B such as; 


eeneih Or) 


B=1+0 (pux)®; (42,43) 
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(A) An illustration of broad beam geometry and the 
scattering processes giving rise to build-up. 
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Figure 5.4 The Mechanism and Effect of Build-up. 
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B=1+0(pux) +8 (pux) 7 > (44) ts (5% 14) 


M 
B= 2) A, exp (a, ux) ri (A'3) vieheneo ero 


where: 
A, a, $8 are constants. 


Lambert's Law is thus modified to include build-up as; 


I=BI exp (-pux) COAES*. 1G) 


There are only two means by which build-up can be 
circumvented and both are largely impractical in an 
industrial setting. Firstly, narrow beam geometry will 
virtually eliminate build-up. However, this will also 
effectively reduce the source strength (I,.) and hence 
the measured intensity (I) giving rise to poor 
measurement statistics. (That is, for a Poisson process, 
the relative error of meaSurement is proportional to 
1/VI'*.) Secondly, since Compton scattering predominates 
and because these scattering events are inelastic, the 
incident photons are always in good geometry. Thus a 
detection system capable of resolving primary and 
scattered photon energies could be used to detect only 
uncollided primary photons. Given the resolution versus 
photon energy characteristics of the detectors commonly 


used, this is probably not practical. 


'4For a more complete discussion on the statistics of 
radioisotope emmisions and photon detection the reader is 
referred to any text on basic X-ray analysis(52,67). 
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In calibrating the pulp density gauge, if it is 
judged that build-up creates sufficient non-linearity'® 
then simple models for B, of the form given in equations 
(5.13)-(5.15) must be incorporated in the calibration 


model (eqn. (5.16)). 


5.3.2.3 Broad Beam Geometry 

In the design of most density gauges the source is 
generally collimated to an extent and the detector is 
left essentially uncollimated. Since the crystal of the 
scintillation counter (the most common detector) is much 
larger than the source collimation aperture these gauges 
operate with broad beam geometry. Neglecting the effects 
of build-up and considering the source to be in point 
form, as opposed to disk form'*, the effect of beam 
divergence can be estimated. 

Figure 5.5(A) sets out the geometry of broad beam 
measurement in spherical coordinates, assuming a point 
source. The fraction, 2, of the primary photons which | 
report to a differential element of area on the detector 
surface, with a source-detector separation of x, and 


assuming a vacuum is; 


_x“cos¢dgdé _ cosodode Fe hoe eda) 


erste At 


‘SA particularly good example of the effect of build-up ina 
radio-gauging application is given by Gardner and Ely (see 
(42) pg .7330)ie 

'*In most cases the path length is long enough to make this 
a reasonable assumption. 
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(A) Geometry of the broad beam measurement system. 
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(B) Theoretical attenuation curves for the broad beam. 


Figure 5.5 Broad Beam Geometry and its Effect on 


Calibration. 
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where: 
Q =fraction of photons to the differential element. 
x =Source-detector separation. 
6, » =the angles as defined in Figure 5.5(A). 
Accounting for attenuation, the photon flux on this 


element of the detector surface, AN, is calculated as; 


= eX 
dN=N exp ( Bosscost ben 5.48) 


where: 

AN =photon flux to the differential element. 

No =Source activity in photons/s. 
Integrating over the positive quadrant of the detector, 
i.e. @>0, @>0, and invoking goemetrical similarity 
yields the expression for the flux on the detector, N; 


Les eh 


b 
yl cos¢exp (- ee a ee Garrat) 
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where: 

6,, ¢, define the limits of the detector in the 

positve quadrant. 

By inspection of equation (5.19) it is evident that for 
very slightly divergent beams, i.e. 6,70*, $,70°*, this 
is essentially equivalent to Lambert's law written in 
terms of flux rather than intensity. To illustrate the 
effect of beam divergence on measured flux suppose that 
the detector window is square (ref. Figure 5.5(A), i.e. 
6,=¢,. Choosing a value for 6, and some typical pux 
values, the variation in N with pux can be studied. 
(Note that equation (5.19) requires numerical solution.) 


The results of this exercise are shown on Figure 5.5(B) 


where the calculated flux values have been reduced by 
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dividing by the flux corresponding to the lowest pux 
value used. Clearly, the logarithmic relationship 
expected from Lambert's law is maintained and the latter 
could be applied as the calibration model. In this case, 
the apparent source strength would increase with an 
increase in beam divergence. The resultant increase in 
mean path length would mean a decrease in the slope of 
the attenuation curves (slight increase in the apparent 
uw, aS illustrated in Figure 5.5(B). It can be concluded 
from the sample calculation that broad beam geometry 
will not introduce curvature in the calibration curves, 
within practical limitations. 

Thuse far) buald-up, effects haves been neglected 
although they can be very important in broad beam 
geometry, In this case build-up has the previously 
discussed effect and the calibration model must be 
appropriately modified if necessary. This raises the 
issue of design compromise where economy of source 
requires broad beam irradiation while reduction of 


build-up requires narrow beam geometry. 


5.4 THE OPTIMIZATION OF PULP DENSITY SUB-SYSTEM DESIGN 

A number of authors(44-48) have dealt with the optimal 
design of pulp density gauges in varying degrees of 
complexity. Apart from Tabor and Lakoski, they have 
neglected build-up as well as beam divergence. In all cases 


filtration was excluded from the analysis. Consideration of 
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all of these possible interferences is theoretically 
possible but practically very difficult. 

Since most of the radioisotopes used in this 
‘" application emit Shimondtademaer ic photons, or have 
secondary emissions which are of sufficiently low energy to 
be totally attenuated in the absorber, filtration can be 
reasonably ignored in any analysis. Furthermore, as was 
illustrated in Figure 5.5(B), beam divergence does not cause 
a departure from Lambert's law other than to increase the 
apparent source strength and slightly increase the apparent 
mass absorption coefficient. Therefore, beam divergence will 
permit the use of Lambert's law in any analysis. However, 
Eheneftiects-ofebuildjupycanrnot lbewsoseasilys excluded. 
Equally unfortunate is the fact that it is impossible to 
estimate, a-priori, the coefficients of a build-up model. 
The design analysis which follows illustrates the methods 
uSing Lambert's law and excluding build-up. General comments 
on the expected effect of the latter will be made at the 
appropriate points. 

The objective in optimizing performance is to design in 
such a manner that pulp density accuracy is maximized. In 
addition to good installation this includes maximizing 
sensitivity. The optimization criterion can be obtained by 


application of the Law of Propagation of Errors., 1.e. 


2 
37 ()=[$2] o7 (I) oot 20)) 


where: o(3) =the variance of variable j. 


Taking the square root of both sides of equation (5.20) and 
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dividing by the pulp density yields the expression for the 


relative error in pulp density measurement'’. 


ae) 2 =| o(T) bes bee24) 
The optimization criterion is simply stated as; 
MINIMIZE  jecen soathae2) 


By appropriate choice of the design parameters, which, with 
reference to equation (5.1) include, source strength, mass 
absorption coefficient and path length, optimal design can 
be achieved. Substitution of equation (5.21) into equation 
(5.22) with the appropriate expression for o(I), depending 
upon the nature of the error, allows solution to the 
minimization problem. 

There are several sources of potential error in the 
intensity measurement. Since the importance of each in a 
particular application is the subject of careful assessment, 
they are best treated separately with respect to optimizing 


Gauge design. 


5.4.1 Statistical Error 
This 1S'the error’ which is due to variation in the 
measured radiation intensity at the detector. Since 


radiation intensity measurement is observed to follow a 


'7An equivalent expression can be obtained using the 
classical definition of instrument sensitivity (see (46) for 
example). Note that the relative error is defined as if the 
the intensity error were a random variable. This is not 
necessarily always the case. 
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Poisson distribution, the expression for o(I) is simply /yI. 


The relative error may be written as; 


Celi pd isimgtrie sales) 


QO 


ereleOee oi) 
YT jeux 


and can be minimized by; 


a) 


b) 


choosing a large I, (source strength) subject to 
shielding constraints, (and detector saturation) 

setting pux numerically equal to 2. (This is established 
using calculus for extremum location.) Clearly, the 
radioisotope (u) and path length (x) choices are 
dependent for the given system. 

Kawatra(47) obtained the same results with a slightly 
different approach. His method does not provide any 
information about the objective function in the 
neighbourhood of the optimum, which is available from 
equation (5.23). In fact, the “objective response" is 
relatively flat in the range 1 < pux s 3. The effect of 
build-up is essentially to translate the the objective 
reponse to higher values for the optimum of pux. This is 
expected since, from Figure 5.4(B), the effect of 
build-up is to show an apparent decrease in u. On the 
other hand, the presence of divergent beam geometry 
tends to lower the value of this optimal solution as 
divergence increases. This is also expected since 
increasing divergence effectively increases the mean 


path length through the absorber. 
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5.4.2 Zero Drift Error 

This is the error which is due to the electronics 
employed in the measurement system. While in the strictest 
“sense it is not possible to assign a o(I) value to this 


error source, the effect can be illustrated by setting, 


o(1I)=K RA Boies 20d 
where: K =a constant. 


Solving for the relative error yields; 


o(p) _ —-Kexp(poux) 


p Ip OUX 


owe on 


and minimization is achieved by; 

a) choosing a large Io. 

b) setting pux numerically equal to 1. 
In this case the objective response exhibits pronounced 
curvature in the neighbourhood of the optimum, however, 
excursions of pyux to higher values are more acceptable 
than the converse. The effects of build-up and beam 
divergence can be rationalized along the lines of the 


previous sub-section. 


5.4.3 Proportional Error 

This is the error which is due to such factors as 
radioactive decay, mechanical instability, and variations in 
detector efficiency and amplifier gain. For mathematical 
convenience the standard deviation of the measured intensity 


will be defined as; 
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The expression for relative error in pulp density is then 


- given by; 
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In this case choosing pux as large as is possible will 
minimize the contribution to measurement error. This 
Statement applies to systems including build-up and beam 


divergence as well. 


5.4.4 Extraneous Absorber 

This is the error which is due to solids accumulation 
in the radiation pathway, variable quantities of entrained 
air, etc... These errors can only be minimized through the 
proper considerations in gauge installation and pipeline 


operation. 


5.4.5 More On Entrained Air 

In industrial applications air entrainment is something 
over which the operator may have little control. If 
Significant variations in the air content do not occur, its 
effect can be removed in the calibration exercise (i.e. 
air+water=single phase). If significant variations do occur 
then the pulp density measurement may be virtually 
meaningless. Where high quality density measurements are 


required, such as in the OSCAA system, even small variations 
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in entrained air are unacceptable. Process pulps often 
contain air and in variable quantity. Kawatra(49) has 
observed that in flotation pulps it is not unusual to have 
entrained air contents of from 5% to 20% by volume. 
Fluctuations of this magnitude can be expected to have 
profound effects on the attenuation of the incident photons. 

The effect of air in estimating the pulp composition 
can be easily demonstrated considering the ASHSCAN device 
(see sub-sec. 5.3.1). AS was observed above, increasing air 
content will serve to lower the apparent pulp density. 
However, the mass absorption coefficient will be largely 
unaffected by the presence of air (ref. eqn. (5.9)). Thus, 
in solving the simultaneous linear system, the apparent 
solids concentration will be lower than expected while the 
mineral matter (ash) content will be higher. If the device 
were calibrated about some mean air content then a decrease 
in this variable would have the opposite effect. This is 
precisely the reason why Lyman built a high pressure pumping 
System into the ASHSCAN device. 

One possible but often impractical solution to this 
problem is to make a parallel density measurement on a 
sample of process water "Standard" which is aerated to the 
Same extent. Subsequently, the measured intensity through 
the pulp can be reduced by dividing by the corresponding 
reading on the standard which will largely remove the effect 
of air. Implicit here is that the process water standard 


sample be can be obtained and that changes in the air 


ne a ag 7 di) 


(OW) ers Quet gat yosep. otdatven ah pas alae 
etree: 0 Litres nies) eet oh gl 1 Robsadgte. nt “dtd 
eae we ROS ee a “Yo seansno> abe 


wed of rege “ae oan hws game ay tor anois 


’ iit git, Io TO} tates wif? no balers: 
t fe 7y ra _ 3 | 
' 4 td ant & oat: Na 1h8 be daotta: of 


OA pels ei aay bere tacomst attnae § 


Ow a4 t : a * irk | a any tre in s it ean ah’ f Pe f,2 err 


ireesisp Al . Sa ee a wold oF evise ttiwis 


Pw SOLD Cae aotagiosds 22 6m sid ‘aa 4 a 

Tate Cf in WO eetoeng add ad: pera6¥ 

| —aepeeee aiesited.iaia ott eaivs ee fe 

LW) iO OAT) sis a atbw Sata Riri,” | 

seg ee ser ioum, ameanoD (das) aoddem | 

btsueh' <2 wel: is so vie aN emo sede petattitn 


SAS SS eet | OMe ant $¥ad Biy ow aldatasys doe 


a 
rs 


oh! bs vetoed name yriw HOB 7 ads errr cl 
apo twady WAIBURA out ‘otat noveua i 


eid? ot dobteled tegter orate dene ud oe id ta@og eno <a ~ 


ve Sse wa wow Cre MeISh Ag har ac A eden od zl iidengy © 
LS tat: Peswied. ah ane ie Dov reds TO76W Bas2o0Ig id’ ‘shqaee 
eI! Lora sr son Kl ony etre buppadus ~soetxe: ‘omee’ 
PVSREAGE OTIC wt et ¢ pA LRS' is ee Bepubex ad noo lu etd! 
natn wii oxone: ytvezail ie Hsibe biabees, peh3!M10 putbeet 7 


beabrate 2osa« eeyootg ey Veet iat “nena gtot Kent ovta to : Pa 


i 
726 ws af 2ep eng Bie as erence ee ad, etainne’ 4a 
ek Pas fi! oe 
s a. = : | a. at oe gh. a : ‘a 
a ge? MTL r ru | ) 7 io pie!” lle TS hens, Ga 
“oe | oe “SOREN se a i Ty cc Abe te § 
by : ——-— » i > a ae 


TE 


content occur on a relatively long term basis. Both of these 


requirements are expected to be seldom realized in practice. 


There are a number of more direct means of attempting 


’ to minimize the air problem, ranging from passive in-line 


techniques to rather expensive mechanical methods. They can 


be generally categorized as follows: 


a) 


b) 


Passive Methods 

These methods typically involve a secondary sample 
splitting device (cf. On-stream X-ray fluorescence 
spectrometer systems). Although the primary function of 
such a device is to extract a small but representative 
portion of the primary sample for anaysis, it can be 
designed to provide for fairly quiescent flow thereby 
promoting deaeration. The system will be ina 
"distributed" equilibrium and very small bubbles will 
not have the Same probability of escape as the large 
ones. Therefore, if small bubbles are present in high 
concentration these methods will prove largely 
ineffective. 

Non-Mechanical Centrifugal Methods 

These methods have been employed by the nuclear and HVAC 
industries for the deaeration of process water. In 
essence, a mechanism is used which imparts a rotational 
motion to the fluid thus enhancing air removal at the 
centre by means of centrifugal forces. The nuclear 
industry in particular(50) has developed the technology 


to the extent that very little water is lost in the air 
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take-off stream while very efficient deaeration is 
achieved. 

Mechanical Methods 

These methods involve mechanical systems which operate 
either on the principles of (b) as, for example, 
Autometrics PSM air eliminator(51) or on pressurization 
to collapse the air bubbles (cf. ASHSCAN). 

The capital and operating costs of the class (c) 
methods make them a last resort proposition in the 
opinion of this author. Class (b) techniques would 
require reasonably careful design and operational 
control to avoid sample bias but, if used in conjunction 
with class (a) methods would seem to provide an 
attractive alternative in both a technical and economic 
sense. The class (a) methods represent the simplest 
possible solution and a variant was chosen for use in 
the OSCAA prototype. It is implicitly assumed that the 
slurry flow streams would not contain high 


concentrations of very fine air bubbles. 


5.4.6 Summary 


In the optimal design of nuclear pulp density gauges 


the souce strength should be maximized subject to the the 


shielding contraints. Depending upon the principal source of 


error(s) the optimal value of pux can take on different 


values. In practice the gauge is usually designed such that; 
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Pes ux 3 
Since p is the property to be measured, the radioisotope (xu) 
and the path length (x) must be chosen to satisfy the 
equation (5.28). Special care must be taken to minimize the 
entrained air (especially fluctuations therein) in the pulp, 


wherever this 1s possible. 


5.5 THE PRINCIPLES OF BACKSCATTER AND FLUORESCENCE RADIATION 
MEASUREMENT 
In this OSCAA sub-system the interaction of the 

incident radiation with the sample is somewhat more complex 
than in the transmission measurement, although Lambert's law 
is the theoretical base. Whereas in transmission gauging, 
one is simply concerned with the overall photon attenuation 
in passage, the backscatter geometry requires more detailed 
knowledge. That is, at any particular level in the sample it 
is necessary to know: 

a) the incident photon intensity, 

b) the fraction of these photons absorbed in this layer, 

c) the fraction of these absorbed photons which interact 
via the various processes and the spatial distribution 
of the secondary photons which result, 

d) the attenuation of the secondary radiation which is 
directed toward the sample surface, 

e) the sum of all of these factors over all layers in the 
sample. 


While most textbooks on X-ray fluorescence provide some 


ati im oa ee ae Jean tT: > sabsaae 2) 


ay 


ay. ae ai ls 9 


% 


ea 


x | A * 


shuts 9) eka na PA: sa "0. exsatoutas werd 


&4niee ) asayendae rv O20 aids’ an 


ravod OF aaa oat ni ‘aphtetbe dan 


| 
rope ; Fs 


Hane su om noige imanss3 ang ad) 


Vea Le 
ct 1 a ine a eee ig @ site anv a ae aed i ut dao ee 7 
‘ ne ’ } TAR oh kee = 2 

roid? | Sinai oo oft) a3 a ‘Bow sarin tigate ate 
K 


Loeb 4 2¢m 2 osha pasieretoed sanz ies: 


Loline eels ae. tee aly iglesia? 4 von te va) ‘$edT -spbeLng 
oma teond: ‘oF yiaeneen aE 


a 


‘Given idous taebion a3 tal 


weuei ett gi an iio Deilig seats ig eohtoas’ ae ro 


“Si 


S2Gq2osat- haw 2 otetic padscede omand Fo aoisoet? ., 


i ‘ 
/ 


iHpe jee? Thee. (2k deol atd Bene tmeaee 1g BOE TRY, act, sky, 
uae ahs Pee ie waiaot e abnooee ott to" 
‘pl Pride noize se IRR An wee ont te no ttaunesss att 


\ 


sjn3 “tue, sane wa $aawo: beioertb , ose 


#)i2 yO Se me | ey: ‘le 2 ave 2" 4. IDBt spade: We ite ye mee ois. {5 ‘ ‘ , 
. ; ‘= : ” 
ial a}: ey |  -siqmen “4 


‘ 11 i | 
segr abi ving «.nesassovii pens? h | eaoaseet Seon os tite Ra 
ne 
rs a ve 
Piri " 


- 


re iin and ; , 
oe / 


ree) |. na ery etl en ee 


ee ne o 


114 


theoretical background to the subject (see (52); Chap. 4 for 
example), the mathematical modelling efforts of Rhodes and 
co-workers (53,54,55)' are moré ‘easily applied to the present 
_ investigation. Although in theory it is possible to account 
for a continuous particle size distribution(54,55), in 
practice, and for the purposes of OSCAA, it is sufficient to 
lump the size distribution into a single average particle 
diameter(53). The mathematical development which follows is 


a consolidation of the theory developed by Rhodes et. al. 


5.5.1 The Theory of the OSCAA Ash Analysis Sub-System 
In developing the theory of operation of similar gauges 

Rhodes et. al. made the following assumptions; 

a) The primary and secondary radiation enter and leave the 
Sample perpendicular to the sample surface. For oblique 
incidence the terms uw; in the equations which follow are 
merely replaced by u;cosecé, where 6 is the angle the 
incident beam makes with the surface. | 

b) All particles have the same average size with an average 


linear dimension d, where d is defined as; 


~ stars 2s) 
d= Ve 
where: 
V =particle volume. 
a =particle area presented to the incident radiation 
averaged over all possible orientations of the 
particle. 


c) In a sample containing more than one type of particle, 
the distribution of the particles throughout the sample 


is completely random but, in any horizontal layer of 
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thickness d, the average relative volume distribution of 
the particles is the same as the average over the whole 
sample. In Figure 5.6(A) the (i+1)th layer represents 
such a typical galore? in the sample. Moreover, the 
physical and chemical properties (p,u, etc...) of each 
individual particle are homogeneous. 

Under the assumptions above, the radiation flux leaving the 


sample from the (it+1)th layer, N,.1, is; 


Niyz (F-11955; adn, F. eo ehoeee On) 


wherest Njareh) sradiat ton et lix of type st rarriving vate«the 

Sample surface from the i+1th layer; t= f for 
fluorescent radiation and t=b for backscattered 
radiation. 

Ipo' =incident intensity of the primary eh Sei. 

7; =number of j type particles in the it+ith layer 
‘where there are a total of k different particle 
types. 

a =average particle area presented to the incident 
radiation. 

J; =average transmission through i layers of the 
Sample for the incident radiation. 

J,;' =average transmission through i layers of the 
Sample for the secondary radiation. 

F,; =secondary radiation intensity per unit incident 
intensity for particle type j. 


(Implicit in this formulation is that the 
fluorescent and backscattered radiation are 
mono-energetic. This is correct for the former when 
considering a particular line, however, since both 
inelastic and elastic scattering can occur this 
assumption is in question for the latter. It will be 
shown later that, for the OSCAA system, this 
assumption is indeed reasonable.) 


To derive the fluorescent and backscattered radiation models 
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from equation (5.30), the terms on the right hand side, 
other than Io, are treated individually. To further simplify 
the analysis each type of the secondary radiation will be 


considered separately. 


5.5.1.1 Fluorescent Radiation Model 


Reaemice. on Factors2nd 4 sed 4 

Consider J, first of all. For generality, suppose 
that the sample consists of k types of particles. For a 
very narrow photon beam the i layers may be divided into 


the set ingens fen..en eg) where; 


k 
Seen, Fi SepetO ou) 
J 
where: 
n; =the number of layers in which the incident 
photons encounters the jth type of particle. 
Invoking the multinomial distribution, the probability 


that the radiation beam will intercept any specific set, 


n, 1S given by; 


k 5 i N05 32) 


where: 

m; =the probability of encountering the jth particle 
type in the sample. This will be taken as the 
volume fraction of this particular species in 
the entire sample, C;. . 


The expected attenuation of the incident photons is then 


calculated as; 
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M;,o; =the mass absorption coefficient and density 
for the jth particle type. 


Employing the following change of variable; 


@x=1 Le 2d 30.. ete 34) 
7 oot rae, 


and the multinomial series identity, the expression for 


J; can be written as; 


L i=n, 
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Deny. Be, (155035) 
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or, with the appropriate substitution; 
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Since the fractional attenuation across a layer is 


constant. By similarity; 
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where: 
uj;' =mass absorption coefficient for the fluorescent 
radiation in the jth particle type. 


Number of Particles: 7, 
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The number of particles of type j, n;, is 


calculated as 


nj=C,Ad =C,a/; ORES 


V 
where: 
C; =volume concentration of the jth particle type in 
the sample. 
A =cross sectional area of the irradiated sample. 


Intensity Factors F; 


With reference to Figure 5.6(B), consider the 
infinitesmal layer dx, located at a depth x below the 
surface of the jth particle type. The fraction of the 


incident intensity absorbed by this layer is; 


d 
=-— I “11,0.-xX)axt=I. UW. Oe -U.0.x)ax 
STF ES au sls td ot 505? 505 BERS 039) 


where: 
Pi; =intensity collided in the layer x to x+dx for 
particle type j. 
Lipa pposacteibutes of jth particle. 


The probability that these photons will produce 
fluorescent radiation is calculated as; 


Tov. W 
be er ere ... (5.40) 


j uy 


where: 

P2; =probability that the absorbed photons will 
produce fluorescent radiation. 

7; =photolectric mass absorption coefficient for 
particle type j. 

y; =fluorescent yield for particle type j. 

Wx; =Mmass fraction of the X-ray fluorescent element 
of interest in the jth particle. 


In particular, 
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Figure 5.6 Conceptual Illustrations of the Slurry Sample and 


an Individual Particle. 
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from the various shells (K,L and M). In most 
applications one is looking specifically for K or L 
Shell X-rays in which case the appropriate values for 
the shell of interest would be input to equation (5.40). 
The attenuation of the fluorescent radiation which is 


directed toward the particle surface is given by; 


where: 
P3,; =fraction of fluorescent photons reaching the 
particle surface from a depth x. 
G =geometrical factor. 
Thus, the total expected fluorescent radiation intensity 
at the particle surface is the product of equations 


(5239), (5,40) fend 35.41) antegrated over the range ot 
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wheres: 
I,; =intensity of fluorescent radiation at the 


surface of a j type particle. 
Dividing both sides of equation (5.43) by Io gives the 


expression for F; as; 
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(5.44) into equation (5.30) and summing over all z 
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layers, where d (z+1) is the sample thickness, yields; 


NoC.T YW. - ie Zz i 
N(£) = —211 6(1-exp(-o.d(uitu.'))  (o.d.") 22.05.45) 
VE et Seer a 


where: 
No =incident flux upon the sample, i.e. No =IdA. 


Of particular importance in this work is the case where 
zrm, i.e. an effectively infinitely thick sample. Using 
the convergence properties of infinite series (see 
(56);pp. 413-415), the transmission factor product 
Summation can be written as; 
Set ge ee meib. £6)) 
00 (1-gJ') 
and the final expression for fluorescence flux at the 


Sample surface is; 
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This model can be further simplified by allowing 
the average particle size to approach zero. Under this 
assumption the sample becomes a homogeneous mixture. 
Taking the limit of equation (5.47) as d+0* gives an 
indefinite result as both the denominator and the 


numerator vanish. Applying L'Hospital's rule yields; 
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where: 
q =a-constant. 
w; =mass fraction of component j. 


and substitution gives; 
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Rhodes and Hunter(55) have done a similar exercise in 
reducing the more general form of the model given by 
equation (5.45) to the homogeneous sample 
representation. The difference lies in the fact that 
equation (5.50) assumes an infinitely thick sample. It 
will be shown later that this is a reasonable assumption 
with the OSCAA flow cell design. 

Figure 5.7(A) compares equations (5.50), the 
horizontal line, and (5.47), the curve, as a function of 
the mean particle diameter, d. The ordinate data is the 


predicted reduced iron fluorescence flux typical of the 
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Figure 5.7 Particle Size Effects in Fluorescent and 


Backscattered Radiation Measurement(Flotation Feed). 
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various flow streams around the McIntyre flotation 
circuit (ref; Figure 3.1)). (Only the curves for the 
flotation feed are shown, but they are representative of 
the tailings and concentrate streams as well.) The 
vertical line at 212 microns represents a typical mean 
particle diameter derived from screen analysis. For the 
purposes of this calculation the irradiating isotope was 
taken to be Cd-109 and the fluorescent radiation of 
interest was the iron Ka X-rays. Since the detailed 
mineralogy of the solids is unknown a "best case" and 
"worst case" are presented. The best case assumes that 
the slurry consists of water and solids particles. The 
latter are assumed to be homogeneous with constant 
composition. The worst case assumes that the slurry 
consists of water and three types of solids particles. 
These are iron-free coal organics, iron-free mineral 
matter, and hematite. The method of estimating the 
component attributes is discussed in detail in Chapter 
he 

The) results “illustrated on Figure 5.7(A) compare 
with similar calculations done by .Rhodes and Hunter (see 
(54); fig. 3). For a particular flow stream there is 
considerable difference between the best case and the 
worst case. As expected, the best case curve is much 
less sensitive to particle size. From a practical point 
of view the "real" curve will lie somewhere in between 


the limiting cases, however, it is expected to lie 
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closer to the best case curve because of the partially 
liberated nature of solid phases. This will be 
especially true for the flotation output streams where, 
by the selective nature of the process, some 
homogenization occurs. The slope of the "true" curve at 
the "true" mean diameter will determine the extent to 
which changes in particle size will affect the 
fluorescent radiation flux. It is expected that this 
will increase in importance in the order of concentrate, 
to feed, to tailings. The other major point of interest 
from this figure is the loss in secondary photon flux as 
d increases. This would be expected to have implications 
for design withprespect (to, the. determination of 
accumulation times. 

Despite the differences in predicted photon flux 
between the homogeneous sample and the best and worst 
cases, equation (5.50) is proposed for calibrating the 
OSCAA analyser. Apart from the absence of explicit size 
dependence (for which there is no means of meaSurement), 
the only difference is that the apparent incident photon 


flux will be lower. 


5.5.1.2 Backscattered Radiation Model 


Transmission Factors 
These are of the same form given in equations 
(5.36) and (5.37) with the only difference arising in 


equation (5.37) where, u;', is the mass absorption 
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coefficient for the backscattered radiation. 
Number of Particles 

This is identical to equation (5.38) given above. 
Intensity Factors | 

Unlike fluorescence, where only one element (iron) 
in the sample was considered, each component in each 
particle will scatter some of the incident radiation. 

Proceeding as in the previous sub-section with 
reference to Figure 5.6(B), consider the infinitesmal 
layer dx, at a depth of x in the jth particle type. The 
fraction of the incident radiation absorbed in this 
layer is given by equation (5.39). The probability that 
these absorbed photon will produce backscattered 


radiation is calculated as; 


Boge ee oly 
where: 

P2; =probability that the absorbed photons will 
produce backscattered radiation in the general 
direction of the detector. 

0; =a composite scattering mass absorption 
coefficient which has contributions from both 
the Rayleigh and Compton processes. It depends 
upon incident photon energy, absorber character, 
system geometry, etc... 

The attenuation of the backscattered radiation is given 
by equation (5.42), with the understanding that u;' is 
now the mass absorption coefficient for the 
backscattered radiation. 

Thus, the expected backscattered radiation at the 
surface of a j type particle is the product of equations 


(5.39), (5.42) and (5.51), integrated over the range of 
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Ciccone) ) J 3 

a) 
where: 
I,; =intensity of backscatter radiation at the 


surface of a j type particle. 
Division of equation (5.52) by Io gives the intensity 
RPactCu ras) 
52 a 
ne esac) 


j na: 
(Hs Hs ) 


Substitution of the appropriate equations into equation 
(5.30) and assuming an infinitely thick sample yields 
the model form of equation (5.54); 
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Before deriving the homogeneous sample model 
another simplification can be made to equation (5.54). 
Rayleigh events produce scattered photons with the same 
energy as the incident photons whereas Compton events 
result in an energy loss. The extent of this energy loss 
depends upon the scattering angle, w, and the incident 


photon energy. Mathematically; 
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AE=E ue eee 5 e 55 
| arene a 


where: 
AE =energy difference between the incident and 
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Scattered photon in kev. 
yw =the scattering angle; 0° < w < 180°. 
For Cd-109, with E~22 kev and using w=120° and w=180°, 
the respective energy losses are 1.3 kev and 1.7 kev. 
With such relatively small changes in photon energy one 
can assume tnat the mass absorption coefficient for the 
primary and secondary radiation are equivalent, i.e; 
No Uri rei Corer On 


and equation (5.54) may be rewritten as; 


NgG Ke o7.Gs i 
N(b) =————._ 9) —-1_-1(1-exp(-29.du.)) ...(5.57) 
(1l-gg') 2u. Jd 


Now, taking the limit as d+0*, employing L'Hospital's 


rule gives the model for a homogeneous sample as; 
Ms BAe), 


Again, as with equation (5.50), equation (5.58) 
implicitly assumes an infinitely thick sample which is 
reasonable for the OSCAA system. 

Figure 5.7(B) compares equation (5.58), the 
horizontal line, and equation (5.57), the curve, to 
illustrate the particle size effect in backscatter 
measurement. The ordinate data is the predicted reduced 
backscatter photon flux. All other conditions are the 
Same as those used in the development of Figure 5.7(A). 
(For the purposes of clarity, only the flotation feed 
curves are presented.) The scattering cross sections are 
computed by the method described in the following 


Subsection. Again, best case and worst case curves are 
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presented. 

The results illustrated on Figure 5.7(B) indicate 
that the particle size effect is much less pronounced in 
backscattering than in fluorescence. In particular, the 
best case curves are very nearly the same as the 
homogeneous sample model for all three flow streams. It 
is expected that equation (5.58) would be useful in 


calibrating an OSCAA type system. 


5.5.1.3 Calculation of the Composite Scattering 
Coefficient 

The composite scattering coefficient, aj, 
introduced in equation (5.51) is not simply the sum of 
the tabulated Rayleigh and Compton coefficients'*® as was 
assumed by Rhodes et. al. in their work. Both of these 
events exhibit a strong angular dependence and to avoid 
bias in the composite coefficient it is necessary to 
appropriately correct om, and O,m. The method amounts to 
Calculatingtfirst order correction factors, rm, and itan 
which, when applied to the tabulated values, provide the 


composite coefficient, o, i.e. 


Oe EC a on nim 2+. (5.59) 
where: 
o =composite scattering mass absorption coefficient. 
Omr Onm =COMpton and Rayleigh mass absorption 
coefficients. 
Linnie Ui COGGec taonmtactons:. 


For the OSCAA system geometry (see Chap. 6) it is 
‘®The tabulated values account for scattering in all 
possible directions. 
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required to calculate the correction factors for 
scattering in the range of 120° < w < 180°. The actual 
intersectionson thegdetector with this band is) accounted 
for by a simple factor reflecting the area covered and 
it can be lumped into G. 

For Compton scattering, the exercise begins with 


the Klein-Nishina equation given by Zeigler et. al.(57). 


zZ 
a -26 E! EB Ee F By, weet 5S. 60) 
oe =3.92X10 E IE + pee sin | 


where: 

Om’ =the Compton cross-section per electron per unit 
solid angle, i.e. cm* per electron per 
steradian. 

E =the energy of the incident photons. 

E' =the energy of the scattered photons which is 
calculated from equation (5.55) as E'=E-AE. 

yw =the scattering angle. 


To compute the overall cross-section it is necessary to 
integrate equation (5.60) over the spherical surface, 


and muliply by the number of electrons per atom, 27; 


Ze 
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Davisson and Evans(58) provide analytical expressions 
which allow this calculation, however, numerical methods 
provide equivalent answers and are simpler to use. The 


correction factor is computed as; 


by 
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For Cd-109 with E=22 kev, numerical solution of equation 
C563) yields ir.=09279: 
For Rayleigh scattering, Zeigler et. al. have shown 


that the cross-section is approximately given by; 


mri 3 
Sam “KEW 27 23 Oa) 


where: 
Onm' =the Rayleigh cross-section per atom per unit 
solid angle, i.e. cm* per atom per steradian. 
k =a constant. 
f(W) =a function dependent upon the scattering 
angle. 


The function f£{))is shown inePigure 5.8. Using a 
digitized representation of this curve in the numerical 


solution of equation (5.65), 
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For. Cad-109Myreldsir: 7. =0). 053) 
In both equations (5.63) and (5.65) the 
denominators representing the total scattering 


cross-sections were compared to the tabulated values 
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Figure 5.8 The Angular Dependence of the Rayleigh Scattering 


Cross-Section(Ziegler(57)). 
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with good agreement. (Before leaving this, and for 
completeness sake, the mass absorption coefficients are 
computed from the cross sections above by multiplying by 
Avogadro's number and subsequently dividing the result 


by the atomic weight.) 


5.6 THEORETICAL IMPLICATIONS IN THE OSCAA SYSTEM DESIGN 

To ensure optimal operation of the OSCAA prototype it 
was important to incorporate the theory in the system 
design. This minimizes calibration problems to the extent 
that are dependent upon design. 

For the nuclear pulp density gauge, optimal design 
procedures were outlined in section 5.4. Provided that one 
has a knowledge of the major potential error sources the 
design is generally governed by equation (5.28). While this 
implies that any radioisotope may be used provided that the 
path length is suitably chosen, there may be practical 
limitations on the latter. A compromise must therefore be 
reached. In a recent publication, Kawatra(59) has observed 
that Gd-153 is better suited for pulp density measurement in 
coal slurries than Cs-137, which is most commonly employed. 
Unfortunately, he does not elaborate on the reasons 
(presumed to be geometric) for this, however, the statement 
does support the choice of Gd-153 for the OSCAA prototype. 

For the ash-analysis sub-system, the optimal design 
procedures are not nearly so straight forward. Since the 


technology of flow cell design has been developed and 
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refined in on-stream XRF applications, this design was 
direclty adopted. Futhermore, since it was desired to detect 
both iron Ka fluorescent and backscattered radiation, as 
well as to intersect a large fraction of this secondary 
radiation a wide-window proportional counter was selected 
for detection'’. Xenon was chosen as the detector gas to 
ensure a high probability of detection for both types of 
secondary photons. The detector window was made of very thin 
Beryllium to minimize the attenuation of iron X-rays. This 
leaves the radioisotope as the principal means of optimizing 
design. 

Considering the work of others employing backscatter 
geometry (see 4.4.2) and the physical constraints of 
fabricating the source and detector assembly, a single 
radioisotope was chosen for irradiation purposes. As 
previously mentioned Cd-109 was selected and the rationale 
for this choice is discussed below. 

A cursory analysis indicated that the best two 
candidate radioisotopes for this application were Cd-109 
(ref. Kawatra(36)) and Pu-238 (15<E<17 kev). (ref. 
SORTEX(11)). Subsequent investigation of their availability 
indicated that Cd-109 was most easily obtained. (In fact, no 
source for Pu-238 was found.) Since Cd-109 had previosly 


been used with reported success and because of the 


'*There was no effort made to evacuate the space between the 
source and detector assembly and the flow cell window. Under 
these conditions and with XRF style irradiation and 
detection equipment the possibility for analysing lighter 
elements, e.g. 165ZS525, would exist (see(60) for example). 
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availability factor it was chosen for use in the OSCAA 

prototype. 

Retrospectively, it 1S possible to look at the 
implications of having made this selection. Boyce et. 
al.(61) provide the general framework for such an analysis 
in examining the design of the SORTEX analyser. To simplify 
the illustrative calculations which follow, several 
assumptions are made: 

a) The homogeneous mixture models for an infinitely thick 
Sample can be applied. 

b) The slurry consists of three components, water, coal 
organics (assumed to be Carbon: Z=6), and mineral matter 
(assumed to be Magnesium: Z=12). All of the iron is 
assumed to be in the mineral matter. 

c) The average attributes of the flow streams are given in 
Figure 3.1 and the mean slurry pulp density is 1.05 
g/cm? in all cases. 

In the design of an OSCAA system there are three 
criteria of major importance: 

a) Maximize sensitivity, 

b) Maximize flux (minimize accumulation time for a given 
accuracy), 

c) Ensure that sample is representative. 


In mathematical terms, instrument sensitivity is defined as; 


Spon 
N (oa Bete 56) 


where: S =sensitivity. 53 
N =flux of secondary photons at the sample surface. 
a =attribute to be measured. 
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The models may be written in the form, 


NiB)= (i=s)6 +s{l-m)o.tsmo 
as is a OOT)) 
(1 s)u+s(1 m)u+smu, 


where: w,0,m =Subscripts denoting water, coal organics, and 
mineral matter respectively. 
Ss =mass fraction of solids in the sample. 
m =mass fraction of mineral matter in the solids. 
(For the purposes of these calculations the composite 
scattering coefficients are simply the sum of the tabulated 


Rayleigh and Compton coefficients.) 
tTsm 


A Cis i pe ore a mgerneT Te ee, Sane Ra 
(1-s) +s (1-m)uotsmu Regn Gsi-# cis.) 


where: @, =Sum of the mass absorption coefficients for the 
incident and iron Ka photons, i.e. #,;=u; +u,'. 
The expression for sensitivity as a function of mineral 
matter content are: 
SS m(bc-ad) 
(c+md) (a+bm) 
with: = las +a¢ = bs Barra sy) ae) 
a=(l S)o. SO, oc eae S)U +su 


b=s (5-707) d=s (uy -u) 


and, 
(Tes) rate 


Ss(£) = —me mm — 
(1-s)u, +s (1-m)uo+smu_, 4s Ace) 
With the exception of the iron Ka radiation, the mass 
absorption coefficients are dependent upon incident photon 
energy and both flux and sensitivity for each secondary 
radiation must be calculated. Figures 5.9 and 5.10 
illustrate this energy dependence for both backscattered and 


iron fluorescence radiation respectively. For sensitivity 


the curves are plotted at extreme values of the mineral 
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BACKSCATTERED RADIATION 
SENSITIVITY 
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Figure 5.9 Sensitivity and Flux in Backscattered Radiation 


Measurement. 
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IRON FLUORESCENCE RADIATION 
SENSITIVITY 


0) 10 20 30 40 50) 
(A) ‘ PHOTON ENERGY (kev) 


REDUCED IRON FLUORESCENCE 
ELUX 
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(B) PHOTON ENERGY (kev) 


Pigureé 5,10 Sensitivity ander.ux in Lronarluorescence 


Radiation MeaSurement. 
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matter concentration, m. The results are consistent with the 
findings of Boyce et. al. and since the general trends are 
the same for all flow streams, only the feed stream data is 
' presented. Figure 5.11 is included to demonstrate the 
relative difficulty in making the measurement on a Slurry 
Sample as compared to dry solids. The vertical dashed line 
represents the Cd-109 photon energy. Pu-238, having 
principal emissions in the range of 15-17 kev, lies just to 
the left of the Cd-109 line. 

On the basis of Figures 5.9 and 5.10 the following 
conclusions can be made. For backscattered radiation, Cd-109 
Sensitivity is somewhat less than that for Pu-238, however 
the-photon flux would be significantly greater, all other 
things being equal. For iron fluorescent radiation, 
sensitivity is about the same for both radioisotopes but 
Pu-238 would provide a significanlty greater flux. In 
summary, Pu-238 would appear to be a marginally better 
choice for a primary photon source than Cd-109. 

With respect to ensuring a representative sample, this 
is simply measured by the maximum depth from which~.secondary 
radiation will return to the sample surface, i.e. the "thin 
Slice" for analysis. Since sample homogeneity has been 
assumed, the appropriate expression for calculating this 
depth is taken from the decay term in equations (5.44) and 
(5.52). Using Clayton and Cameron's(46) definition of 
Saturation depth as the "depth from which 95% of the total 


secondary radiation is obtained", the appropriate 
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Figure 5.11 Backscatter Sensitivity Differences in Dry and 


Slurried Solids Samples. 
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expressions are; 


SENS i aa Mila OAM BS aE Re el es a) 
2o{(1-s)u +s (1-m)uotsmu, } 


rea Seats) etn ee eo ee eo | 
p{ (1-s) +s (1-m)ii,tsmil,} eer tou 72) 


where: d,(sat.) =the saturation depth in cm. for secondary 
radiation type j. 


Figure 5.12 illustrates the relationship between saturation 
depth and photon energy for both types of secondary 
radiation. From this figure it is evident that a 
Significantly thicker slice of the sample is analysed using 
Cd-109 for the backscattered radiation. This is true in the 
case of iron fluorescence as well but to a lesser extent. In 
this case Cd-109 would be preferred to Pu-238 because of its 
Sampling characteristics. In considering the saturation 
depths for the secondary radiation in comparison to the flow 
cell depth of 2.54 cm, the infinitely thick sample 
assumption in the modelling work is justified. 

This semi-quantitative analysis indicates, at least to 
the author, that there is little to choose between the two 
radioisotopes. In fact, the final choice would logically be 
made on the basis of price, availability, half-life, etc..., 


with the appropriate weighting factors. 
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Figure 5.12 Saturation Depth for Secondary Radiation. 
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5.7 CALIBRATION EQUATIONS 

Two types of calibration equations may be used with an 
OSCAA type system, either a set of semi-empirical equations 
-derived from the theory or a set of empirical equations. For 
the purposes of process monitoring the most important stream 
attributes are the solids content of the slurry (% by 
weight), and the mineral matter content of the solids. The 
iron content, while important in compensation, does not 
require to be reported, unless it is on an alarm basis (e.g. 


tramp heavy medium). 


5.7.1 Semi-Empirical Calibration Equations 

The physical model of the slurry employed in 
calibration includes water, iron free coal organics, iron 
free mineral matter and hematite in an infinitely thick 
homogeneous mixture. Designating the mass fractions as w,, 
Wo, W3, and wy respectively,the necessary models are written 


in the form; 
Pulp Density 


Cp 2 
bom. | ow] = A ee eee LOCKS +BY 


Yow 


L/P; reer) 


where: Cp =counts through pulp in the accumulation time, 
Cw =counts through water in an equivalent time 
period. 
w; =mass fraction of component i in the slurry. 
p; =density of the ith component in g/cm?*. 
A,, B; =calibration constants where, in theory, 
A,=7"UXx and B,=P;,uUx. 
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SPenie ee ae) 


where: Cb =backscatter counts in the accumulation time. 

a, =the composite (see eqn. 5.59) scattering 
coefficient. 

Az =calibration constant including incident flux, 
geometry, etc.. 

Bz =calibration constant accounting for offset due 
to detector choking, radiation leakage, 
backscatter from flow cell housing, etc... 


Iron Fluorescence 


uw, e575) 


where: Cf =fluorescence counts in accumulation time. 
AgveGimi lar: to-A.>. 
B; =Similar to B:. 
From calibration experiments each of the three models can be 
parameterized (A,,B,) using linear regression. The component 
attributes p;, u;, and, o;, are obtained by experiment, by 
curve fitting or from published tables. 

Subsequent to the individual model calibration 
exercise, the three equations can be combined to forma 
linear system which may be used to predict the mass 
composition of the slurry from the radiation measurements. 
In vector-matrix notation the system looks like; 


2a ber ie aa i 


Pa 


Ww =| D(us-H4) - (05-04) D(U-H4)- (03-05) D(u,-H,)-(9,-05) 
ro ae ek ie eG) 
Wy (HoH) (Uz-H) (U,-Hy) = = 
q. 
ee 
r) 5, -Du, 
and, 


mihi a od i pa eR eeetoad= dd 
‘oP ec.¢, .pS sea Ms begs ene oO. 
PASTEL I2 995 

e smn: ivide id! (6a* cf 7 ) 

, - OUS. (eRFowg s¢ yt 
Nop! ieteies Horses Lass @ ) 

sAgel bat. : » s6RH0gS6 .od 
: | Dts of}. moi? tedgenrzoed 


, 
Wk x 


sane okasou lie YD 

» Z ThA Nags af 

40 OF talintia= 28 

P ‘ BR ‘ e ia 

zfenur 2¢3¢4 efs {Lo doBe winemiiegze noles3ci tao wo 


& np 
enil oalew (,8,,A) beritsdsmaml 


do 2%) 0 ,bne, a 1.0 eoUeioae 


i Ce nt ood 
-, seload bacwhiduq meat ae gaisat yan 


i 7 > peor = % { 

) Labor Lev Sawin! edt os tndupeddta oy 

Dos 6090 ad poy engi gee Baste oft ,saiotemee 
; hehe a a ~ 

iy 29theag ee dee id fet doidw wegeye wenkly 

; ‘ J ) ry ed ; ae) 

aii nor? yrieie exit Yo nod tieequoay: 


- 
7h! 


> 
~~ 
~ 
a4 
’ 


We ea pre 
9 (8Ng. wie Os patom Xizism-sotoev. At y 


r is n 


Ge «Gzen 


145 


i PALS a bale = a ee UG. 27a} 

pai 

Cb- 
speed EOL S) 

Re 

f- 
a za eee Ot oO) 

Be 


whence; 


S=1003% (w. +w +w 


37W 4) e2e5. 80) 
where: s =% solids in the slurry by weight. 


and, 
means [ “34 | 


Ss (5.81) 
where: m =% mineral matter in the solids by weight. 
and, 
eg 
ae ok gai 2 82) 


where: a =% ash in the solids by weight. 
Rma =mineral matter to ash ratio for the coal. 


5.7.2 Empirical Calibration Equations 

Cooper(62) has observed that, in general, calibration 
of on-stream analysers using semi-empirical models provides 
better predictive accuracy than the more frequently 
encountered general polynomial (empirical) calibration 
equations. Since the latter was entertained in the OSCAA 


protype calibration work the formulation of Smallbone and 
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Gurvich(60) is presented here; 
nN nN n 


where: Cj; =the counts of the jth radiation type in the 
accumulation time. 
n =the total number of radiation measurements made. 
a jie, ay. (= sbest £1 Ei ecoeriicrents. 
P =the slurry property to be predicted. 
Implicit in this formulation is that only those terms which 
are statistically significant would be retained in the 
calibration equation. Thus, stepwise multiple linear 


regression would be a suitable method of parameterizing the 


equation. 


5.8 SUMMARY 

The theory of operation of the OSCAA sub-systems has 
been discussed with emphasis on potential sources of error 
and complications in the measurements. The theory has been 
used in some aspects of the prototype design as well as in 
the development of calibration models. The content of this 
chapter should equip the reader to understand the operation 


of any similar type of on-stream coal ash analyser. 
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6. OSCAA LABORATORY EQUIPMENT DESIGN AND EXPERIMENTAL 


PROCEDURE 


6.1 INTRODUCTION 

The major effort in this thesis project was to design a 
prototype on-stream ash analyser and study its performance 
in the laboratory. Based upon the results of this portion of 
the project a limited plant trial was to be considered. 

It is generally observed that equipment of this type 
will perform better under well controlled laboratory 
conditions. The measures of performance made in the 
laboratory are estimates of what can be expected in the 
field. By attempting to closely simulate "real-world" 
operating conditions the estimates so obtained are more 
reasonable expectations of plant performance. 

For the purposes of the laboratory investigation four 
calibration runs were deemed necessary for performance 
evaluation. These were: | 
a) Calibration on Flotation Feed, 

b) Calibration on Flotation Concentrate, 

c) Calibration on Flotation Tailings, 

d) Calibration on Flotation Feed With Varying Iron 
Concentration. 

The sections which follow describe the final equipment 

configuration, the experimental design and the pre/post 

calibration run sample handling methods. The discussion of 


the results of the calibration tests are treated separately 
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in Chapter 7. 


6.2 EQUIPMENT DESIGN 

The OSCAA equipment layout for the laboratory program 
is illustrated in Figure 6.1. For the purposes of a more 
detailed explanation it is convenient to consider the system 
in three parts: 
a) ‘The Flow Loop, 
b) The Pulp Density Sub-System, 
c) The Ash Analysis Sub-System. 


These are discussed below. 


6.2.1 The Flow Loop 

A SALA model SPV 180 3.8 cm (1.5 in.) vertical pump 
with a 0.028 m® (1 £t?) sump and a 5 horsepower variable 
speed drive was used to pump the slurry to the secondary 
sampling device. The primary ne eee ery rate was 
nominally 9.1 m?/h (40 USGPM), which is typical of most 
on-stream XRF systems. The primary sample slurry flow rate 
was estimated from a calibrated "V-notch" weir in the 
collection launder which delivered the recycled pulp to the 
sump. The launder also provided for access to the flow 
streams for sampling purposes. 

The pump sheave size was calculated for a primary 
sample flow rate of 13.6 m?/h (60 USGPM) of water, 
recognizing that the corresponding delivery of aerated 


Slurry would be somewhat less (ref. (63)). The head for this 
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Figure 6.1 A Schematic Illustration of the OSCAA Prototype 


Test Rig. 
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flow rate was estimated by solving the macroscopic 

isothermal energy balance and these two values were used in 

conjunction with the manufacturer's pump curves to specify 
impeller speed and drive horsepower. A variable speed 
controller was included on the drive for fine adjustment of 
flow rate. 

The function of the secondary sampling device in the 
flow loop is threefold: 

a) To split a representative secondary sample from the 
primary sample input. This sample was to be provided to 
the OSCAA components at a substantially constant static 
head to prevent window warping in the ash analysis 
sub-system, which would affect measured radiation 
intensity 

bye To remove, via a screen, tramp material (e.g. woodchips) 
which may report in the primary sample stream and which 
could potentially block or otherwise disrupt the 
secondary sample flow. ES 

c) To remove, at least in-so-far as was possible, the air 
which was entrained in the primary sample. This. was 
required to eaeore accuracy in the.pulp density 
measurement. 

Since the secondary sampler design, although commonly 

employed, is not often given in the published literature, a 

sketch is included in Figure 6.2. 

The secondary sample from such a device is expected to 


have a flow rate of nominally 1.6 to 2.3 m?/h (7-10 USGPM). 
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Figure 6.2 A Sketch of the Secondary Sampling Device. 
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With the configuration ef the’ 2.54 cm it in) ‘and 1.3.em (0.5 
in) schedule 40 PVC tubing used in the fabrication of the 
OSCAA sample line, it was estimated from the macroscopic 

' isothermal energy balance that the flow rate would be 
approximately 1.7 m*/h (7.5 USGPM). Subsequent testing of 
the system in operation proved the typical value to be 1.9 
m?/h (8.4 USGPM), which is good agreement considering the 
method of analysis. 

Since the secondary sample flow in the OSCAA sample 
line is essentially vertical a homogeneous flow was 
expected. Using the correlations of Lowenstein(64) it was 
shown that for all three flotation product streams to be 
pesred, the critical velocity“) was ‘exceéded ‘byat! tleasti‘a 
Bactor- "of 2. 

Finally, size degradation studies were performed by 
pumping a slurry through the system for an extended period 
of time. The principal change appeared in the coarsest size 
fraction where, in a four hour period, the mass frequency 
retained on a 0.84 mm (20 mesh) sieve decayed from an 
initial value of 6.57% to a final value of 4.9%. This 
magnitude of change is expected to have a negligible effect 
on the backscattered and fluorescent radiation intensity 
considering the relative mass frequency. In addition, 
comparisons of sample attributes for various slurry mixtures 


were made between the reject and secondary sample lines. 


2°The critical velocity is the minimum velocity ina slurry 
pipeline which will ensure homogeneous distribution of the 
solids in the axial plane. 
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Good agreement was obtained with a typical result shown in 


Table 6.1. 


6.2.2 The Pulp Density Measurement Sub-System 

The equipment used to make the pulp density measurement 
is illustrated in Figure 6.3 and Plate 6.1. 

The radioistope employed in the pulp density gauge was 
Gd-153 which has a principal photon energy of ~100 kev. This 
relatively low energy photon source allowed a gauge design 
with practical length and minimized shielding requirements. 
The "Z-type" geometry was adopted to provide the path length 
required for optimal performance while maintaining 
homogeneous flow. In the design of the pulp density gauge it 
was assumed that the principal source of error would arise 
from the statistical fluctuation of the detected radiation. 
In this case, and as was shown in Chapter 5, pux=2 provides 
the minimum relative error in measurement if Lambert's law 
holds. Depending upon the nature of the slurry, the density 
may vary in the ‘range of-1.0 9s ps1. 1, which, when 
substitution for wu is made, yields an optimal path length of 
11.4 cm < x(optimal) < 12.5 cm. Unfortunately, the physical 
constraints of fabricating the gauge resulted in a mean path 
length of 14.5 cm., close to, but not within the theoretical 
optimum range. (Note that because of beam divergence and 
minimal detector collimation some build-up would be expected 


and hence a Slightly longer path length is desirable.) 
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Table 6.1 Comparison of Sample and Reject Line Slurry 


Attributes 


{ATTRIBUTE REJECT STREAM SAMPLE STREAM 
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% Ash 
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Figure 6.3 A Schematic of the OSCAA Pulp Density Sub-System. 
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Plate 6.1 The OSCAA Pulp Density Gauge 
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The scintillation counter and the necessary support 
electronics are shown on Figure 6.3. A Single >cChannel 
analyser (SCA) was used to examine the energy region of 
interest (ROI) in the amplified detector output (i.e. 50 
kevsSE<100 kev). The SCA provides two parallel output 
channels. One of these signals was connected to the pulse 
counter circuitry in the remote I/O interface. The other was 
connected to a simple scaler for cross-checking purposes. 
The remote I/O was under the real time control of the 
computer program supervising data acquisition. (A more 
complete decription of this equipment is left to the next 
section.) 

The functional specifications of the: pulp density 


Sub-system are given in Table 6.2. 


6.2.3 The Ash Analysis Measurement Sub-System 

The equipment used for ash analysis, or, more properly, 
to obtain the backscattered and iron fluorescence radiation 
measurements, is illustrated in Figure 6.4 and Plate 6.2. 

The radioisotope employed in this system was Cd-109. 
This choice was more fully evaluated in chapter 5. The 
detector was a Xenon filled side window proportional 
counter. The source detector geometry illustrated in Figure 
6.4 is not a conventional configuration (ref. (52) pg. 369 
for example). The geometry used in the OSCAA ash analysis 
gauge was required for a number of reasons which include: 


a) size of source and detector relative to flow cell, 
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Table 6.2 Functional Specifications of the OSCAA Pulp 


Density Gauge Sub-System 


CATEGORY ITEM DESCRIPTION 
Pipeline Piping (a) 2.54 cm > Sch. 40 PVC 
(b) Total Vertical Pass Line Distance - 
Source to Detector ~ 36.8 cm 
(c) Vertical Pass Line Distance Through 
Slurry ~ 14.6 cm 
(d) Inflow/Outflow elbows at 45° 
Source Radio- (a) Gadolinium-153 
Camera isotope (b) Activity when Purchased 130 mCi 
(e@) Half life 241.5 days 
(d) Principal radiations at 97 Kev 
(e) Disk Source 15.2 mm¢d x 8.3 m 
Collimator (a) Length 10.3 cm 
(b) Diameter 0.95 cm 
Casing (a) Overall Outside Dimensions 
16.35cem x) 1261 cm’x 10.2 cm 
(b) Outer case Constructed with 1.27 cm 
Steel Plate 
(c) Case core is aluminum 
Detection Collimator (a) 2.15 cm $ in 1.27 cm Steel Plate 
Detector (a) Model 802-3 Canberra’ Scintillation 
Detector 
(b) Crystal 5.08 cm > by 5.08 cm 
Electrical Preamp- (a) Model 805 Canbercen Scintillation 
Hardware lifier(s) Detector Preamplifier 
(b) Model 802.5 Canberra Tube Base 
Single (a) Model 2015 Canberra’ Amplifier/ 
Channel . Timing Single Channel Analyser 
Analyser 
High (a) Model 3102 Canberra’ High Voltage 
Voltage Power Supply 
Power 
Supply 


Scaler/Timer (a) 


+ 
Canberra Industries Inc. 
45 Gracey Avenue 
Meriden, Connecticut 


06450 


e 


Model 1776 Canberra Dual Counter 
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Figure 6.4 A Schematic of the OSCAA Ash Analysis Sub-System. 


ee 


e*rirs 
Gow 


* 
Ss 
i 


hoe ~ 
iy, e » 
i rv “ 
be * feet 45 
‘ Asi 
a oh f 
ro ei ORR 
= ry al} -. 
aa p 
2.5 te, 


=—vT' a of 
j Ag 
Awe uw fy: Pee ye 
a 
Li 
Yan | 
‘ ) 
* 
‘ee 
we i 
oa” 
Ai 
ert 


& van 


Ce 
‘ Vy ty 
, : 


oe 


: Ce ae 
Bis - 


gue yews: eiaylank ee ania fo vis 


er . 


oe 


160 


The Flow Cell The Assembly For On-Line 


Measurement 


The Ash Analysis Head 


(Source & Detector) 


Plate 6.2 The Ash Analysis Gauge 
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b) the requirement for adequate shielding between the 
source and the detector, 

c) the requirement to minimize the distance between the 
detector and the window to prevent excessive air 
attenuation of the iron Ka X-rays. 

Since the geometry 1S unique a cross-sectional sketch 
is included in Figure 6.5. 

The source and detector housing was equipped with a 
rotating four position shutter mechanism. One position 
provided for maximum attenuation of the primary photons 
required for safe storage. The other three positions 
provided for primary photon attenuation by factors of 0.0, 
0.5 and 0.75. The latter two positions used aluminum filters 
for attenuation purposes. Experience had shown that the 0.5 
position was the best from an operating point of view since 
it prevented detector choking while maximizing the secondary 
photon fluxes. 

The flow cell itself was a 3.2 cm (1.25 in) diameter by 
2.5 cm (1 in.) cylindrical cavity with the sample flow 
perpendicular) to the Jongitudinaliaxis (refi, fig..i6.5). It 
is common practice in XRF applications to provide for leak 
detection and thus prevent sample spillage if the thin 
window ruptures. For the OSCAA sub-system this was 
accomplished by using a dual Mylar™ film covering. Two wires 
were placed between the films such that when the inner 
window ruptured an electrical circuit was closed which 


generated a process interrupt to the remote I/O. 
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Figure 6.5 A Sketch of the OSCAA Ash Analysis Gauge 
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Subsequently, the secondary sample flow was stopped by means 
of a solenoid valve on the secondary sampler discharge. 
While this circuitry was included in the design it was never 
used since the operating life of a single window was more 
than sufficient to complete a given laboratory or plant run. 
Finally, by proper positioning of the flow cell in the 
sample line, the window was kept under slight negative 
pressure thus minimizing warping which could occur with 
changes in the static head. 

The materials of construction of the flow cell and 
source and detector assembly were aluminum and steel 
respectively. These were chosen because of their mechanical 
Propertiestand, in, thescasesotothe datters,, forpits: Shielding 
ability. Interference effects which would result from the 
backscattering of primary photons from the aluminum flow 
cell housing were minimized through collimation. Similarly, 
interference effects which could arise by primary photon 
excitation of the iron in the steel assembly were minimized, 
but not eliminated, by wrapping the detector in a 3 mm (1/8 
in.) lead sheet. 

The electronic components required to support the 
detector and analyse the output are shown on Figure 6.4. A 
multichannel analyser (MCA), under control of the data 
acquisition program was used to acquire the secondary photon 
energy spectrum. When the data acquisition period expired 
the energy spectrum was transferred to computer memory and 


thence to cassette tape by means of a high speed serial 
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interface and a fast buffer. Plate 6.3 shows the primary 
elements of the electronic support and data acquisition 
Systems for both of the measurement sub-systems. 

The functional specifications of the ash analysis 
sub-system and electronic support equipment are given in 


Table 6537. 


6.3 EQUIPMENT SETUP AND OPERATION 

As mentioned in sub-section 5.1.1, the theory of 
operation of the nucleonic equipment employed in the OSCAA 
prototype hardware configuration is adequately described in 
most textbooks dealing with X-ray analysis. The equipment 
operating settings, i.e. detector voltages, amplifier gains, 
etc..., were determined using methods described by Bertin 
(ref. (52): pp. 330-338). No explicit deadtime corrections 
were Hees, these were effectively embodied in the best-fit 


parameters determined on the calibration models. 


6.4 EXPERIMENTAL DESIGN 

The raw material used to generate the samples required 
for the laboratory experimentation was obtained from 
McIntyre Mines Ltd. The material consisted of two large 
Samples of flotation concentrate and one of flotation 
tailings. All three of the samples were individually mixed 
and sub-sampled for bagging and storage. A representative 
Sample of each was analysed for ash content, iron 


concentration of the ash and solids specific gravity. 


MA i } hy pe deel 1 
ey AAP EE NM a NS 
Hie ee 
bso PY fa) 


vy at? awlede 6.3 ened sees ge on 


" Ae as eis wisish> Sire SAOgEIS! 


unovawanduh 9. Hes 7) a eee ans 
bey lpia dee) mee vo -anotysohiigagae radish 


eon shmey eRe dxeogaue oénonapate & 


ics ayy wi a 
i 

fT é. 

at 5 


fi) Wewl 


t 


~ ies 5. Wij ali ae hea 
HOITANRED joie Gorne. 

to wiowde ofa) thie opbavaa~dise Loa keno aa 
CASO oth AE beatae. vein Faye rhc f sistt pegs 
‘i. Setipedes “ys etaupeos 8.2) aos¢9 iupitnes. esentinnd 9 
comykepe at? .ciarlana yep doew onit ied ood te: 23 
yehts iq , cap szioy yotos?eS sek epndsss é 
si yiet ed bedi rres® ehodvene paige henimra anh 9° 7 us | 
eres i oenuen enltybeahlaks i figaten* “ah eee sa (sa 

j-49ged 902.0 than role 0 Dey te ape pxaw ad di: 
= 2hat eae ae brie 


| WaTeae, 
boxtopet eslgme od3 stmenee oo basw let ie san’ aniseed 
mous bea sido: naw noltacapeiseges yiosesadss) 
ere wet Fo sedetanes: ts} 193 om oat .B2d4 
feltescti3 io ened a} a2dNEDAOD ne iansol3”. 
boxte QhanuOivibal oie shlowoe add to setds LEAs 5 
evideationeiye: + ,ogelmie Sap pei pred OF ‘eataan . | 
es ‘He den soon due So? Bery/ sane canis penne 


: sivary 11 Lomgailaitios: bak das a ahi tone 
7 | cam | oe sh aes 


165 


Interfaces Real-Time 


Preamplifiers HPIB/ZRSZ32¢C Clock 


NIM Bin Multichannel Remote 
(Power supplies, 
SCA, scalar/ 
timer) 


Analyser impuc /OutpuUe 


Printer Plotter 


Cassette Control 


Lape Computer 


Plate 6.3 The OSCAA Electronic Support and Data Acquisition 


Systems 
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Table 6.3 Functional Specifications of the Ash Analysis 


Sub-System and the Electrical Hardware 


CATEGORY ITEM DESCRIPTION 
Pipeline Piping (a) 1.27 cm > Sch. 40 PVC 
Flowcell (a) 3.18 cmd x 2.54 cm 
(b) Window Material 10 yp Mylar 
Analyser Head Radio- (a) Cadmium 109 
Source isotope (b) Activity When Purchased 35 mCi 
(c) Half Life 453 days 
(d) Principal radiation at 22-26 
kev(Ag k x-rays) 
(e) Disk Source 11.4 mmd x 5.4 mm 
Collimator (a) Length 2 cm ‘ 
(b) Diameter 0.64 cm 
Housing Casing (a) Fabricated from Steel 
(b) Approximate Overall Dimentions 
15 cm x 8 cm 
Detection Detector (a) Model 45415 LND* Sidewindow 
Prop. Counter 
(b) Xe Filled at 1 Atmosphere 
(c) Berylium Window of 
(i) thickness ~ 0.13 m 
(ii) size 2.54 cm x 5.1 cm 
Electrical Preamp- (a) Model 1406 Canberra Proportional 
Hardware lifier Counter Preamp. 


(Energy Disper- Multichannel (a) 
sive Analysis) Analyser 


High Voltage (a) 
Power Supply 


MCA Control (a) 


Model 8100 Canberra MCA with 
Nécessary Optional Features 


Model 3003 Canberra HVPS 


Model 8100-14 Serial Computer 


Interface Interface - 
Electrical Computer (a) Model HP9825 Hewlitt Packard 
Hardware : Desk Top Calculator With Necessary 
Features < 
(Data Acquis- Remote (a) Model HP6940B Multiprogrammer 
itions, Stor- Input/Output with Necessary Optional Features 
age, Presenta- (b) Model HP59500A H.P.I.B. Interface 
tion and To The Multiprogrammer 
Processing) 
Interfaces (a) Model HP98036A RS232C 
Serial Interface 
(b) Model HP98034A H.P.I.B. 
Parallel Interface 
Clock (a) Model HP59039A ASCII Clock 
Plotter (a) Model HP9872A Graphicoplotter 
Printer (a) Model HP9871A 30 Baud Impact Printer 
Mass Storage (a) Model HP9162-0061 Tape Cartridges 
Medium 
*LND Inc. 


Oceanside, New York 
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For the purposes of organizing a particular performance 
evaluation on a given product stream, orthogonal design 
principles were adopted. For the first three runs (feed, 
concentrate and tailings) a 2-ftactor (%solids, t%asnhy, 
5-level design was used yielding 25 calibration points. The 
fourth run was a’3-factor (%solids, Yash, %iron), 3-level 
design yielding 27 calibration points. The advantage of 
uSing an orthogonal design is that it provides for the best 
estimates of calibration model parameters by minimizing 
covariance among the factors. From an analysis of the 
preliminary testwork data described in chapter 3 and through 
discussion with the operators, the design factor levels 
given in Table 6.4 were identified. 

To make maximum use of the the raw coal material which 
was available it was necessary, wherever possible, to move 
from point to point on the design by simple "make-up" 
additions. (By the nature of the experiments there was no 
need to randomize the order.) Toward this end, an 
interactive computer program was used to solve the system 
constraint equations’' and determine the water, clean coal 
(concentrate), reject coal (tailings) and magnetite 
requirements of the slurry. Based on component inventories 
after sampling for the preceding design point, the 
incremental additions of the various components were 


calculated. If negative increments were encountered the 


2'The constraint equations included %solids, %ash, %iron and 
the total system volume. The latter was targetted at 38 
litres after the standardization run with water. 
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Table 6.4 The Orthogonal Design of the Laboratory Runs 
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component exhibiting the largest relative negative change 
was fixed at its current level, the volume constraint 
removed and the additions recalculated. If this calculation 
provided a slurry volume in excess of 46 litres, a new 
Slurry was mixed for this design point. Calculations for 
"bleeding" slurry from the system and adding water to reduce 
only the solids concentration were done by hand since this 
happened only infrequently. 

Having calculated the necessary water and solids 
additions for each design point in a particular run, the 
solids components were obtained from the bulk raw samples, 
mixed in proportion and bagged with appropriate labelling 
fOreuse ins the calibration run. For illustrative: purposes 
the experimental design for the flotation feed calibration 
is presented in Figure 6.6. The points at which a new slurry 
was mixed, where the system was bled and the moves involving 


incremental additions are shown. 


6.5 EXPERIMENTAL PROCEDURE 

The steps involved in a calibration run may be divided 
into two parts, operational and sample handling. The former 
involves all activities up to and including the stream 
Sampling and experimental data storage and the latter 


involves all subsequent sample handling and analyses. 
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6.5.1 Operational Methods 

Each run was divided into a series of design points. 
For example, in Figure 6.6 there are four series of points 
(labelled as @). Each series began with the preparation of a 
new slurry and testing within a series involved incremental 
additions and bleeding. 

Prior to the start of a new series the system was 
allowed to "run-in" on 36 litres of water. This provided for 
instrumentation warm-up, room temperature equilibration 
etc... Steady state was established by monitoring the the 
secondary photon counts. When this conditon was achieved, 
the appropriate readings were taken on the water standard 
and stored for future reference. Although this particular 
experiment does not appear in the orthogonal design it was 
conducted for each series. 

Once the water standard had been measured the necessary 
incremental additions were made to the sump, slowly, to 
facilitate good mixing. The system was then allowed at least 
10 minutes of circulation to ensure slurry homogeneity. Upon 
completion of this stage, the data acquisition period was 
initiated by means of the computer program. When the program 
suspended acquisition and began data transfer and storage 
operations, a 2 litre sample of the slurry was taken from 
the secondary sample line. When this was complete the 
incremental additions for the succeeding test were made and 
another. 10 minute mixing period allowed prior to data 


acquisition. 
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The decision to take 2 litres of sample was made on the 
basis of particle size distribution requirements. Although 
particle size cannot be explicitly accounted for in the 
‘' calibration models présented in chapter 5, a knowledge of 
the relative changes was deemed necessary for possible 
correlation analysis with calibration results (i.e. analysis 
of residuals)??. Using the flotation feed for design 
purposes, specifying 95% confidence limits and using 
Barbery's(65) equations for sample size determination, at 
least 100 g of solids were required. For a slurry of 5% 
solids, 2 litres were required to provide this weight. It 
was decided to fix the sample size at 2 litres for all 
Peses;esince, withvthesexception of esome ofethe design 
points in the tailings run, the 100 g requirement would 
always be met. Furthermore, 2 litres of sample is relatively 
easy to handle and is small enough to ensure the maximum 
usage of the raw coal resources. 

Between each series the system was thoroughly flushed 
with water. The next series was begun with the addition of 
36 litres of water and the warm-up period. 

After each calibration run equipment maintenance was 
performed and the Mylar™ window on the flow cell replaced. 
The hard copy output and the data tape were removed, 


labelled and stored for the subsequent analysis. 


72From the theory, modest changes in the particle size 
distribution were not expected to significantly affect 
results. 
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6.5.2 Sample Handling 

The 2 litre sample from each design point was obtained 
in two 1 litre beakers. The wet weight of the slurry was 
recorded and it was then transferred to a drying pan. Rather 
than use filtration to expedite drying, evaporation under 
infra-red lamps was employed to avoid potential sample 
losses. The dry solids weight was recorded and the solids 
transferred to a 0.6 mm (28 mesh) sieve. Using hand 
Screening and gentle pressure the fine coal agglomerates 
were broken down. The solids were then thoroughly mixed and 
a representative 100 g sample was cut for size, ash and iron 
analysis. In the Special cases of runs 1 and 2, 20 g samples 
werienalso; cuts fora specific, gravaty,.determinations. Any 
reject material was then bagged and placed into storage. 

The 100 g sample was screen analysed (dry) on the 
following sieve nest; 0.84 mm, 0.6 mm, 0.42 mm, 0.3 mm, 
O.21mm, 0.15mm (20/28/35/48/65/100 mesh). The results were 
recorded and the screen fractions recombined. Minimal weight 
losses were observed in this step. (In the case of run 3 
where some of the slurry samples contained only =20 g, 
screen analyses were not performed.) 

The recombined size fractions were then reduced to -0.1 
mm (-150 mesh) in a coffee grinder. (Independent studies had 
shown no iron contamination from the stainless steel 
blades.) From this very fine sample, duplicate 3 g cuts were 
obtained for ash analysis. The ash analysis was performed by 


heating the coal in a muffle furnace at a temperature of 750 
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C for a period of approximately 6 hours. Inspection of the 
uncombusted residue showed no trace of organic (black) 
substance. The reject was bagged, labelled and placed in 
storage. | 

The ash residue from the ash analysis step was weighed 
and the ash content calculated. Except in the very 
infrequent case of error, where one sample quite obviously 
gave incorrect results, the mean of these two assays was 
reported as the solids ash content. 

The duplicate ash residues were combined and submitted 
for iron analysis. This analytical work was done ina 
commercial laboratory. Despite the fact that both a 
method(65) and the equipment for performing these analyses 
were available, consistent results could not be obtained. 
(This is thought to be due in large part to the state of 
repair of the atomic absorption analyser.) Thus, a local 
laboratory with a plasma spectrometer was chosen to do the 
work. To reduce analytical costs only selected samples from 
calibration runs 1, 2 and 3 were submitted. It was argued 
and later shown (ref. Chapter 7) that, since iron was not 
independently (externally) controlled in these calibrations, 
it would be strongly correlated with the ash content. 
Therefore, a correlation constructed on the basis of a few 
samples could be used to predict the iron concentrations for 
those samples not analysed. All of the ash samples for run 4 


were submitted for iron analysis. 
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The results of the attribute analyses were tabulated 


and stored for the subsequent calibration studies. 
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7. THE OSCAA SYSTEM LABORATORY CALIBRATION RESULTS 


7.1 INTRODUCTION 

This chapter describes the methods employed in reducing 
the data derived from the laboratory experiments and the 
results of the subsequent calibration studies. Details on 
the MCA spectrum analysis, semi-empirical model parameter 
estimation and data preparation prior to calibration are 
provided. The criteria for assessing the performance of the 
OSCAA prototype is defined and finally the results of the 


calibration studies are presented. 


7.2 MCA SPECTRUM ANALYSIS 

The MCA spectra acquired for the different experiments 
were analysed to provide measures of the backscattered and 
iron fluorescence counts required for the calibration 

Studies. This analysis was complicated by two factors: 

a) the position of the spectrum with respect to energy was 
a function of the count rate. ( This is consistent with 
the count rate effect.) 

b) due to the relatively poor resolution of proportional 
counters and the choice of the detector gas (Xe), there 
was a background discontinuity which underlay the iron 
peak. 

These effects were only important in the case of the iron 

peak since peak shift effects were negligible and there were 


no background discontinuities in the case of the backscatter 
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peak. A typical MCA spectrum for an experiment is presented 
in Figure 7.1. It is convenient to consider the methods of 


peak resolution separately. 


7.2.1 Backscatter Peak Resolution 

Studies on the nature of the background both on the low 
energy and on the high energy side of the backscatter peak 
indicated that it could be modelled via the simple linear 


expression: 
B=b)+b,C 
where: B =background counts in channel number C. 
C =channel number 
bo, b, =zcoefficients of the linear model. 
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Since Xe has no absorption edges in the energy region of 
interest (ROI), there was no reason to suspect a 
discontinuity. To estimate the linear model coefficients in 
equation (7.1) a "quick" computation was used rather than 
the more mathematically (statistically) correct linear 
regression analysis. The quick method invalved the following 
steps: 

a) Define a region on the spectrum to the low energy side 
of the peak where the contribution of the peak is 
negligible. Calculate the mean background count for this 
region and the centroid channel. 

b) Repeat the calculations on the high energy side of the 
peak. 

c) Using the two points defined in (a) and (b) above, 


estimate the coefficients of the background model. 
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Figure 7.1 A Typical Experimental MCA Spectrum. 
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The magnitude of the counts in the channels which 
define the backscatter peak preclude the need of any 
mathematical manipulation of the raw data such as smoothing. 
- Once the low energy ane high energy peak boundary channels 
were defined, peak integration over the ROI is performed 
Simply by subtracting the predicted background from the 
observed counts for each channel within the boundaries, and 
summing the results. This procedure is illustrated 


schematically in Figure 7.2. 


7.2.2 Iron Peak Resolution 

The problems associated with analysing the iron peak 
occur because of spectrum shifting and an underlying 
background discontinuity. 

Peak shifts result when the gross count rate at the 
detector changes and the eesciacicn of the individual pulses 
is not complete (ref. (67) pg. 75). Although these shifts 
generally have a range of only a few channels in the MCA 
Spectrum for a particular stream, they must be accounted for 
to ensure maximum accuracy in deriving the iron fluorescence 
count data. This does not represent a particularly difficult 
problem since all of the peak analysis methods studied were 
adaptive (i.e. capable of responding to peak movement). 

The discontinuity in the background which underlies the 
iron peak represents the more serious problem in resolution. 
It was strongly suspected that the discontinuity was 


correlated with the Xe detector gas absorption edges (ref. 
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Figure 7.2 A Schematic Representation of the Method Employed 


in Resolving the Backscatter Peak. 
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(52) pg. 305). The mass absorption coefficient for Xe over 
the energy range of interest is shown in Figure 7.3. The 
principal iron fluorescence peak, Fe Ka at 6.4 kev, is 
observed to lie just to the high energy side of the Xe L 
shell absorption edges. Given the resolution of the 
proportional counter it was expected that there would be 
some overlap resulting in a discontinuity which would 
underlie the iron K fluorescent radiation peaks. 

A special experiment was designed where it was possible 
to directly irradiate the detector with the Cd-109 source, 
both with and without iron fluorescence. A typical result 
for the experiment is shown in Figure 7.4. It is readily 
apparent from this figure that overlap does indeed exist. 
Calculations on the spectral data indicated that the 
beginning of the discontinuity on the high energy side 
corresponds to the Xe L1 absorption edge (5.5 kev). 

Clearly, the method used to extract the iron 
fluorescence counts should be somewhat more sophisticated 
than that which was employed to resolve the backscatter 
peak. Since it is desirable, from a statistical accuracy 
viewpoint, to use the area under the peak as a measure of 
the fluorescence counts a peak modelling approach was 
intially adopted?®. 

The modelling approach, while more mathematically 
elegant, proved to be less precise than a simple resolution 


?32This is the preferred method of measurement in energy 
dispersive analysis as opposed to the single energy 
(wavelength) measurement common in wavelength dispersive XRF 
analysis. 
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Figure 7.3 The Mass Absorption Coefficient for Xenon as a 


Function of Photon Energy. 
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Figure 7.4 Identification of the Discontinuity Underlying 


the Iron Fluorescence Peak. 
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technique, and for this reason was not used in the OSCAA 
Studies. However, in view of the potential application of 
peak modelling in the next generation of instrument 
' development, the methodology and the results of the 


resolution studies are described in the text which follows. 


7.2.2.1 Modelling the Iron Fluorescence Peaks 

The iron fluorescence radiation peak model is given 
in equation (7.2). The model assumes a double Gaussian 
function to account for both the Ka (6.4 kev) and KB 


(7.05 kev) contributions. 


rages jor E faecal" snes E [ss ey eebie see) 


where: vy, Ssiron hk Ltluorescent radiation counts. in 

channel 1. 

Y' =iron Ka fluorescent radiation counts for 
"a-peak" modal channel. 

k =numérical constant (=41n(2)). 

1 =channel number. 

Ca =modal channel for the Ka radiation. 

Wa =full width at half maximum (FWHM) for the Ka 
peak. 

CB =modal channel number for the Kf radiation. 

R =ratio of the KS to Ka peak heights. 

WB =FWHM for the Kf peak. 


The unknown parameter set in the equation (Y', Ca, Wa, 
CB, R, WB) can be reduced from the theory. Jenkins and 
Devries((67) pg. 175) show that R is approximately equal 
to 0.15. Using a steel sample in a scanning electron 
microscope equipped with energy dispersive XRF analysis, 
R was observed to be 0.143. Moreover, since the MCA 
Spectrum was linear in energy (channel number) with a 


zero intercept, the ratio of C8/Ca was calculated to be 
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1.102 and therefore; 
CB=1.102Ca aah «3) 
In addition, the theory states that; 
wap??? ee op eA) 


and, therefore; 


Os 
fs 05 
we [2-95] Wa=1.05Wa Saas. 5) 


Equation (7.5) was corroborated by the scanning electron 
microscope work. Substitution of equations (7.3) and 
(7.5) and the independently obtained value of R into 
equation (7.2) produced a model in three unknown 


Parameters, namely Y', Ca, and Wa. 


7.2.2.2 Modelling the Background 

The background model was assumed to take on one of 
two forms. The first form considered the background to 
be constant on the high energy side of the 
discontinuity. This meant that only the data in this 
region of the spectrum could be analysed, however, as 
shown in Figure 7.4, the fluorescence peak is relatively 
well defined in this region. It is possible, in 
principle, to fit the combined background and iron 
fluorescence models over this region. If Cm is taken as 
the channel corresponding the maximum counts in the 
observed spectrum, the constant background model can be 
estimated as follows; 
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where: B =estimate of the background counts on the high 
energy side of the discontinuity. 
n =the averaging increment (taken as 2 in this 
work). 
Y,; =observed counts in channel j. 
Cm =modal channel in the observed spectrum. 

The second form of the background model was derived 
from inspection of plots of the first derivative of the 
"pure" background spectra from the special experiment 
(ref. fig. 7.4). Since the resulting curve resembled a 
Gaussian density distribution it was argued that the 
background could be partially modelled as an integrated 
form of this function. This required that the background 
be considered as a three part function. It would be 
constant both below and above the discontinuity and be 


equivalent to the integrated Gaussian over the 


discontinuity (refietig.. 724: -*pure: background”). 


Mathematically; 
it+n (7.7) 
B — ° ;— eee e 
L sasaon J 5 *3/ant1) et O0y: cm 
j=i-n 
itn 
BU=MINIMU dee Calan te OS 1 es ONG) 
oa Stan on x j (ane) ty ay oo) 


where: BL, BU =constant background on the lower and 
upper energy side of the discontinuity, 
respectively. 

Since the total counts under a Gaussian peak, I, is 

related to the FWHM, Wb, and the counts in the modal 


channel, Y'', as; 


T=Y' 'Wh/2T 
2.35 Syoont <9) 
wheres I =total counts under the background "peak". 
Y'' =counts in the modal channel for the 


background "peak". 
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Wb =FWHM for the background "peak". 
then, if BL and BU are fixed according to equations 


(Zohn sone. (a8 bs 


es a (re, 40) 


and combining equations (7.9) and (7.10) with 


rearrangement yields; 
_2.35(BU-BL) 
ViotoaYe! 


Thus the unknown parameters for an integrated Gaussian 


Wb el aT 
function describing background would be Cb, the modal 
channel, and Y'', and the final model for the background 


is given in equation (7.12). 


i See 
B.=BLt+ ee exp ~« [Ae] ocala) 


j=Hz 


where: B, =background counts in the ith channel. 

4 =channel corresponding to BL. 

Cb =modal channel for the background "peak". 
7.2.2.3 Modelling the Spectrum 

Considering the discussion above, there were two 

methods of modelling the photon energy spectrum in the 
region including the iron peak. The first was to examine 
the spectrum on the high energy side of the 
discontinuity. In this case equations (7.2) and (7.6) 
were used, where, for the latter, B was also considered 
an unknown whose intial estimate was obtained from 
equation (7.6). This was termed the "1/2 Peak Method" 
Since it examined a little more that one-half of the 
peak. The other method examined the entire peak using 


equations (7.2) and (7.12) and was termed the "full peak 
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method". 

For comparative purposes a very simple extraction 
technique was also studied, analogous to wavelength 
dispersive XRF ayaa st In this case the background was 
estimated from equation (7.6) and the aenes in channel 
Cm, corrected for background, B, were taken as the iron 
radiation measurement. This was termed the "simple 
method". 

Parameterization of the models was accomplished 
using an orthogonal multivariable search algorithm (an 
adaptation of the Coggin's algorithm; (68)). Initial 
estimates of the parameters were available from previous 
work. The objective function minimized in the 
Parameterization exercise was the standard x? criterion 


having the form; 


b eee, 
Z ere 
F=2 Y, ewe td al Ss) 


a 
where: Yj, =predicted counts in channel i. 
a, b =channel limits defining the range of fit. 


Once the model parameters are known, the area under the 
iron K fluorescent radiation peak is analytically 
available from expressions of the form given by equation 
(7.9) (summing a and 8 contributions). 

With regard to the spectrum itself, the raw data 
showed considerable statistical fluctuation as was 
expected at the relatively low count levels. As an 
alternative to a direct fit on the raw data, a smoothing 


algorithm was developed as described by Marr(69). The 
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spectrum smoothing procedure is recommended by Bevington 
(ref. (70) pg. 250) to improve the fit by reducing the 


x? 


value, but not without warnings on potential problems 
related to inereated uncertainty on the parameters. 
Briefly, the smoothing algorithm requires that the raw 
count data be transformed by taking logarithms. Then, 
uSing linear regression’?‘*, and assuming a quadratic 
relationship, the transformed counts in a particular 
channel were estimated by fitting the quadratic over 
(2n+1) points, n on either side of the channel of 
interest. For this work n was taken to be 12. In 
accordance with Marr's suggestion, the smoothing 
operation involved two passes with the algorithm, the 


results of the second pass were re-transformed to 


produce the smoothed spectrum. 


7.2.2.4 Results of the Resolution Studies 

The three extraction methods were then tested on 
both raw and smoothed spectra to determine which would 
provide the most precise results (accuracy could not be 
tested). For the purposes of the extraction methodology 
testing a series of 4 spectra on pure water samples were 
used. The water spectrum shows a pronounced iron peak 
due to the steel used in the fabrication of the source 
and detector housing. It is, however, the minimum peak 


height obtainable and thus is most susceptible to 


74The regression is done using convolution integers because 
of the equispaced channel data. 
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Statistical fluctuations. For this reason the water 
Spectra make the best test-bed for methodology 
evaluation. It was argued that the method selected must 
be reliable for all Spectra acquired and therefore must 
be accurate for any subset, including the pure water 
spectra. 

A selection of the results of the evaluation 
Studies is presented in Table 7.1. In Table 7.1, (A) and 
(B) compare the "1/2 peak method" on both the raw and 
smoothed spectral data. For comparative purposes both Y' 
and I are reported. In both cases Y' was a more stable 
estimator of the fluorescent counts than I, as reflected 
in the relative error of estimation. This is the result 
of the additional uncertainty in the FWHM, which is 
required in addition to Y' to compute I. It is clear 
from the differences in the relative error figures that 
the smoothed data gave more stable estimations than did 
the raw data. (B) and (C) compare the "1/2 peak method" 
and the "full peak method" as applied to smoothed data 
demonstrating the superiority of estimation for the 
batter. This: is) thought.to be: the result of the “full 
peak method" including the entire range of the iron 
fluorescence peak thereby improving parameter 
estimation. Finally, (D) demonstrates the results for 
the "Simple method" which was judged to be the best 
method for resolving iron count data. (The trends shown 


here were later substantiated in the calibration 
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Table 7.1 Results of Iron Peak Resolution Studies 


(A) 


(B) 


(C) 


(D) 


"1/2 Peak Method" on Raw Data 


Std. Dev. 
Rel. Err. 


"1/2 Peak Method" on Smoothed Data 


60442.9 
62660.2 
61908.6 
67450.2 


Mean 
Sta. Dev. 
Rel eaoie 


"Full Peak Method" on Smoothed Data 


counts 


igs 
ve ae aaa 


65256.8 
74752.6 
61109.4 
73487.2 
Mean 1378.1 68651.5 
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exercise for run 1 on flotation feed.) 

On the basis of these findings the "simple method" 
was selected for the subsequent spectral analysis work. 
Figure 7.5 illustrates the mechanics of the "simple 
method" for extraction of the iron count data. Briefly, 
the raw data is first smoothed, the modal channel 
identified and the background estimated whence the 
corrected counts in the modal channel, Ymax, is 
calculated. 

It was indeed unfortunate that a reliable estimate 
of peak area could not be developed within the time 
frame of this project. A reliable estimator of this form 
would have provided statistically better estimates of 


the iron counts. 


7.3 PARAMETER ESTIMATION FOR THE SEMI-EMPIRICAL CALIBRATION 
MODELS 

The four component semi-empirical OSCAA calibration 
model given by equations (5.76) through (5.79) requires that 
various mass absorption and scattering absorption nuclear 
coefficients be available for each component in the assumed 
physical model of the slurry. In addition, estimates of the 
component densities are required. In principle these 
parameters could be estimated by including them as unknowns 
in the model calibration exercise (together with the A;,,B;). 
However, it is considerably simpler to calculate the nuclear 


coefficients from the data available in the literature, 
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Figure 7.5 Resolution of the Iron Fluorescence Peak by the 


Simple Method. 
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knowing the approximate elemental composition of the 
components. Furthermore, the densities can be estimated by 
fitting the assumed physical model of the coal solids to 
"actual specific gravity data derived from the analysis of 


experimental samples. 


7.3.1 Calculation of the Nuclear Coefficients 
An interactive computer program was written (Appendix 

1) and the appropriate data base constructed from the 
tabulated values(39) such that upon entry of: 

(i) the ash analysis, 

(ii) the ultimate analysis of the coal, 
the appropriate nuclear coefficients were estimated from the 
elemental composition. Included in the output from this 
program were the mass absorption coefficients at 6.4 kev 
(principal Fe fluorescence photon energy) and 22 kev 
(average principal Cd-109 photon energy), as well as the 
combined scattering absorption Seerriciencss (The later were 
estimated by calculating o, and oym, correcting them with 
lm;r Unm (Sub=section!’ 525.1.3)>)°and adding the results 
together). These various coefficients are computed for: 

(i)the organic phase, 

(ii)the overall mineral matter phase (including the 

iron component), 

(iii) the iron free mineral matter phase. 
Since the composition of the other two possible components, 


namely water and hematite, are fixed, these calculations 
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were done only once. A sample of the relevant portion of the 
computer program output is presented in Figure 7.6. 

Several ultimate and ash analyses were available on 
McIntyre Mines Ltd. coal (see, for example, (71), pg. 55). 
These were done at different times and by different groups 
and included results obtained on the raw coal samples used 
in this work. Interestingly, the ultimate analysis changed 
very little and the most significant differences between the 
ash analyses were the concentrations of iron, This suggests 
that, in this particular case at least, the nuclear 
coefficients will be reasonable for the deposit as a whole 
thus the four component model assumed for calibration 


purposes was quite reasonable. 


7.3.2 Estimating Solid Component Densities 

It was explicitly stated in the discussion in 
sub-section 5.7.1 that the coal solids were assumed to 
consist of least two phases, minimally an organic phase and 
a mineral matter phase. The coal solids density, in terms of 


these two phases, may be written as; 


iste eee 
p= A 
S| Glam) tM Per fetes) 
Po ae 
where: p, =coal solids density. 


P2, P3' =coal solids organic and mineral matter 
phase densities. 
m =mass fraction of mineral matter in the coal 
solids. 
To estimate pz and p3;', the specific gravity data derived 
from run 1 was used as it provided the widest range of 


values. A non-linear regression program(72) was employed to 
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ORGANIC PHASE CALCULATIONS.... 


CARBON HYDROGEN NITROGEN SULFUR OXYGEN ASH 
WT. % (RAW) 83.70 4.44 1.26 0.49 P) NE) Se 
WT, % (ADU) SORon 4.80 1) 2 Shs 0.53 22,30 0.0 
TOTAL AC(6.4)= 9.5469 TOTAL AC(22)= 0.4069 INCOH. SCATTER (22)= 0.1704 


COMpeSGCANTER (22) == O07 O57. 4a TONAL SCALLER  (22)i= ~O7O0SO6 


INORGANIC PHASE CALCULATIONS.... 


OXIDE OF WEIGHT % WEIGHT % EXCLUDING 


RAW ADJUSTED HEMATITE 
SILICON 52 16 53.98 58.63 
ALUMINUM 23.65 24.48 26.58 
CALCIUM 3.27 3.38 3.68 
MAGNES. 0.74 ORT 0.83 
TITANIUM 1-89 1.96 Pa 22 
SODIUM 0.84 0.87 0.94 
POTASS. be CU 1.42 1.54 
PHOSPH. Oly 1.00 IO? 
SULFUR 4.08 4.22 4.59 
IRON T5695 Te O2 0.0 


DATA FOR ASH FORMING MINERALS: 


TOTAL AC(6.4)= 79.4664 TOTAL AC(22)= 2.2408 INCOH. SCATTER (22)= 0.1503 
COH: SCAT TERT (22) =O, 1507 TOTAL, SCATTER (22)= 0.0499 


OXYGEN FRACTION IN THE IRON FREE INORGANICS IS 49.87 % 


DATA FOR MINERAL MATTER: 


TOTAL AC(6.4)= 74.0140 TOTAL AC(22)= 2.0941 INCOH. SCATTER (22)= 0.1524 
COH. SCATTER £022) 2 40, d4430 FOTAL SCATTER (22)) = .0.0502 


DATA FOR HEMATITE: 


TOTAL AC(6.4)= 55.9200 TOTAL AC(22)= 13.8680 INCOH. SCATTER (22)= 0.1301 
COH:. SCATTER (22)= 0.3500 TOTAL SCATT€R (22)= 0.0549 


DATA FOR COMPOSITE MINERAL MATTER; 


wee ee a a a a a a a a a we ee ewe we www ww ee ee 


TOTAL AC(6.4)= 72.5814 TOTAL AC(22)= 3.0263 INCOH. SCATTER (22)= 0.1507 
COH. SCATTER (22) 0.1606 TOTAL SCATTER (22)= 0.0506 


Figure 7.6 Sample Program Output for Calculating Nuclear 


Coefficients in the Semi-Empirical Calibration Models. 
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calculate the estimates. The raw data and the regression 
results are presented in Table 7.2. 

For the purposes of the four component model, where it 
"was necessary to provide a p3; value for the iron free 


mineral matter, equation (7.14) was used to transform p;3'; 


Ci) 
£|\hiea- 2 
P3 Rar eam Pee) a sy, 
3 4 
where: p3 =density of the iron free mineral matter. 


py =Gensity of hematite. 
f=mass fraction of iron bearing mineral in the 
mineral matter phase. 


7.4 DATA PREPARATION 

This section deals with the methods employed to 
manipulate the raw data prior to performing the calibration 
studies. It is divided into two sub-sections, the first 
dealing with slurry attribute analyses and the second with 


the radiation count data. 


7.4.1 Preparation of the Slurry Attribute Data 

As described in sub-section 6.4.2, three slurry 
attributes were measured; 

a) % solids in the slurry, 

b)” % ash in’ the solids, 

GJ ee Lronin the ash. 
The analyses for (a) and (b) were done in-house but for (c), 
the analyses were performed by a commercial laboratory. In 
ordér to minimize’-costs for the latter, insall but run 4, 


where the iron concentration was independently controlled, 
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Table 7.2 Estimates of the Solid Phase Densities in the 
McIntyre Coal 
RAW DATA FOR FLOTATION FEED CALIBATION RUN 
; . 4 ++ 
Experiment Mass. Fraction Sodas 


Number of Mineral Specific 
Matter (m) 
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only a few of the selected ash residues were submitted for 
iron analysis. It was expected that a systematic 
relationship between iron content and ash would be observed 
for the simple raw eon blending used in runs 1 through 3. 

A quadratic function can be used, over a limited range 
of ash values, to characterize the iron versus ash 
relationship. A typical result of this approach is 
irustrated=in- Tabte-7~>s Lor“run tirana a’ portion of run 2 
where the same raw coal samples were used for blending. The 
quadratic approximation is a relatively accurate predictor, 
and through substitution of the ash values for particular 
samples the corresponding iron content can be estimated. 
These estimated values were used in the calibration studies. 
In the cases where the regression was insignificant the 
model predicted relatively constant iron contents, as would 
be expected, and these values were used. 

Calculation of the mass fractions of the various slurry 
components was left to the computer programs which perform 
the calibration model parameterization. Tables 7.4 through 
7.7 present the attribute analyses for the four laboratory 


runs. 


7.4.2 Preparation of Radiation Count Data 

The backscattered radiation data derived in the manner 
described in section 7.2 were not manipulated prior to the 
calibration, however, both the pulp density and iron 


fluorescence counts were modified. The method of handling 
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Table 7.3 Estimating Iron Content From Ash 


RAW DATA FROM FLOTATION FEED CALIBRATION (RUN #1) 
AND FLOTATION CONCENTRATE CALIBRATION (RUN #2) 
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Table 7.4 Slurry Attribute Analyses for Laboratory Run 1 


EXPERIMENT ho SOLIDS IN *% ASH IN 4 IRON IN 
. | NUMBER SLURRY SUELDS ASH 
S26 0725 
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Table 7.5 Slurry Attribute Analyses for Laboratory Run 2 


EXPERIMENT % SOLIDS IN % ASH IN *% IRON IN 
NUMBER SLURRY SOLIDS ASH 


{ 22.08 Tan 5.69 
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Table 7.6 Slurry Attribute Analyses for Laboratory Run 3 


EXPERIMENT % SOLIDS IN *% ASH IN % IRON IN 
NUMBER SLURRY. SORBTOS ASH 
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Table 7.7 Slurry Attribute Analyses for Laboratory Run 4 


EXPERIMENT % SOLIDS IN *% ASH IN % IRON IN 
NUMBER SLURRY, SOLEDS ASH 
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the pulp density counts was a two step process. 

a) Subtract Background: In this case the background was 
estimated by closing the source camera shutter, which 
effectively totally shields the Gd-153 source from the 
detector. A series of 300 second counts were the 
recorded and the mean value (typically 40,000 counts) 
was scaled up to the corresponding data acquisition 
period and subtracted from the raw count data. 

b) Reduction by Standardization: The corrected count values 
obtained in (a) were divided by the corresponding values 
for the pure water experiments run prior to each series. 

The corrected and reduced count data from this process were 

used in the calibration studies. 

For the iron fluorescence count data, each experimental 
value was corrected by subtracting the corresponding value 
for the pure water experiments. 

In summary, all of the radiation count data was 
corrected for background at some point er pncehec in the 
spectral analysis or data preparation stages. It is common 
to reduce density count information as well as other 
Spectral data, utlizing standards, to account for equipment 
drift etc... In this case only the density count data was 
reduced since this has the effect of partially compensating 
for varying entrained air content between series. The iron 
count data was corrected from the pure water standards to 
remove the contribution due to the fluorescence from the 


Steel assembly housing the source and detector in the ash 
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analysis sub-system. The radiation data used in the 


calibration studies is presented in Tables 7.8 through 7.11. 


7.5 THE CRITERIA FOR ASSESSING CALIBRATION RESULTS 

There are a variety of means of assessing the 
performance of an on-stream analysis device intended for use 
in process control. It is the opinion of the author that the 
best criteria for evaluating OSCAA is a modified version of 
that defined by Howarth et. al.(73). If the unit meets the 
following criteria for both % solids and % ash prediction it 


is acceptable for process control applications. 
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eee N nea 19) 


where: N =number of experiments. 

S(y) =an uncorrected estimate of the standard error 
of prediction. 

y =the mean value of the measured attribute. 

y; =measured value of the attribute for experiment 
iv 

Fi =predicted value of the attribute for experiment 
1s 
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Table 7.8 Radiation Count Data for Laboratory Run 1 


EXPERIMENT REDUCED PULP BACKSCATTER CORRECHED 


NUMBER DENSITY COUNTS| COUNTS IRON COUNTS 
Px} E06) PR Sk) 


0.95809 6.117952 one 
D 0.94000 6. 102803 2.115 
3 0.91419 6.086581 9.502 
4 0.89005 6.075626 3.110 
5 0.86558 BOs 7437 4.123 
6 0.85429 6.013170 4.906 
7 0.84192 5.961919 6.193 
8 0.81751 5.884961 eae 
9g 0.80255 5.819540 7.952 
10 0.96104 6.094012 1.274 
11 0.93173 6.080016 2.133 
12 0.90955 6 Os6e73 2.569 
13 0.87998 6.038117 3.651 
14 0.87161 5.991722 4.490 
ales 0.86107 | 5.942789 suilele 
16 0.85212 | 5.889042 6.142 
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Table 7.9 Radiation Count Data for Laboratory Run 2 


EXPERIMENT 


NUMBER 


WN] OWOONMDUBWNH— 


CORRECTED 
IRON COUNTS 
CAT de cue 


BACKSCATTER 
COUNTS 
Pin 6-06} 


REDUCED PULP 
DENSITY COUNTS 


0.84267 6.006070 ne 
0.83124 5.996590 4. 
0.80688 5.991073 4. 
0.78777 5.986919 5. 
0.77289 5.987024 6. 
0.76768 5.954420 6. 
0.76403 5.912414 Z 
0.76903 5.877636 id 
0.76588 5.837851 8. 
0.85075 5.986594 =) 
0.83438 5.982572 4. 
0.81811 5.967861 4. 
0.80539 5.961365 5. 
0.80315 5.926213 5. 
0.79414 5.898715 6. 
0.78765 5.860857 of 
0.83899 5.966444 ms 
0.82373 5.948627 a, 
0.80551 5.948843 7 
0.80103 5.911218 4. 
0.80082 5.882074 5. 
0.84226 5.946631 a 
0.82465 5.938064 A 
0.82538 5.907781 4 
0.84476 5.922986 a 
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Table 7.10 Radiation Count Data for Laboratory Run 3 


EXPERIMENT REDUCED PULP BACKSCATTER CORRECTED 
NUMBER DENS LIY COUNTS] GEUNT:sS IRON COUNTS 
(cet E061) VA (Ee?) 


1 0.9972¢ Me a5 03uS 0., 529 

Zs, 0). 97090 10 .492260 1.696 

3 0.97238 10.419597 2.496 

4 0.96037 10.351724 3.864 

5 0.95048 WOeeo oto 4.790 

6 Ui96839 10.532543 0.803 

7 0.98382 10.462132 2.146 

8 0.96922 10.384149 3.041 

g 0.96012 10.320749 3.943 

10 0.94432 Wo Zorg 09 5.414 
11 0.98750 102506076 12096 
12 O29s0ua. . 10. 442373 1.747 
ts 0.96485 10.349414 3.406 
14 0.95249 W227 3168 4.450 
15 0: GISh6 10.204092 oe she 
16 O33ceS 10.508805 1.424 
vi, 0.97214 10.414876 2.854 
18 0.88250 10.336678 82928 
0.95410 10.236924 D ero 

G. Sock 10.144000 Sy pow 

0.29203 10.472705 1. 608 

QO. SU ag] 10.381948 Pao eo 
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Table 7.11 Radiation Count Data for Laboratory Run 4 


EXPERIMENT REDUCED PULP BACKSCATTER CORRECTED 
NUMBER DENSITY COUNTS}COUNTS TRON COUNTS 
Ox 97E406)) iX40 1 vbe02)) 


.291798 
264365 
216802 
-031330 
92833685 
peeoeue 
. 136001 
a Sholene(ae. 
SS eS 
BPA PASTS 
no09070 
eeg6d2 
SZ 
211968 
UGS 2 
052209 
OT66.16 
. 082540 
. 243785 
Wao ios 
. 086588 
plays (Ecls. 
.630489 
. 144875 
.000495 
5841708. - 
.033064 


— 


—_ —-A 


WO OTM NIM AIM CO VO1O NOW 4AW OD & O10 010) — OO O10 PO 


1 
2 
3 
4 
5 
6 
yer 
8 
fe) 
10 
Eel 
12 
13 
14 
ifs) 
16 
17 
18 


OOO OO @O @O@ 2 OOO OOM 
NW — — 


OOS 1 OS Q'S Ole CeO OO Oe Ce C'o Oa ae 


“4 nom yaodeveted 362 e490 sated aedeaihee re 
al sy apenas 
it 


Sa a ea 
annoo |.) SeRTADS OMe) | RUS eS 
vO9 > (0D. /BSRUGS yTL2mag 
at (G0es°h Alt : 


° —< o ote sf henner 
— 
a~ , o ' 
Tn 5b ‘ eh (iT i 
oot. oe 
“ hha eo © Rem eit 
4 ' “% r 
i 8 Bs CY 3 c iv 
; moth = Me te j 
<9 <4 a | ¥ ) / 
~< 4 a 
> : ic i 
j * | 
Aw ' , 4 
4 es ~~ 5 | 
'« t ( , 
» ra] ' cy S 
Ay ; 
~ ‘" z 
hy 
a# he at we lee j { 
é é ‘a Tie) As { 
14 ‘ei "| « | 
“ + } 
| A ma ] 
} =! i 


AT 
Fe i rye 
ap ‘ oa) J ’ 
r ‘ \/ 
' > at | } ’ } 
Baer 
\ : | 
" 6 \. . ° 
es BhOGa «2 4 
| 2 C #. ' i 5 ; . 
n ~abi 


i 


a - 
—_ “ - os ee ED ee a Mia eS ee i - 
iy 
\ 
; 
P| 
! 
i 
i 7 q 
i. » 
i 
Bail 
Pn 
- 2 
: : ie 
a _ 
va 
we -_ 
‘ 
if be 
2 
4 
i] tn) 


Zila 


performance. The addition of criterion (2), the range error, 
ensures a reasonable spread of calibration data when 
compared to the standard error. The modification herein is 
that the denominator in equation (7.17) does not include a 
correction for the number of parameters estimated from the 
experimental data. This was done to facilitate direct 
comparison of the results obtained with OSCAA to those 
reported by Lyman(32) for the ASHSCAN system. (To account 
for coefficients the 10% figure in equations (7.15) and 
(7.16) should be reduced to approximately 9%: worst case or 


most coefficients.) 


7.6 CALIBRATION STUDIES ON THE LABORATORY DATA 
The calibration models entertained in these studies 
were as follows: 
a) The four component model. 
This model is described by equations (5.76) through 
(5.79). ‘¢ 
b) The three component model. 
This is the model which is the OSCAA analog to the 
ASHSCAN calibration model, where there is no attempt to 
correct for varying iron concentration. The calibration 
equations were a subset of those in (a), with the 
appropriate algebraic modifications. 
c) The Kawatra formulation. 
This was somewhat similar to the form stated by Kawatra 


and Dalton(36) where the % ash attribute is predicted 
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directly from the radiation measurements via an 
empirical model. In this case, however, the general 
Lucas Tooth model (eqn. (5.83)) was employed with a 
stepwise multiple linear regression program(72) for 
parameterization. 

The combined formulation. 

Earlier work by the author on a bench scale version of 
OSCAA had shown that linear expressions, of ene form 
derived from fitting the general Lucas Tooth model, 
could be obtained for direct prediction of both the % 
solids and % ash in the slurry. These models can be 


combined to calculate the % ash in the solids. 


For each model, the measures of performance defined in the 


previous section were computed and reported. Appendix 2 
includes the source listing for the four component model 
calibration computer program and the output for the four 


laboratory runs. 


While the performance criteria were calculated for the 


iron prediction in the four component model they were not 
considered to be relevant in evaluating the overall 
performance of OSCAA. The important process control 
variables are the % ash, and to a lesser extent, the % 
solids, and these were used in determining the system 
Wt but yes 

To minimize the amount of descriptive material 
presented for each of the four laboratory runs, only the 


four component model performance was illustrated in any 
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detail. This model was chosen over the other three because 
it compares favourably with them on a performance basis and 
because it contains a theoretical component, which makes it 
a more sound pyecictie tool. (It is also capable of 
predicting both of the attributes of interest.) It is worthy 
of note that, in the author's opinion, the results of the 
four component model could be improved to some extent if the 
fixed model parameters (section 7.3) were also treated as 
unknowns (cf. Lyman's approach). However, for the purposes 
of this work it was judged unnecessary, and, indeed, the 
goodness of the calibrations represent experimental support 
for the model. 

The four laboratory tests, because of their independent 


nature, are discussed separately below. 


7.6.1 Run 1 Results: Flotation Feed 

For run 1, the data acquisition period was 10 minutes. 
This period was judged short enough to be acceptable For 
process control at McIntyre (see sub-section 3.1.3) and long 
enough to minimize statistical fluctuation in the radiation 
count data. The results of the calibration studies are 
Summarized vinelabless.u2. pigure 7./ 1llustrates the 
predictive accuracy of the four component model for both the 
% solids and % ash attributes. (The + 10% relative error 
envelope is included on the figure to demonstrate the degree 
of precision and distribution of residuals in the 


estimation.) 
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Table 7.12 Summary of Calibration Results for Laboratory Run 


1 


TEST: Flotation Feed (#1) 


Semi-Empirical Models 


RELATIVE ager RANGE ERROR 
is Solids Ash Ss Solids % Ash 


4 Component 
Model Dierili2 


3 Component 
Model Dito 4.27 


Empirical Models 


RELATIVE ERROR RANGE ERROR 
Solids Ash So0l1Lds 3 Ash 


Kawatra > 
Formulation 


Combined ** 
Formulation 


+ } : : : : " ; 
Calibration equation contains 4 "independent" variables 


a= : ; 3 
Calibration equations contain 3 "independent” variables 


cy r : i 


‘> oe 
° aii 


: 2 =“ Faded v by 
‘ a Ps : > a 
tut wretereds? to) eelvead oolzesdsie? ge yiamest oe 


y 
= a ’ ays 
>is a c 
wi abo. Leortiquges 
i vw bos a | —_ 
Sierras 
LOO OE AGEL * A: ——- a =~ —— eS saat 


to <3 | eae ! ae ey .& 


ee oa 


4 
i 
4 = ad 
— te ge = ; 
i pi we 
? * : 
at My cere 7 
yh. «% 
. . a. | 
_— = 12 = ee oe 


: | va _y if fs. : — Mf 


. ot ae 
seidetzev. “snshbneasbeal*- Gaiasne: sos taupe: rolsardaiee a6 


j 
j 
t 


7 


! i 


~ 


‘ca tyra s ‘(pu aha ogo bird e ALANS -e\1G2 ZR Pe. TOU seidpvind re 
= ' = t ‘s | 


ue 44 


4 
~ a 
‘ 
‘A af 

a od a ; f as 

s ~ 

Nv 7 2 er iw 5 

= Ta 

Ley ' 7 


23 


% SOLIDS 


PREDICTED 


oF ELLOr 
Envelope 


MEASURED 
Ze 


PREDIGLTED 


Envelope 


9) 
MEASURED 


Figure 7.7 The Predictive Accuracy of the Four Component 


Model for Laboratory Run 1. 
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A rather extensive residuals analysis was performed on 
the four component model for each calibration run. Since the 
findings were generally the same for all runs, the nature of 
the investigation will be summarized for run 1 only. As 
mentioned above, scatter plots similar to Figure 7.7 provide 
a graphical representation of this information. 

It was mentioned in chapter 6 that size analyses were 
performed on the coal solids obtained from each experiment, 
where quantities permitted. Defining the mass fraction 
passing 0.5mm (35 mesh) as the reference size (reflecting 
coarse particle concentration), the residuals from the 
linear regression analysis on equations (5.74) and (5.75) 
were plotted against size. The range of the size variable 
was Significant, i.e. 68% -0.5mm > 78% -0.5mm. The residuals 
versus Size scatter plots were examined for trends but in 
both cases random behaviour was apparent. As a further check 
on the size dependency of the iron and backscatter systems, 
the residuals associated with the measured and predicted 
attribute values were plotted against size. Again, random 
behaviour wasS apparent. It was concluded from this and 
comparable studies for the other runs that the error 
introduced in calibration as a result of changing particle 
Size distribution was insignificant. 

The errors quoted in Table 7.12 represent an overall 
picture of prototype precision. It is important that these 
values be homogeneous or at least within the acceptable 


limits over the range of the calibration if the device is to 
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be of any practical utility. To study the homogeneity the 
relative errors of prediction for a particular nominal 
solids concentration or a particular nominal ash quality 
"were computed. The results are presented in Table 7.13. It 
is evident from this table that, with the exception of the 
first entries (low attribute value), the errors are 
approximately homogeneous. The apparent problem in the case 
of the low % solids value relates to the performance of the 
density gauge at these levels and will be discussed more 
fully in section 7.6.3. In the case of the low % ash 
computation the magnitude of the relative error is due 
Primarily to one very poor predicted ash value in this set. 
It was concluded from this work that the relative error in 
prediction is approximately homogeneous, or at least within 
the acceptable limits, over the range of the attributes 
studied. 

From Table 7.12 it is clear that all of the models meet 
the performance criteria. With respect to the prediction of 
the ash attribute, the semi-empirical models are observed to 
be slightly more precise. Moreover, the three component 
model exhibits slightly better precision than does the four 
component model. The latter is thought to arise as a result 
of the inclusion of the iron count data in the four 
component model. Since the iron content in the ash is 
relatively stable, this source of information contains more 
noise than it does statistically valid information. 


Nevertheless, the four component model is quite Suitable for 
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Table 7.13 Homogeneity of Relative Error over the Attribute 
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calibration purposes. 

By comparison with Table 7.12, the relative errors on % 
solids and % ash for the ASHSCAN system in similar service 
were 5% and 11% respectively. While these values were 
determined on a plant test scale, the ranges of the two 
variables were about the same. Moreover, since Lyman could 
not employ a unique experiment orthogonal design approach 
under these conditions, there exists the possibilty of. 
multicolinearity?* and/or repeat experiments in his data. 
These would be expected to enhance the performance measures 
to some degree. Finally, Lyman utilizes best fit techniques 
to calculate all of the parameters in his calibration model 
which is also expected to enhance precision. In summary, 
OSCAA compares favourably with ASHSCAN in monitoring a 
flotation feed stream and provides the potential for 
adjusting to changes in ash mineralogy (specifically iron 


concentration) without recalibration. 


7.6.2 Run 2 Results: Flotation Concentrate 

For run 2 the data acquistion period was also 10 
minutes for the reasons given above. Unfortunately, in the 
preparation of the solids samples for the experiments two 
different clean coal bulk samples had to be used. While this 
did not substantially alter the orthogonality of the 
experimental design with respect to the principal 


75Multicolinearity describes the situation where there is a 
high degree of correlation between two or more of the 
independent variables used in designing an experiment or in 
the construction of a linear model. 
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attributes, it did produce a change in the iron 
concentration of the mineral matter. (See Table 7.5 
experiments 1-16 and 17-25.) The net effect was the 
introduction of iron as a pseudo-factor in a non-orthogonal 
manner. For this reason, the nuclear coefficients used in 
the semi-empirical models were taken as the arithmetic 
average of those calculated for the two clean coal samples. 
The results of the calibration studies are summarized in 
Table 7.14. Figure 7.8 illustrates the predictive accuracy 
of the four component model for both the % solids and % ash 
attributes. From Table 7.14 it is clear that, with the 
exception of the three component model, all other 
calibration models meet the performance criteria. However, 
in this case the empirical models exhibit slightly better 
precision with the combined formulation yielding excellent 
results. The failure of the three component model arises 
directly as a result of the eee iron concentration in 
the mineral matter. This illustrates the necessity of 
including iron in the calibration model to preclude frequent 
recalibration work. The comparatively poorer performance of 
the four component model relative to the empirical models is 
thought to be predominantly the result of the compromise in 
calculating the nuclear coefficients. However, the results 
are still quite good when compared to the performance 
criteria and, more importantly, the variation in the 


pseudo-factor was reasonably accommodated. 
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Table 7.14 Results of Calibration Studies for Laboratory Run 
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Figure 7.8 The Predictive Accuracy of the Four Component 


Model for Laboratory Run 2. 


; if ry on 
¥ es yay 
* 14 a! 
: 
UP = 
to ns 
i ay 
: sly 
i = 
=) 
' ny a 
2 bh ort 
7) 
i 
i = 
= — ~ a ee ?_ —- . - Ae et a AD A LS A ca 
- = oe tb 
- ey - in 
* 
* : é —< 
% ae 
™* <2 
& 
: oe 
* ~~ 7 
+ 
“AS 
~ he ™a, 
a 
= ~ » 
' =) 
™ * ~~ 
* ‘ " 
ie a * 
ae ay ? 
, ” td ba 
8 . i , 
in. * 


ea 
Be 
~ 
1 

hee ment — he i —— - —_ ———iprom- re 

rs >» Avan 

Stel “Lit SA- @ 
Ra 
——— a ee aa — oe apt mehr ee =e) MA. ne 


— 


pea a 


+oOx3e att pall ~ ™ ° . 


A ™ ~ 
Sie Loud « % be : ~ 
5 A ‘ x ~ 
7 * ~ 
Me 
2 ; Ny — 
c ~ — —— en Tees me 1° ae Vn 


 gaireaas | ae 
ia. 
41s) roqgnes 176% ads oe Yawavno4 avizoiteee oat 8. a 


soup” 
AN ba Me; . . Ss wa xvoee ones eel 


papers i =) : Ye Dis eo ‘| rae _ 


223 


Again, by comparison, the relative errors on % solids 
and % ash for the ASHSCAN device in similar service were 3% 
and 3.8% respectively, although the ranges were slightly 
_ larger than those in run 2. The same comments on 
experimental integrity hold. In summary, OSCAA compares 
favourably to ASHSCAN in monitoring a flotation concentrate 
Stream and has demonstrated its ability to account for 
changes in ash mineralogy, something of which ASHSCAN is 


incapable. 


7.6.3 Run 3 Results: Flotation Tailing 

A considerable amount of time (=7 months) elapsed 
between the completion of run 2 and the commencement of run 
3. As a result of the relatively short half life of the 
Gd-153 radioisotope (242 days) the data acquisition period 
was doubled to 20 minutes. This provided for approximately 
the same number of counts in the pulp density sub-system as 
fornmsimilan conditions in run 2, but had the effect of 
improving the statistical accuracy of the backscatter/iron 
count data because of the longer half life of the Cd-109 
radioisotope. The results of the calibration studies are 
presented in Table 7.15. Figure 7.9 illustrates the 
predictive accuracy of the four component model for both the 
% solids and % ash attributes. 

BLOM able: i slogiteis cléar that none? of) the 
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Table 7.15 Results of the Calibration Studies for Laboratory 


Run 3 


TESS: PlLOovarron-Tarlingmro).. All Data 


Semi-Empirical Models 


RELATIVE ERROR RANGE ERROR 
MODEL S Solids % Ash ZeSOLEGS 2 aS 
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46.87 oe 147.7 
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Figure 7.9 The Predictive Accuracy of the Four Component 


Model for Laboratory Run 3. 
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solids concentrations and would thus be prone to large 
errors in prediction based on radiation statistics alone. 
This was demonstrated using the theoretical three component 
‘model. Equation (5.1) characterizes the pulp density 
sub-system and equation (5.58) characterizes the ash 
analysis sub-system. Using the theory and assuming a typical 
flotation feed pulp composition to estimate the other 
unknowns in the equations, I» and No were calculated to 
provide total detected counts (excluding background) 
comparable to those measured in run 1; i.e. ~400,000 and 
~6,000,000 respectively. Using Monte Carlo techniques a set 
of simulations were performed at different levels of % 
solids and % ash, and for each the relative error of 
prediction for each attribute was calculated. The results of 
the simulations are summarized in graphical form in Figures 
7o10gand 7.11 a krom the plots, it 1s apparent that as % 
solids decreases the relative errors increase, this being 
particularly evident at low levels of the former. Moreover, 
while the relative error on % solids is effectively 
independent of the ash content, the relative error on % ash 
does exhibit a slight dependence on the attribute value. 
Calculations at very low solids concentrations, ~1%, were 
unstable as expected, however, the relative errors appeared 
to be in excess of 70% in all cases. This is consistent with 
the trend of both relative error curves at low solids 
concentrations. In this light, the values of the performance 


criteria reported in Table 7.15 seemed quite reasonable. 
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RELATIVE ERROR ON % SOLIDS PREDICTION (3) 


% SOLIDS 


Figure 7.10 Relative Error on % Solids Prediction as a 


Function of the Attribute Levels. 
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ASH PREDICTION (3) 
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Figure 7.11 Relative Error on % Ash Prediction as a Function 


of the Attribute Levels. 
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To partially remove this statistical effect in 
calibration the data set was modified by excluding all 
design points at 1% and 2% solids. The modified data set was 
then re-analysed and the results are summarized in Table 
7.16 and Figure 7.12. From Table 7.16 it iS apparent that, 
with the exception of the three component model, the 
relative error for % ash is met in all cases but the range 
error is unacceptable. The latter is thought to be the 
combined effect of statistical fluctuation coupled with 
Sampling and analysis errors. The % solids predictions fail 
both performance criteria due to statistical fluctuations 
together with the low mean value of the attribute and the 
Narrow range of study. In examining the results for the 
three component model it seems evident that, in this case, 
the iron count information is necessary to stabilize the 
results. 

Considering these results it is recommended that for 
monitoring a flotation tailings stream of low solids 
concentration, either extended count times would be 
required, or preferably, a stronger source be employed to 
improve upon the statistics. Monte Carlo techniques could be 
used to obtain a, firstwapproximation on “Source Strength 
uSing the method described above. 

For the sake of completeness, the relative error on the 
% solids and % ash for the ASHSCAN system in similar service 
were 9% and 7% respectively. While the comments on 


experimental integrity are valid here as well, an important 
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Table 7.16 Results of Calibration Studies for the Modified 


Data Set in Laboratory Run 3 


TRST: Flotation Tailing (#3): Reduced Data Set 


Semi-Empirical Models 


RELATIVE a ak RANGE eee 
ae Solids Ash Solids Ash 


Component de erat’, LU a34 
as 
3 Component 
Model oe. 1 7] Ole 209.20 24 


Empirical Models 


RELATIVE peo RANGE ee 

MODEL Solids Ash ee Ash 
+- 

Kawatra 

Combined ** 

Formulation 12.54 2070147 


+ 
Galibration equation contains 1 "independent variables 


aa : one ‘ 
@aliibracion equations contain 3 “independent! variables 


) CVS eal a) 
bel 43 BOM sty tak" ieee nol Sane Lae esis at 


} . 
toe Ase Fans 2 § 
te ~ewe Bh 
~— 
oe os oo — — 2 eee. ead ong anne oe eh 


57 ‘eau BW) 209% “SVY CALE 
ied Fant ces .- thiice # 


i 
e 
= a ~ 
; i” \ t. yt i ee 
‘ i % 7 
—o yetewe rem — oe se lc A ALLO AA ALLE LLL 
4 
“ ~ 
‘ 
i 
e * 
Soe 
2 = 


Si) We | RAINS ers 

fis ea kia ROBE LOS 8 

- ee ee a, ee Lealandesl ee a. 
eta | a5 | 


ea .eS Meio 2... eee { Aus 


— —— 1 oe 


‘ ) Ts 


7 Ne ve 


rm ay 4 ‘i 5 hyaSqaba ** =f at) a + — 4 ae | 6 ry Ac serdte = Ts 


‘;e i. 
ze 


We 

aL 

~) 
3 


2 4 
> iv 4 od 


Se1CE628Y, 3h hisce efnd*® & niéeinay e6eg2nupe abieewtdsese 


j ~ i 7 
? af cn = 
is ; 
‘et Fi 
: i] j a4 
iat 
i. - : : 
2 j ® _ ' y af i 7 
i] v 
; ) é 7 i, are a 
a] “ae a4 ; 
ae ’ ir. y 
Vi ’ a 
Af : ~ : ’ — 
Ae Ni 


251 


5 
% SOLIDS 
fa) 
i 
a 
‘) 
= 
Q 
| 
[oe 
A, 
t+ 10% Error 
Envelope 
2 
MEASURED Z 
49 
Q 
fx 
a 
Ss) 
He 
Qa 
ica} 
~ 
Ay 
Envelope 
29 


29 
MEASURED oe 


Figure 7.12 The Predictive Accuracy of the Four Component 


Model for Laboratory Run 3:Modified Data. 
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difference was the range of study. The range of solids 
concentration for ASHSCAN more closely resembled that of the 
flotation feed stream in run 1. This would be expected to 
Stabilize density gauge performance and hence the overall 
performance of the OSCAA system. In summary, OSCAA, in its 
present form, is not well suited to monitoring flotation 
tailings streams at low solids concentrations. However, 
increasing the statistical accuracy of the pulp density 
Sub-system together with the minimum value for the solids 
concentration would be expected to produce results 


comparable to those reported for ASHSCAN. 


7.6.4 Run 4 Results: Flotation Feed with Varying Iron 

Concentration 

The data acquisition period for this run was 20 minutes 
for the reasons given above. The material used to control 
the iron content of the ash was fine magnetite, not hematite 
as assumed in the four component model. This was done to 
reflect contamination by tramp heavy medium which can occur 
in those plants which employ heavy medium washing circuits. 
The results of the calibration studies are summarized in 
Table 7.17. Figure 7.13 illustrates the predictive accuracy 
of the four component model for both the % solids and % ash 
attributes. 

From Table 7.17 it is apparent that, with the exception 
of the three component model, as expected, the other 


calibration models met or were very close to the performance 
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Table 7.17 Results of Calibration Studies for Laboratory Run 


4 


REST: Flotation Feed with Varying Iron (#4) 


Semi-Empirical Models 
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Figure 7.13 The Predictive Accuracy of the Four Component 


Model for Laboratory Run 4. 
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criteria. The comparatively poorer performance in % solids 
prediction for this run when compared to run 1 is due, in 
large part, to the unexplainable difference in pulp density 
Sub-system performance (see density gauge calibration 
statistics in appendix 2). Never-the-less, the results for 
the four component model make it an acceptable predictive 
BOOL, 


There was no similar test for the ASHSCAN system. 


7.7 SUMMARY OF LABORATORY INVESTIGATIONS 

From the discussion above it is evident that an OSCAA 
type system would be feasible for process control purposes 
Meche fine coal circuiteat Mcintyre Mines. Ltd. The tailings 
stream would have to be excluded from the analysis based on 
current prototype design. Rather simple design changes would 
be expected to correct this problem and allow all three 
process streams to be monitored. Furthermore, it is expected 
that with improvements in design (specifically more counts 
in a shorter time), the data acquisition period could be 
reduced to 5 minutes, which is more in line with the process 
response time. 

It is apparent from run 2 and run 4 where the iron 
concentration of the mineral matter was a factor that a 
three component slurry model was not satisfactory for 
attribute prediction. This is important as evidenced in this 
project where two clean coals from the same deposit provided 


Significantly different iron values, and in Lyman's work 
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where he observed that the coal could have an iron 
concentration in the ash in the range of; 5% < [Fe] < Sy. 
This would require frequent instrument recalibration or at 
‘the very least frequent updating of the calibration model 
parameters. This kind of system requirememt is incompatible 


with the those of the industry. 
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8. THE OSCAA SYSTEM PLANT CALIBRATION RESULTS 


8.1 INTRODUCTION 

Upon completion of the laboratory calibration runs it 
was decided to proceed with a plant trial of the OSCAA 
prototype device. McIntyre Mines Ltd. personnel agreed in 
principle to cooperate, pending the finalization of the 
experimental plan. 

Based on the laboratory results and the difficulties 
attendant to flow stream sampling and multiplexing, a sooo 
flow stream was chosen for monitoring, the flotation 
concentrate. For practical reasons it was necessary to 
choose the concentrate stream from the number four flotation 
cell bank. With this general objective in mind a 
reconnaissance trip was made to the minesite for preliminary 
design purposes. Subsequently, an experimental plan was 
developed and the responsibilities of the operating 
personnel clearly defined. McIntyre agreed to the proposal 
and installed the necessary piping and valves required for 
primary sampling. In concert with this work the balance of 
the materials and equipment were prepared in the University 
laboratories. One of the paramount concerns was to provide 
adequate housing for the instrumentation since no facility 
waS available on-site. Toward this end a small trailer was 
rented, and the required modifications performed including 
the installation of an air conditioning unit. Provision of 


the latter is standard practice for industrial on-stream 
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analytical instrumentation to control airborne dust and the 
ambient temperature. On June 8, 1981 the OSCAA prototype was 


transported to the minesite. 


8.2 EXPERIMENTAL PROCEDURE 
The experimental procedure may be divided into three 
parts, equipment set-up, data acquisition and sampling and 


analysis. These are discussed below. 


8.2.1 Equipment Set-Up 

For a period of 5 days after arrival on the minesite 
the OSCAA prototype was prepared for the plant calibration 
work. This in¢luded™tratlerelocation; anstaliation of the 
electronic equipment in the trailer, plumbing the sampling 
System, installing the power lines and situating the 
analytical facilities. 

Plates 8.1 and 8.2 illustrate the installation. The 
primary sample was obtained from the number four flotation 
cell bank concentrate collection box and flowed by gravity 
to the primary sample pump?‘. The sample was then pumped to 
the secondary sampler where a representative sub-sample was 
cut for analysis by the OSCAA unit. The reject sample from 
the secondary sampler and the discharge from the OSCAA 
sample line were combined in the collection launder. This 
launder was designed such that either open or closed circuit 


7*Pipeline sizing for the primary sample, which would ensure 
sufficient flow while minimizing the need for valve 
correction was done by solving the appropriate macroscopic 
isothermal energy balance. 
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Plate 8.1 An Illustration of the On-Site OSCAA Equipment 
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Plate 8.2 The OSCAA Measurement Sub-Systems for the 


Treva l 


The ash analysis and pulp 
density measurement sub- 
systems of the OSCAA proto- 
type, as installed for the 
plantGeenaadl. 
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sample flow was possible. In the open circuit mode the 
launder discharged to the primary sample return line which 
flowed by gravity to the main sump in the preparation plant 
' basement. Ultimately the primary sample was recirculated to 
the preparation plant feed. This recirculating load was 
negligible. 

Once the flow loop had been fabricated and tested with 
water, a batch calibration was attempted using solid samples 
prepared in the University laboratories. A two factor (% 
ash, % solids) three level orthogonal design was used which 
wasS intended to provide estimates of the calibration 
equation parameters. It was hoped that by using these 
estimates it would be possible to provide the operators with 
on-stream analyses of the concentrate at a relatively high 
frequency... Thisswould demonstrate) the, utalitysofysuch.a 
device in daily operation. The resultant count data was 
consistent with previous laboratory findings, however, the 
Samples which were required for attribute analysis were 
accidentally discarded by a member of the McIntyre quality 
control laboratory staff. As a result, it was impossible to 
proceed with the demonstration which had been planned. 
Notwithstanding this unfortunate circumstance, the equipment 
appeared to have operated satisfactorily although the 
ambient conditions at the time of this calibration work were 
wet and cold, thereby minimizing dust and temperature 


problems. 
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One point which is worthy of note here is that it was 
impossible to obtain a process water standard for data 
reduction purposes. The only water supply available was the 
fresh water make-up which is pumped to the plant from the 
nearby Smoky river. This water was very highly aerated and 
Silty. In practice this source contributes only a small 
portion to the process water requirements, the bulk 
consisting of water recirculated from the closed plant 


dewatering circuit. 


8.2.2 Data Acquistion 

The data acquisition period was defined to be twenty 
minutes which was consistent with laboratory calibration 
runs 3 and 4. The computer programs controlling data 
acquisition and data processing had to be modified to 
operate on a continuous basis. Because of the limited fast 
access memory of the computer an extensive amount of tape 
transfer was involved to allow the two programs to work in 
sequence. Output from the programs included the time of day, 
gross density counts and the corrected backscatter and iron 
fluorescent radiation counts. 

The objective of the plant trial was to run as long as 
was possible subject to two constraints, namely; 
a) the maximum number of samples to be processed should be 

6nhPthe ordér tof 42009 +and, 
b) rental equipment was to be returned by the end of June. 


Since the operation of the equipment was the responsibility 
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of the author, daily shift maximum operating periods of from 
8 to 12 hours were defined. Moreover, since the biggest 
changes in attribute values were expected to occur during 
‘day shift, experimental work was to take place in the 
interval of 8:00 am to 8:00 pm. Operating under these 
constraints together with scheduled and unscheduled plant 
Operating interruptions, calibration work was only possible 


on five operating days. 


8.2.3 Sampling and Analysis 

To allow one man to adequately monitor the OSCAA 
system, obtain the samples and perform the first steps in 
their analysis, it was necessary to resort to grab sampling 
over the more accurate composite sampling. On the basis of 
the fluctuations in the concentrate stream attributes 
observed in Figures 3.4 and 3.6 this was deemed an 
acceptable sampling methodology. Furthermore, since the flow 
Stream is relatively high in slurry solids content, only 1 
litre of sample was taken from the OSCAA sample line. (While 
the method was thought to be acceptable it should be noted 
that the representativity of such a sample is, on the 
average, about 1 part in 2650, cut at one instant in time. 
It will be shown later that this method appeared to have 
contributed significantly to the poor measured performance 
of the pulp density gauge.) 

Samples were cut at approximately the midpoint of the 


corresponding data acquisition period. In some cases, 
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Sampling was done at both the midpoint and at the end of 
this interval as a check on the assumptions in the sampling 
methodology. Duplicate samples were also cut to provide 
estimates of sampling and analysis error. 

Immediately after a sample was cut its wet weight was 
recorded and the slurry filtered using a portable, vacuum 
pump driven, Buchner apparatus. The moist filter cakes were 
bagged, labelled and packaged for shipment to the University 
laboratories. where they were subsequently dried and 
weighed. Size analyses were not performed on these samples. 
A representative 50 g sub-sample was pulverized to -0.15 mm 
(100 mesh) and duplicate 3 g samples extracted for ash 
analysis. On the basis of cost and the fact that the iron 
content of the ash appeared to vary negligibly over the 
calibration trial, iron analyses were excluded. The latter 
meant that the four component semi-empirical calibration 
equation including the iron correction could not be tested 


on the plant data. 


8.3 PROBLEMS ENCOUNTERED IN THE PLANT TRIAL 

On the first day of the full calibration testwork the 
weather was sunny and warm and it quickly became apparent 
that the air conditioning unit was faulty, although it had 
successfully completed trial runs at the University. The 
electrical demand of the air conditioner on start-up was 
determined to be 35 amps, which, because all instrumentation 


waS On a common circuit, caused some of the nucleonic 
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instruments to fail. After consultation with the plant 
electrical department and reference to the appropriate 
operating manuals it was decided to continue the work 
without providing air conditioning. It waS arranged to 
record the temperature inside the trailer and make periodic 
background measurements with the pulp density sub-system 
(closed shutter). This data would be used in case some 
instrumentation correction was required although it was not 
expected. 

Another problem which arose was the result of poor 
communications between plant personnel. The author was to be 
notified in advance of any scheduled plant shutdowns to 
facilitate calibration work re-scheduling. Unfortunately, on 
two occasions, this information was not relayed and, as a 
result, considerable time was wasted. It should also be 
noted that at the time of the plant trial, McIntyre was 
somewhat behind on production commitments and the 
preparation plant had been running with minimal planned 
preventive maintenance. This resulted in the plant being in 
a sub-normal state of repair and frequent short (=1 h) 
emergency stoppages were experienced. 

The calibration work was prematurely halted when it was 
observed that the density counts obtained on recirculating 
water were virtually the same as those recorded for the 
concentrate slurry. This suggested that there was a serious 
problem in the pulp density measurement (e.g. air 


entrainment) and effort was then taken to study this. The 
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net result was that good?” data was obtained for only three 


Operating days. 


8.4 DATA PREPARATION 
The data preparation was divided into attribute 
analyses and radiation count handling, similar to the 


laboratory calibration runs. 


8.4.1 Preparation of the Slurry Attribute Data 

While attribute analyses were done for all of the 
samples taken during the plant trial, only those 
corresponding to the calibration period are presented. In an 
effort to obtain better estimates of the> attribute value 
ever the 20 minute period than the point estimate provided 
by the grab sample, a form of data smoothing was employed. 
The method was to fit a quadratic to the sample point under 
consideration and its two nearest neighbours, in time, on 
either side. From the quadratic, using analytical 
integration, a mean attribute value was calculated. The 
procedure is illustrated in Figure 8.1. In cases where one 
of the neighbours was sufficiently distant (>0.5 h) so as to 
be potentially unrelated, it was replaced with a dummy point 
having the same value as the centre point and 0.5 h removed 
from it. The results of the attribute analyses and the 


smoothed data are presented in Tables 8.1, 8.2 and 8.3 and 


27In this context "good" implies complete in the sense that 
temperature, density background counts and resonable 
operating times were obtained. 
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the corresponding figures 8.2, 8.3 and 8.4. The effect of 
the smoothing algorithm is clearly illustrated in the 
figures where the smoothed data are connected by a solid 
‘line. These figures were included to provide further 
evidence that the trends in Figures 3.4 and 3.6 were typical 
of process behaviour (i.e. process fluctuations appear 
amenable to process control). 

Not all of the data presented in the tables were used 
in the calibration equation parameterization exercise. The 


reasons for this are given in the next section. 


8.4.2 Preparation of the Radiation Count Data 

Since a process water standard was not available it was 
impossible to reduce the pulp density counts and correct the 
iron fluorescence counts in the manner described in the 
previous chapter’*. Thus, the data is presented in the form 
where only background corrections have been made. 
Furthermore, data which was clearly abnormal was excluded 
from the calibration data set. Abnormal data was identified 
by unusually high pulp density counts and correspondingly 
low backscatter counts. This was the result of a temporary 
interruption in the primary sample flow generally 
attributable to an emergency plant stoppage. The calibration 
data set is presented in Table 8.4. The temperature data 
included in this table were generated from point estimates 


7®Since the thin mylar windows were replaced each day, as a 
precautionary measure, the absence of standardization would 
be expected to allow the effect of bias due to window 
positioning to enter the calibration data. 
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Table 8.2 Attribute Analyses for June 20, 1981 
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Table 8.3 Attribute Analyses for June 22, 1981 
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Figure 8.2 Attribute Analyses for June 19, 1981. 
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Figure 8.3 Attribute Analyses for June 20, 
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Figure 8.4 Attribute Analyses for June 22, 
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Table 8.4 Radiation Count Data for The Plant Calibration Run 


* Solids} % Ash Density Backscatter Tempe: 

Counts Counts Deg.F 

(1xE-5) (1xE-6) 

31.15 7224 4.57343 10.388398 6.3736 92.6 

32.23 6.81 4.51423 10.408299 7.0382 93.2 

31.50 7.08 4.73937 10.351318 7.0461 95.4 

34.47 6.37 4.69215 10.487373 7.01469 96.6 

S202 7-68 4.72298 10.355154 6.3174 97.4 

32.90 7-42 4.68882 10.4446918 6.5315 98.1 

31.58 8.46 4.69010 10.340450 5.7865 98.8 

29.07 10.77 4.96664 10.112243 LALIT LI/ 99.7 

26.80 11.58 4.95584 10.108125 765554 100.3 

26.2 11.56 3.02029 10.111031 7 +2696 POT? 

eae ee 11.69 DlOeeeny 10.114601 6.8410 101.8 

25.64 10.71 53-01043 10.1473506 6.92495 102. 

24.76 10.46 3+044642 10.1468697 6.9175 102. 

26.67 10.16 3+14163 10.157407 46-8279 103.4 

24.86 Heder 3+17278 10.1018465 6.8933 103.7 

27.76 11.70 Deore 9.994504 7+4028 103.8 

23.35 DS cares IeOds es 9.962301 7.7018 104.1 

27.42 11.88 3.31728 10.018292 769339 104.6 

Bbsyes UG? 11.90 9+27665 10.007783 7.6978 104.7 

27.22 10.90 5.38057 10.072563 7.8258 104.5 

30.13 Gee CAE 5-24708 10.193346 7+3009 103.5 
32.16 Pare & =] 93-15765 10.220228 7A 103.3 

Lil CAs, 9.24 eG, 10.242427 7+4754 103.6 
24.75 11.74 Siem 10.088490 7 +8554 103.5 
26.99 10.90 Deleon 10.115785 73053 103.7 
A os 9.50 53-04081 10.143489 7 +3243 103.35 

25.49 10.50 3+07205 10.139630 7.1547 102. 
2 24626 10.94 3-00054 10.145406 7+1506 102.3 
29 35.00 Tee 4.465055 10.412752 TAS RVR 93.8 
30 35.06 6.98 4.58692 10.404745 669329 94.0 
31 35-34 6.70 4.60629 10.451601 &-8009 94.5 
Se 32.00 LOS 2 4.61127 10.179662 7 +6086 96.2 
33 32.40 10.38 4.62596 10.165187 7-5892 96.7 
34 32.78 9.81 4.63606 10.242674 7.0897 Tiel 
35 Slice 9.44 4.73058 10.247135 7.0242 Diet 
36 30.24 10.83 4.71622 10.120563 7 - 6383 98.5 
37 29.89 shale Zhe) 4.73850 10.10185 7+3852 OSS 
38 30.47 10.82 4.85252 10.171263 7 +3063 100.8 
39 34.34 8.64 4.88778 10.330417 6.9058 101.4 
40 35.41 7Al 4.95682 10.432425 6.3553 101.7 
41 30.77 8.14 53-0421\9 10.457833 $4133 100.7 
42 34.86 6.49 4.99100 10.478192 6.6882 100.3 
43 34.41 7-00 4.99742 10.471340 6.4082 100.0 
44 33.20 6.93 4.31606 10.449181 6.8990 87.7 
45 32.66 6-90 4.21734 10.428157 6.9321 88.2 
46 30.43 Cae, 4.31206 10.197273 We Boel 89.8 
47 30.16 yo 4.35408 10.238573 Bg Set 90.6 
48 29.84 isciizg 4.37133 10.204775 7-2288 vil aes 
49 30.39 9.20 4.36165 10.241688 7.0615 91.8 
50 30.16 9.38 4.38158 10.196981 6.9368 92.8 
myst 29.74 9.09 4.41910 10.140953 Pile 94.5 
S2 27a 9.34 4.64718 10.229208 6.6122 SUE 
5 24.24 9.68 53+096700 G6, 215297 6.7017 100.9 
-54 so Fat § Ps 4S 53701687 10,241673 7s bh Oe SOA ay 
aS 29.29 8.98 5.06190 10.272856 7.0417 TO2R2 
36 30.45 8.98 4.98296 10.265250 7.0721 103.1 
S 30.02 OC AES 5-05247 10.194013 6.72935 103.2 
38 AM de FAs 10.66 3.05342 10.142391 7.2761 102.6 
=? Pada ges Peal 5.14169 10.183740 74153 102.4 
40 29.83 9.86 5.04899 10.1929860 Pee 102.0 
61 29.85 9.46 5.13097 10.204108 7.0594 101.9 
62 29.83 9.87 Gp Les LO 1 9LS77 7+1490 101.7 
63 PAROLES 10.01 5.19442 10.180995 7.1969 102. 
44 29.83 10.04 3+05396 10.185493 7 pleO 102. 
65 29.69 ORS 5.09294 10.170326 7.3074 LOSeeh 
56 30.04 10.31 Deke 2 JO 1O/F17 Gad 2 Ts 2IAOO 103.6 
67 29.07 10.44 5.14008 10.169977 Zeno oe 103.7 
68 Petia Wer) 10.12 3.01636 10.193524 7.90982 103.0 
290356 10.13 4.98099 10.190334 7.0374 102.9 
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and the quadratic smoothing algorithm described in the 


previous section. 


8.5 PLANT CALIBRATION RESULTS 

It was evident in conducting this testwork that the 
electronics section of the pulp density sub-system exhibited 
a strong dependence on ambient temperature. The extent of 
this influence was not expected on the basis of the material 
contained in the operating manuals for this equipment. (This 
problem was independently addressed in subsequent testwork, 
the results of which are discussed later.) The net effect is 
that an accurate pulp density count was not available and 
the three component semi-empirical calibration model could 
not be tested on the data set without rather radical 
modification. However, the Kawatra formulation of an 
empirical calibration equation was directly tested and the 
results are presented in Table 8.5 with an illustration of 
predictive accuracy given in Figure 8.5. The calibration 


equation had the form shown below: 


3 Ash=b +b, (BACKSCATTER COUNTS) +b 


0 (BACKSCATTER COUNTS) - COR) 


2 (TEMPERATURE) 


The contribution of the cross product term is only barely 
Significant and its exclusion has little impact on the 
evaluation criteria (which become 4.47% and 6.21% 
respectively.) The presence of temperature in the 
calibration equation may indicate a slight temperature 


dependence in the ash analysis sub-system electronics. The 
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Table 8.5 Results of the Plant Calibration Studies 


TEST: IN-PLANT CALIBRATION STUDIES 


Semi-Empirical Model ~ 


RELATIVE ERROR (3%) RANGE ERROR (3) 
$ Solids 


RELATIVE ERROR (3%) 


RANGE ERROR (3) 
% Solids % Ash 
Kawatra * 
ay ake N/A dcaly N/A 5.80 
Formulation 


+ Calibration equation contains 2 "independent" 
variables. 
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accuracy provided by this calibration indicates that the 
OSCAA unit has definite potential for application as an 
industrial process control sensing element. 

In connection with the simulation work discussed in 
sub-section 7.6.3 (ref. Figures 7.10 & 7.11), a Separate set 
of Monte Carlo simulations were performed where the error in 
obtaining the pulp density count, e=/VIo, was multiplied by 
an integer valued variable, m, i.e. me, with m=1,2,...,5, 
prior to perturbing the expected counts for computation of 
the attribute values (m=1 corresponds to the Poisson 
distribution). It was evident from these simulation 
calculations that, while the relative error on % solids was 
Slightly more sensitive to an increase in the pulp density 
measurement error than was the relative error on % ash, both 
remained within acceptable limits. Accepting that the 
prediction of % ash is therefore rather insensitive to the 
pulp density, the three component model was modified and the 
calibration performed, the results of which are included in 
Table 8.5. The modification was to fix the reduced pulp 
density counts at a level corresponding to the mean slurry 
composition. It is readily apparent from Table 8.5 that, in 
so doing, the % solids prediction fails to meet the range 
error criterion as expected, however, the % ash prediction 
criteria are well within the acceptable limits and do not 
differ substantially from those quoted for the empirical 
model. A plot of the % ash predictions using the modified 


semi-empirical model strongly resembles Figure 8.5 and thus 
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has been excluded. 

Considering the accuracy of the the ash analyses 
(empirically smoothed point estimates) the evaluation 
“criteria are remarkably good. One might not expect such 
impressive results for the flotation tailings and feed 
streams, however, for initial control activities the 
concentrate grade is the most important control variable. In 
the opinion of the author, one of the principal reasons that 
a pulp density measurement was not required in the 
calibration equation was that, in practice, pulp density and 
ash in the concentrate stream are broadly correlated (see 
also (74) pg. 234). Figure 8.6 illustrates this correlation 
for the attribute analyses given earlier. Such colinearity 
works in favour of the empirical model which needs only a 
Single meaSurement under these circumstances. 

It must be kept in mind that, for the period of 
testing, the mineralogical composition of the solids species 
appeared to be relatively constant and, as such, should 


provide for the best calibration results. 


8.6 TEMPERATURE DEPENDENCE OF THE SCA 

Upon the return of the equipment from the plant trial 
the temperature dependence of the pulp density Bebe pecan 
was further investigated. Initially, with water in the 
density gauge flow line, the various components of the 
Sub-system were heated with a hot air gun and the effect on 


count rate observed. The only component displaying any 
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temperature sensitivity was the single channel analyser 
(SCA). 

The NIM bin, containing the SCA, was placed in a glove 
box with a heater and a temperature controller. Under these 
conditions the average count rate through water was obtained 
as a function of the ambient temperature. The results of 
this experiment are shown in Figure 8.7. A linear 
temperature dependence was observed. 

Having demonstrated the temperature dependence of the 
SCA output an effort was made to remove the associated bias 
from the observed pulp density counts. A model of the form 
given in equation (8.2) was fitted to the appropriate data 


from Table 8.4. A 


C)=Dexp (-Cp) +AT+B 
ees (BZ) 


where: Ch =predicted observed pulp density counts. 

D =theoretical incident photon intensity (cf. 

Lambert's Law). 

C =the density attenuation coefficient (=ux). 

A =temperature coefficient. 

B =background counts. 
The constant B in the equation was estimated from closed 
shutter measurements and was fixed for the purposes of 
parameterization. The value of p was calculated from the ash 
and solids content assuming a two phase solid (ref. Chap. 7) 
in water. Using an unweighted least squares objective 
function and a multivariable search algorithm(75), estimates 
of the parameters D, C and A were obtained. The values of C 
and D compared closely to those obtained in the laboratory 


experimentation. The index of multiple determination for the 


fit was 0.90 and the temperature effect was observed to make 
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the biggest contribution to the explicable variance in the 
observed counts. The contribution from the term in equation 
(8.2) representing Lambert's law was substantially lower. 
Nevertheless, the "expected observed" counts (temperature 
bias and background removed) were generated from this term 
and a modified data set prepared for analysis. Both the 
empirical (Kawatra) and semi-empirical (three component- 
unmodified) model were tested on this data set. For the 
former, no change in model form was discovered. For the 
latter, the density correlation remained extremely poor and 
the model was incapable of meeting the performance criteria. 
The conclusion which can be drawn from this exercise is 
that the pulp density measurement was strongly correlated 
with the ambient temperature of the SCA. However, even when 
the effects of this bias were removed, there was essentially 
no correlation between pulp density and the expected 
radiation counts from the measurement sub-system. While it 
was impossible to absolutely quantify the contributing 
factors it was strongly suspected that poor estimates of 
Slurry attribute data over the measurement interval, 
inadequate source strength/count time, entrained air and 
electronic noise were the dominant causes of poor sub-system 
performance. This aspect is treated in more depth in 


sections 8.7 and 8:8. 
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8.7 ASSOCIATED PLANT TESTWORK 
There were several other tests which were conducted in 
conjunction with the plant OSCAA calibration trial. These 


are discussed below. 


8.7.1 Estimation of Sampling and Analysis Error 

The results of the calibration have shown that the 
Sampling and analysis techniques provided reasonable point 
estimates. However, quantitative evidence was gathered to 
demonstrate the reliabilty of the samples. 

In order to deduce the error associated with the 
sampling and subsequent analysis of the material obtained 
from the OSCAA sample line during the calibration work, 
duplicate cuts were made for the final 8 calibration points 
on June 22, 1981. Assuming that these errors are normally 
distributed and that the variances are homogeneous, the mean 
attribute value, pooled variance, standard deviation and 


relative error of meaSurement are as follows: 


% SOLIDS 
Mean Value = 29.7% 
Pooled Variance = 0.0144%? 
Standard Deviation = 0.12% 
Relative Error = 0.4% 

% ASH 
Mean Value = 10.14% 


Pooled Variance = 0.004%? 
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Standard Deviation = 0.063% 

Relative Error = 0.6% 
The relative errors are considered to be quite good 
indicating consistency in sampling technique and the 


Subsequent sample preparation and analysis. 


8.7.2 Flotation Concentrate Sample Integrity Check 

To check on the representativity of the sample 
delivered to the OSCAA prototype, relative to the slurry in 
the flotation concentrate collection box, a sampling 
campaign was conducted. Composite samples were obtained (6 
cuts at 5 min. intervals) at three points in the flow system 
as shown in.Figure 8.8. The analytical results for the three 
points are included in the figure. 

Inspection of the sample attributes indicates good 
agreement between the reject and sample lines in the OSCAA 
System. However, an apparent bias exists for the solids 
concentration and particle size Beer ree between the 
collection box sample and the other two. The differences are 
expected because of the flow system design and the sampling 
technique. It is impossible to make a case for either sample 
1 or sample 2/3 being "more correct" and, in fact, the true 
values of the flotation concentrate attributes are pEobably 


some intermediate of those reported. 
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Figure 8.8 Concentrate Sample Integrity Check. 
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8.7.3 Preliminary Studies on Flotation Concentrate Air 

Entrainment 

It was noted earlier that calibration sampling was 
suspended when it was thought that aeration was playing a 
Significant and variable role in pulp density measurement. 
In order to try and study this effect an entire day was 
spent monitoring the concentrate stream attributes, making 
measurements with the OSCAA system over short intervals (12 
min.) and varying the frother concentration in the flotation 
feed. The effect of increasing frother concentration was 
expected to stabilize the froth bubbles. It was thought that 
in so doing more air would be entrained in the primary 
Sample and hence produce interference in the pulp density 
sub-system. Substantial changes in the frother addition rate 
were made (minimum=40cc/min; maximum=180cc/min) and 
essentially no changes were observed in the pulp density 
counts. Although these results are somewhat inconclusive 
there is some indication that the degueesos aeration of the 
flotation concentrate pulp is Substantially fixed and 


therefore will not create a measurement problem. 


8.8 A FINAL COMMENT ON THE PULP DENSITY MEASUREMENT PROBLEM 

In the light of the results described above in section 
8.6 and subsection 8.7.3, neither temperature nor entrained 
air appear to substantially account for the poor performance 
exhibited by the pulp density measurement sub-system. Since 


electronic noise was impossible to quantify retrospectively, 
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but would be assumed to be small, the problem seems due, in 
large part, to statistical fluctuations in radiation count 
data and poor interval estimates of the slurry attributes. 

Considering the data in Table 8.4, the mean pulp 
density was calculated to be 1.0842 g/cm? with a standard 
deviation of 0.0081 g/cm’. Assuming a Gaussian distribution, 
it 1s estimated that 95% of the samples would have a pulp 
density in the range of; 1.068 < p < 1.1 g/cm*. Using the 
estimates of C (=2.8686) and D (=5,648,900) from equation 
(8.2), the expected difference in counts at these limits, 
computed from Lambert's law, would be about 
Df{exp(-1.068C)-exp(-1.1C)}=23,145, a relatively small range 
for calibration purposes. Statistical fluctations in the 
base radiation count data (i.e. p=1.08425> I1=251,912) at a 
99% confidence interval (+30) are approximately 13% of this 
range and thus would make a minor contribution to the poor 
pulp density correlation. Given the limitation on counting 
time, a stronger source would be required to reduce this 
eontrabution. 

Duplicate sampling has demonstrated that reliable point 
estimates of the attributes can be made, however, the 
interval estimate was expected to contain some error despite 
the smoothing operation. (The value of this error was not 
measured, an unfortunate oversight in the experimental 
plan.) Assuming for the moment that all of the variation in 
pulp density is due to poor attribute estimates and further 


that the relative error of estimation is the same for both 
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the % solids and % ash allows a solution for this quantity 
at 9%. (That is, if the mean attribute values had associated 
relative errors of 9%, then the expected standard deviation 
of the calculated pulp density, from the Law of Propagation 
of Errors, would be 0.0081 g/cm*.) It is expected that the 
the relative errors of estimation were less than this value 
but may have been as high as 5% and thus poor interval 
attribute estimation can make a significant contribution to 
masking the expected pulp density correlation. Furthermore, 
errors in the estimation of attribute values at extreme 
points (most distant from mean) can have a 
disproportionately large impact on on the calibration by the 
nature of the parameterization method (linear regression). 
Referring to Figures 8.2, 8.3 and 8.4 it is evident that the 
smoothing operation tends to minimize these fluctuations. 
The failure of the pulp density measurement sub-system 
to provide reliable results for calibration purposes is a 
function Of: all Of (ther contriabuting avatars previously 
mentioned. It is the opinion of the author that, of these 
factors, poor interval attribute estimation in conjunction 
with a narrow calibration range was the most significant. 
Perhaps more importantly, and with regard to OSCAA design, 
entrained air did not appear to play an important role in 


the plant trial. 
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8.9 SUMMARY 

The plant trial was judged to be a success since a 
calibration was obtained which met the performance criteria. 
' However, only one of the three process streams was monitored 
and for this stream the mineralogical composition appeared 
to remain relatively constant. 

The results indicated areas for improvement in the 
OSCAA system design (pulp density sub-system) as well as in 


experimental methods (attribute estimation). 
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9. CONCLUSIONS AND FUTURE WORK 


9.1 CONCLUSIONS 

The conclusions which can be drawn from the work 
completed in this project are presented below. 

For western Canadian coal preparation plants, the 
Steady state optimization objective is to have each circuit 
within the plant producing a product having the same 
incremental ash value, such that the overall ash 
Specification is met for each shipment from the plant. This 
will maximize production from the coal reserves, which, in 
turn, Should maximize before tax operating profit, since the 
Operating costs associated with the additional production 
are expected to be small. The difficulty in measuring the 
incremental ash content precludes the direct implementation 
of this philosophy. In general, the operator selects circuit 
product ash set points which (hopefully) correspond to the 
optimal operating scenario. Short term fluctuations around 
these set points result in production losses. Experimental 
evidence from McIntyre Mines Ltd. and communication from 
other operators indicates that these fluctuations do indeed 
occur. Moreover, they occur in a manner which should make 
them amenable to simple regulatory feedback control systems. 
Spectral analysis of experimental data can be used to help 
the operator assess the technical and economic potential of 
Such simple control systems. Better control of both 


intrinsic and especially extrinsic disturbances will reduce 
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direct and indirect production losses and should allow 
closer adherence to the steady state optimization objective 
in the longer term. 

Thusfar the sensors required for the development of 
these process control systems have received relatively 
little attention. In particular, ash sensors have only 
recently been commissioned for on-line process monitoring of 
essentially dry coal solids. Development projects for the 
on-stream ash analysis of fine coal slurries have recently 
been completed in Australia (ASHSCAN) and in Canada (OSCAA). 
The introduction of these sensors into the plant should 
provide the operator with an useful tool both in improving 
his process. understanding as well as in forming the basis 
for the development of process control schemes. 

The dabor ater calibration of the OSCAA prototype was a 
qualified success. A set of four orthogonally designed 
calibration experiments were conducted on simulated 
flotation product flow streams for the Me pityae flotation 
circuit, including one in which the) iron concentration of 
the mineral matter was independently controlled. A 
semi-empirical calibration model, derived from the theory, 
was used for data reduction. In all but one case, the 
performance criteria for range error and for relative error 
in attribute (% solids, % ash) prediction were met. The one 
exception was the flotation tailings stream where a 
combination of low accumulated density counts from the pulp 


density system (high relative error) coupled with a low 
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solids concentration in the slurry prevented successful 
calibration. It is felt that some relatively simple 
instrument design modifications would correct this problem. 
The performance of OSCAA in the laboratory compared very 
favourably with ASHSCAN in roughly equivalent service. 

The plant calibration exercise was also a qualified 
Success. An acceptable calibration was obtained for the 
flotation concentrate flow stream, even though the pulp 
density sub-system failed to operate properly. Entrained air 
did not appear to be a major source of calibration error in 
this work. On the basis of the laboratory and plant testwork 
which has been completed, OSCAA seems to be an acceptable 
instrument design for application in process monitoring and 
process control of fine coal preparation circuits. 

The OSCAA instrument design is superior to that of its 
chief competetor, ASHSCAN, because it is capable of 
accommodating changes in the mineral matter mineralogy 
(specifically iron concentration) without requiring 
recalibration. Moreover, it iS compatible with existing 
on-stream XRF equipment used in the base metal industry and 
therefore could be supported by vendors currently marketing 
in this area. From the author's discussions with 
representatives of both the coal industry as well the XRF 
vendor companies it is clear that there is considerable 


interest on both sides in continuing OSCAA development. 
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9.2 FUTURE WORK 

The work envisaged for the immediate future of this 
project may be broken into two stages, second generation 
equipment design and evaluation strategy. 

Second Generation Equipment Design 

It iS recommended that the OSCAA prototype design be 
modified prior to further evaluation testwork. The minimum 
set of modifications is given below, however, other 
modifications may also be included (e.g. computer hardware, 
dual isotope - dual detector ash analysis sub-system,...). 

a) include a passive air eliminator in the secondary 
sampling device. 

b) provide for flow stream multiplexing or multi-stream 
measurements. 

c) modify both measurememt sub-system source and detector 
designs to enhance statistical precision, simplify 
Spectral analysis (iron peak) and reduce the analysis 
time. | 

Succesful modifications in these areas should appreciably 

improve system performance and render it suitable for a more 

rigorous evaluation in the plant environment. 

Evaluation Strategy 

To more fully evaluate the second generation of 
instrument design it should be tested for a fairly long 
period of time in the operating plant environment, in what 
might be imagined as normal service. One of western Canada's 


more complex fine coal preparation circuits would be 
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preferred since this should adequately identify multiplexing 
problems, particle size effects, etc..., while giving 
industry a fair indication of instrument utility. For the 
purposes of this evaluation simple process monitoring should 
be undertaken with efforts in process control being reserved 
for late in this evaluation program or for some subsequent 
evaluation program. 

It is the opinion of the author that there are two 
reasonable means of reaching the end of this stage. The 
first would be to encourage a major manufacturer of 
on-stream XRF equipment to pursue the development using the 
information generated from this project and in-house 
expertise. Given the experience of these manufacturers and 
the performance of the first generation prototype it is 
expected that the evaluation could be completed quite 
quickly, barring political problems between the parties 
concerned. The second route would be to have some reserach 
organization, such as the Coal Mining Research Centre, carry 
out the second generation design and evaluation. This would 
probably mean a longer time to completion of these second 
generation activities, but it would allow a more equitable 
and complete transfer of the technology of on-stream ash 


analysis to the industry. 
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